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MECHANICS WASHOUT ERODIBLE FUSE PLUG 


SYNOPSIS 


Laboratory and field experiments are described that were undertaken 
demonstrate the feasibility afuse plug inthe spillwayof major From 
the data and from tractive force theory, similitude relationships are derived 
that permit the application model test results rates washout. 


INTRODUCTION 


The Idaho spillway the Oxbow project, the Snake River between Idaho 
and Oregon, will have fuse-plug control instead conventional gates. The 
prototype fuse plug designed andconstructed for controlof this spillway will 
provide economical, safe, and effective control for the maximum designed 
flow 150,000 cfs. The combined design capacity the fuse-plug spillway and 
the conventional spillway 300,000 cfs. 

The fuse plug simply rockfill dam withspecial arrangement the im- 
pervious core and strict control the rockfill and filter materials. hasa 
factor safety equal that forordinary rockfill earth dams when func- 
tions however, will washed out the reservoir water auto- 
matically suitable rate when the spillway capacity required. 

The general arrangement the fuse plug controlled Idaho spillway shown 
Fig. The fuse plug, located upper end the approach channel, 


open until October 1961. extend the closing date one month, 
written request must filed with the Executive Secretary, This paper part 
the copyrighted Journal the Hydraulics Division, Proceedings the American So- 
ciety Civil Engineers, Vol. 87, No. May, 1961. 
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WASHOUT 
concrete sill. 

the upstream end the embankment, pilot channel, wide with 
invert lower than the crest, provided start the washout the fuse 
plug predesignated location. The location this pilot channel chosen 
that the flow water, during washout, will smoothly channeled into the 
spillway chute. 

When required pass extraordinary flood exceeding the capacity 
the Oregon spillway, the fuse plug would washed out provide additional 
capacity. initiate the fuse-plug breach, the Oregon spillway gates would 
closed slightly raising the reservoir water surface the required predeter- 
mined level, one foot above the invert the pilot channel. After breaching, 
the fuse plug would wash out progressively from the upstream the down- 
stream end orderly manner. During the washout, the Oxbow reservoir 
would kept more less constant level closing the Oregon spillway 
gates, thus shifting the discharge from the Oregon spillway the Idaho spill- 
way with appreciable variation outflow. 

Fuse-plug controlled spillways have previously been designedand construc- 
ted. The Aluminum Company America has seven such spillways service 
North Carolina which were constructed under the directionof James Grow- 
don, ASCE, the originator and pioneer inthe use the fuse plug for the con- 
trol spillways the United States. Fuse-plug sections have been used 
levees for passing extraordinary floods major rivers this country and 
abroad. case unusual floods, the fuse plug designed wash out either 
automatically artifically minimize flood damages. 

most existing fuse-plug spillways have not been service long enough 
pass their design floods, thereby demonstrating their washout performance, 
literature describing the mechanics washout fuse plug insufficient 
provide practical data substantiating design assumptions. Although scattered 
information known from observations washout cofferdams and levees, 
would difficult arrive any reliable statistical deductions the 
washout characteristics. project large asthe Oxbow project the Snake 
River required more thorough investigation. 

Before the fuse-plug spillway could approved for construction, several 
fundamental engineering aspects had thoroughly investigated, sothat there 
could doubt that the fuse plug would function designed. 


How the pilot channel breach? What the water depth over the 
pilot channel before breaches 

Whatis the nature the washout process after the initial breach? What 
the speed washout? 

How complete the washout? Will any portion the plug remain in- 
tact? 


The fuse plug was designed according theoretical computations based 
mainly upon bed load transportation principles and modern concepts soil 
mechanics. However, the assumptions made the design engineers required 
verification large scale model, Furthermore, some basic information, 
which could only found through large scale model test, was needed fi- 
nalize the design. After considerable planning, was decided conduct 
1/2-scale model test the project site where the same materials could used 
for both model and prototype, the construction equipment personnel 
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would readily available, and alarge quantity water andconstant 
provided the Oxbow reservoir upstream from the diversion cofferdam. 

Because the large expenditure involved the construction and testing 
the 1/2-scale model, well interference with upstream power generation 
Brownlee Dam during the test, was impractical repeat the 1/2-scale 
model test. other words, the information which was desired from the large 
scale model had obtained one test. Therefore, preparatory tests 
small scale models the laboratory, and large scale flume models the 
field were deemed necessary. 

From experience with river models which usually are small scale and 
with small particle sizes, known that model results are, best, only 
indicative. However, for large particle sizes and large scale models, model 
test results are surprisingly accurate, demonstrated the Fort Randall 
closure. 

has also been established experimentally that the bed-load transport 
non-cohesive material hasa definite order 
avoid the troublesome problem simulating the particle size, and 
compaction, was decided touse the same materialand the same construction 
control for both the large scale field model and the prototype and establish 
the similitude relationship using the tractive force, the only predominant 
variable, the basis. was also considered important have scale ratio 
largeas possible, minimize not only the scale effect, but also the influence 
critical tractive force. 


LABORATORY INVESTIGATIONS 


Laboratory investigations, including experimental tests and theoretical der- 
ivations, were performed determine (a) the mechanics the washout pro- 
cess, and (b) the similitude relationship for the rate washout. 

For convenience, the model tests were made the downstream portion 
existing model. concrete pad and concrete blocks were placed form 
non-erodible base and vertical wall was constructed on-the right (looking 
downstream) bank. 

The fuse plug models were constructed with axes angled 64° from the face 
the right wall. The typical cross sections the models are shown Fig. 
These are for 1:20 and 1:40 scale models the proposed Oxbow fuse plug. 
The mechanical analysis the materials usedin each model are shown Fig. 

attempt was made duplicate the density other characteristics the 
materials which would used the prototype. However, considerable effort 
was made constructing the different 1:20 and 1:40 scale models achieve 
identical placement and compaction like materials. 

all, ten models were built and washed out shown Table The nor- 
mal core and thin core, referred previously, represent the two types 
cores used, the thincore section, the rockfill material was increased its 
maximum possible the purpose testing its influence washout time. 
The dimensions these cores are given Fig. 


“Dredge Fill Closure Missouri River Fort Lorenz Straub, Pro- 
ceedings, Minnesota Hydr. Convention, September, 
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the 1:20 model, the invert the pilot channel was in. lower than the 
crest the fuse plug, and the width (parallel tothe axis plug model) 
was in. the 1:40 model these dimensions were and in., respectively. 
The operation the tests was follows: 


The forebay was filled the level the invert the pilot channel. 
The forebay water level was then raised rapidly the prescribed amount 
(approximately in.) achieve initial breach action. 
The forebay level was held more less constant this elevation dur- 
ing the washout process continually increasing the inflow. 
record was kept the time which the top and bottom rows 
beads spaced 1-ft intervals along the model were washed away. 
Photographs were taken and the time recorded for each 
movie was taken Test Nos. and 


FLOW 


Embankment 
Material 


Dimensions Model Inches 


FIG, 2,—DIMENSIONS LABORATORY MODELS 


Test Results. 
Description Washout Process.—The washout process observed these 

tests consists essentially two phases, (1) the breaching the pilot chan- 
nel, and (2) the eroding embankment face. The breaching action starts with 

initial slip the downstream end the pilot channel (Fig. 4). This slip- 

page rapidly increases until the clay core exposed. The flow then cascades 

over the clay core and begins erode the material downstream nearly the 

base the model. approximately this time the clay core, which left 
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standing the removal the supporting downstream material, 
fails under the water loadand the breaching action complete. The flow then 
begins attack the base the embankment, carrying away the material 
slides into the flow and producing lateral advance the eroding face along 
the length the fuse plug. When sufficiently large portion the core ex- 
posed, breaks off. This process proceeds most regular manner grad- 
ually increasing the opening and washing out the fuse plug. 

The shape the eroding face depends the nature the materials. For 
uniformly graded material without fines, the slopes the eroding face are 
gradual and the material rolls continuously into the flow (Fig. 5), but for ma- 
terials that have some apparent cohesion mechanical interlocking, the face 
stands nearly vertical (Fig. 6). With the latter materials, the face recedes 
sudden large slips the material and sudden failure portions the clay 
core, 

Experimental Data.—The dataon the washout for eachof 
are plotted graphs Figs. and which show the eroded length the 


TABLE 


Coarse sand normal core 
Pea gravel normal core 
Crushed rock normal core 
Filter sand normal core 
Crushed rock thin core 
Pea gravel thin core 
Crushed rock normal core 
(Repeat Test 


model afunctionof time. Fig. the rate recession the eroding 
face, feet per hour; and denotes the the rate recession for the 
1:20 scale the 1:40 scale model. The average values the top and bottom 
washout times were used. Rows affected breaching action were not included. 

the interest determining laboratory test accuracy, the 1:20 crushed rock 
model test was repeated. The the two tests were remarkably close. 
The differences are well within the expected accuracy measurement this 
type test (Fig. 10). 

Analysis Experimental Results.— Analysis the laboratory washout tests 
results indicates the following: 


The washout afunctionof the grainsize the material; the larger 
the size, the slower the rate washout. However, the effect not pronounced. 

The washout rate for fuse plug between 1.05 and 1.20 times faster 
than the corresponding half-scale model. 

Increasing the volume rockfill (reducing the thickness the clay 
core) decreases the washout rate slightly. 


Theoretical Analyses.—After the pilot channel has been breached, the main 
embankment washes out regular manner, described previously. The 
following treatment the mechanics washout has been made establish 
the similitude relationship. 
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FIG, 6.—VERTICAL FACE WITH MATERIALS HAVING APPARENT COHESION 
MECHANICAL INTERLOCKING (TEST 
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Referring the definition sketch (Fig. 11) the rate recession the 
eroding face the main embankment will 


which Ris the rate recessionof the eroding face (L/T), denotes the rate 
erosionof material (L3/ T), and Arefers the cross-sectional area fuse 
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EROSION FUSE PLUG, MODEL FEET 


TIME MINUTES 
CRUSHED ROCK 


FIG, 7.—RELATIONSHIP BETWEEN 1:20 AND 1:40 MODELS 
IDENTICAL MATERIAL 


plug 


The rate erosion material can obtained from Boys equation 


which denotes the rate erosion material weight per unit time/unit 
width eroding face, isthe sediment transportation characteristic (Fig. 12) 


Thin Core, 


Thin Core 
Test No.9 


Normal core washout 
fates are shown dotted 
lines for comparison. 


Steeper slope normal 
core washout rates 
indicates inversely 
proportional the rock- 
area. 


CRUSHED ROCK 


EROSION FUSE PLUG MODEL-FEET 


TIME MINUTES 


FIG, WASHOUT RATES FOR DIFFERENT 
ROCKFILL ZONE AREAS 


for riverbed material, describes the shear tractive force the eroding 
material, refers the critical tractive force required move material, 
(Fig. 12), and the unit weight water. Now 
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which the dry unit weight material, and denotes the width the 
eroding face. 

The shear tractive force, for uniform flow open channels given 


Bottom washed out first, then 
top portion fell slices. 


Test No. 
Fine Sand 


TIME MINUTES 


FIG, 9.—WASHOUT CHARACTERISTICS FOR MATERIAL 
WITH APPARENT COHESION 


which the depth flow, and denotes the slope energy grade line (or 
water surface for uniform flow). 

For flow around the eroding face the fuse plug, however, more general 
equation required: 


which denotes the coefficient drag, and the mass density fluid. 
For uniform flow can shown that 
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10,—COMPARISON WASHOUT TESTS FOR IDENTICAL MODELS 


maximum velocity available the face, the maximum tractive 
force 


which the friction factor the Darcy equation for head loss. Using the 
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WASHOUT 
Combining Eqs. and 


R = 


which the constant proportionality inthe formula for the cross-sec- 
tional area the fuse plug 


Generally speaking, proportional that is, 


Now, proportional and turn, proportional the 1/3 power 
the relative roughness the eroding face. This last relationship follows 


FLOW 


WATER SURFACE 


PATH ERODED MATERIAL 


PLAN DURING WASHOUT 


W.S. UPSTREAM 


SECTION A-A 
FIG, SKETCH FOR THEORETICAL CONSIDERATIONS 


from the fact that Manning’s “n” proportional the 1/6 power the grain 
diameter for flows over unriffled sediment beds and from the fact that pro- 
portional follows then that 
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which denotes the mean diameter the material. Inserting Eqs. and 
into Eq. yields 


Finally, for conditions which much smaller than (that is, small non- 
cohesive material large fuse plugs) 


Eq. provides basis for determining the similitude relationship for rate 
recession between model and prototype. Let the subscripts and refer 


CRITICAL TRACTIVE FORCE, LBS./SQ.FT. 


MEAN DIAMETER MILLIMETERS 
FIG, FOR TRACTIVE FORCE THEORY 


the words “ratio,” “prototype,” and “model,” respectively. Then, for geo- 
metrically similar models using identical materials (which must relatively 
small both model and prototype) follows that 

1/3 


4 
all other ratios are unity. (Lr the scale ratio the model. 
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For two fuse plugs identical materials with scale ratio 1:2, the larger 
one will therefore erode 1.26 times faster. This relationship roughly 
verified the laboratory tests described previously for which the 1:20 model 
eroded 1.05 1.20 times faster than the corresponding 1:40 model. Although 
the experiments indicate rate erosion for the larger model that some- 
what slower than the theory indicates, should remembered that neither 
the theory nor the experiments are precise and, least froma practical stand- 
point, this agreement between theory and experiment acceptable. 

Summary.—From the laboratory tests and theoretical analyses the me- 
chanics washout fuse plugs, the following has been found: 


The mechanics washout can explained the basis the laws for 
sediment transport. 

The rate recession the eroding face varies approximately with the 
1/3 power the scale ratio for fuse plugs constructed with identical materials 
(non-cohesive) and methods. 


LARGE SCALE FIELD MODEL TESTS 


General tests small scale models inthe laboratory 
had supplied the preparatory information necessary correlate the results 
the large scale model test prototype performance. order finalize 
the design the Oxbow Idaho spillway, the construction the 1/2-scale field 
model was commenced June 1959 with the test scheduled for late July. 

Purposes Large Scale Field Model main purposes the large 
scale model test were verify the design assumptions, and seek additional 
information for the final design the fuse plug. 

The assumptions verified were follows: 


Mechanism the breaching action the pilot channel; 
Depth water required effect successful breach; 
Completion breach before lateral erosion; and 

Washout process for lateral erosion. 


The following additional information was needed: 


Time function the breach process; 

Time function the lateral washout process; 

Effect the core rate washout; 

Completeness washout; and 

Anything unforeseen which could cause malfunctioning the fuse plug. 


Model Scale and Similarity.—The selection a1:2 ratio for the large scale 
field model was prompted number considerations. The primary factor 
was the existence the project readily adaptable site where 1/2-scale 
model could constructed with suitable length-to-height ratio. 

The predominant material fuse plug the rockfill that occupies more 
than 75% the volume (in the design shown), has relatively much larger size 
particles than the filter, and, therefore, controls the rate washout. de- 
termined the laboratory investigations, the similitude relationships between 
the prototype and scale model withidentical materials and construction 
methods would the same between the 1:20 and 1:40 models. 
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Because the impervious core material was have the same inplace density 
and moisture content for both the model and the prototype and, thereby, the 
same unit strength, the similitude relationship the impervious core had 
determined according its failure action. design, the core ex- 
pected break cantilever load, both breaching and lateral washout. 
Therefore, the core thickness was proportioned according 


which the thickness ratio, and denotes the scale ratio. order not 
disturb the similarity the rockfill, the reduced thickness the core (more 
than the 1:2 ratio) was made the material, which little 
effect washout rate. 

For testing the initial breach action, the prototype head and dimensions 
the pilot channel were used because the materials and slopes are identical 
both model and prototype. 

Flume Model Clay Core.—After having studied the general breaching and 
washout actionon small scale models the hydraulic laboratory, was more 
less certain that the initial breach would start the downstream edge 
the pilot channel with very small water head. However, because the imper- 
vious core material, the fine filter, and even the rockfill material may behave 
differently after being compacted maximum density, was proposed test 
the breaching action wood flume the site before proceeding with the 
large scale washout model test. The specific purpose this model was fol- 
lows: 


observe closely the erosion phenomena the downstream rockfill 


and filter material; 
verify the break-up action the impervious core. 


avoid all complications, the model the flume was have the same scale 
that the large scale field model, with the same materials and compaction. 
wooden flume was constructed the project site which was able ac- 
commodate high crest section the pilot channel portion the 1/2- 
scale fuse plug. 
Two flume tests were performed, from which sufficient information was 
obtained for the construction the large scale model. 


Test No. 1.—The flume model had section shown Fig. 13. The im- 
pervious core was supported processed fine filter material the down- 
stream side, and covered the same material the upstream side. Flume 
floor and walls are in.-by-6 in. tongue and groove with polyethylene sheet 
lining. The mechanical analyses are shown Fig. 14. 

The model was constructed obtain the same compaction the core 
and filter zones the prototype. After turning the water and removing 
the downstream bulkhead, the downstream processed fine filter material ero- 
ded only the extent its contact with the lower nappe falling over the im- 
pervious core. Due the apparent cohesion the compacted, well-graded 
material, the processed fine filter material stood almost vertical for the en- 
tire height the model. 

After mechanically removing the processed fine filter material, the im- 
pervious core was tested determine the point failure under cantilever 
action. (The bond between the core and the walls the flume was carefully 
broken the top portion the model.) The first break the core occurred 
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about in. from the top with less than foot water over the top. The sec- 
ond break occurred about 4.5 from the top. 


Test No. second flume model had section shown Fig. 13. 
The impervious core was identical tothat Test No. but the processed fine 


FLUME 


Water supply pipe 


Temporary bulkhead for construction 
removed before test 


SECTION FLUME MODEL 


Temporary for construction LEGEND 
removed before test 


Processed fine filter 
LO’, 1.75' 08' 106.5 Processed concrete sand 
we 


SECTION FLUME MODEL NO.2 
FIG, 13,—LARGE SCALE FLUME MODEL TEST 


filter material was replaced with downstream rockfill zone concrete ag- 
gregate and filter zone composed thin layer concrete sand over pit- 
gravel. were moreor less uniform materials. Their 
gradations are shown Fig. 14. 
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WASHOUT 
The materials were placed the same manner for Test No. but with 
provisions prevent bond and wedge action between the impervious core and 
the walls the flume. 

After the end bulkhead was removedand started, the material down- 
stream the core was gradually washed out from the top. the water hit 
the downstream rockfill material and washed out the flume, the rockfill 
and filter material beneath the nappe gradually receded surface sliding. The 
first break the unsupported core occurred about in. from the top. The 
washout cycle was then repeated. The second break the core occurred 
height about 2.5 below the first. 


the basis the flume model tests, was concluded that the impervious 
core the large scale model was satisfactory with respect thickness and 

These tests also indicated that the compacted, well-graded riverbed sand 
and gravel, because its apparent cohesion, would not satisfactory for use 
the pilot channel section. was demonstrated that clean, uniform material 
was necessary for the downstream filter and rockfill effect rapid and suc- 
cessful breach. 

Half Scale Field Model Tests. 

Description Half Scale Model.—The location chosen for the construction 
the model was natural rock saddle the left abutment the Oxbow cof- 
ferdam. Here large water supply was available from the natural river flow 
and the Oxbow reservoir above the cofferdam could used regulating 
pool keep the water surface almost constant throughout the test. The pool 
between the cofferdam and the main dam could store all the water passing the 
model without any flow over the Oregon spillway. 

The first step preparing the model site was remove part the exist- 
ing cofferdam, and expose the rock foundation. Because the rock surface was 
quite irregular, concrete sill was poured for the cutoff under the impervious 
zone. The rock onthe right side sill sloped upto form 
abutment. However, the left side, the rock was lower and very irregular, 
vertical steel sheet pile wall was constructed form the abutment for the 
model. The piling was set excavated rock and backfilled with 
concrete. The top the piling was securely anchored logs buried well back 
the Hand-tamped clay was placed back the piling. 

The pilot channel the model was adjacent the sheet pile abutment. 
the prototype, the pilot channel will adjacent aconcrete-faced rock 
ment, excavated steeply possible. difference breaching action 
foreseen due the probable minor difference abutment slopes. 

compacted clay foundation pad, thick, was placed under the en- 
tire model the level the concrete sill. 

The length the base the fuse-plug model, measured along the concrete 
sill, was 51.6 ft. The crest length was ft, and the height was 13.5 mini- 
mum, except the pilot channel. The crest was 1-ft elevation the 
right abutment thanit was adjacent the pilot channel. The height the pilot 
channel was 11.5 ft. Typical sections the model constructed are shown 
Fig. 15, and picture the complete model shown Fig. 16. 

the lower the model, the rockfill materials were dumped directly 
place trucks, and spread bulldozer. For the remainder the model, 
the materials were stockpiled near the site, transferred the embankment 
clamshell required, and spread hand. 
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Compaction the riverbed sand and gravel, and the processed fine filter, 
was accomplished vibrating roller. Liberal quantities water were ap- 
plied tothese materials. The concrete aggregates and fine gravel were heavily 

The clay core was spread 3-in. lifts, sprinkled required, and com- 
pacted pneumatic hand tampers. order get satisfactory bond between 
lifts was necessary roughen the contact surface and moisten before 
placing the next lift. general, compaction was between 85% and 90% peak 
modified American Association State Highway Officials (AASHO) density, 
and the moisture content was within 1-1/2% optimum. 

Material Selection and Laboratory Tests.—The model was constructed 
the same materials economically available for usein the prototype. The model 
was designed acomposite section, using several rockfill materials, order 
that washout rates for the various materials could established and thereby 
ensure selection the most advantageous materials for the prototype. 


FIG, 16.—SANDBAGS JUST CLEAR PILOT CHANNEL 


The flume model tests had indicated that pit-run riverbed sand and gravel, 
when placed wet and compacted, had considerable apparent cohesion. Because 
this property might slow the erosion process, was decided use this ma- 
terial where would the last wash out; that is, the end opposite the 
pilot channel. 

the pilot channel was necessary have easily moved, cohesionless 
material. Concrete aggregate passing 3/4-in. screen and retained No. 
screen was sure meet these requirements and was therefore used. Be- 
was anexcess the larger concrete aggregates available, 
was decided test section using rock between the 3/4-in. and 1-1/2-in. 
sizes. This section, long, was placed adjacent the pilot channel sec- 
tion. 
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The processed fine filter used the main dam showed the same property 
apparent cohesion the riverbed sand and gravel. Therefore, this filter 
material was restricted use only inthose sections using the sand and grav- 
el. For the sections using concrete aggregates, pit-run fine gravel witha 
minimum thickness concrete sand onthe face contact with the impervious 
core, was used for the filter. The fine gravel showed cohesion under any 
conditions. The impervious core was constructed the same material the 
main dam impervious core. The upstream and downstream slopes were faced 
with the 1-1/2-in. concrete aggregate. 

For the results standard laboratory tests performed the various em- 
bankment materials, see Table and mechanical analysis curves, Fig. 14. 

addition tothe standard laboratory tests referred previously, series 
special flexural tests were performed thecore material. The tests were 


TABLE TESTS FUSE PLUG MATERIALS 


Material 


(1) 


Clay core 


Fine Filter 


Riverbed Sand and Gravel 
Processed Concrete Sand 
Pit-run Fine Gravel 
Concrete Aggregate 
Concrete Aggregate 
3/4 in, 1-1/2 


Maximum density determined by: (1) Modified AASHO method for clay, (2) field 
density test after compaction for fine filter, and (3) dry, rodded method for all other 
materials, 


required quantitatively evaluate the maximum cantilever projection the 
core that could occur the pilot channel section during removal the down- 
stream supporting rockfill. This value, turn, permitted rough computa- 
tion the time required start initial breakup the core, and the effective 
breach the section. Flexural strength values were desirable for both the 
as-placed material and the saturated material get acomparison between the 
fuse-plug model, which would betested with the core still inthe as-placed con- 
dition, and the prototype, which would certainly subjected prolonged 
period submergence tothe fullreservoir levelat El. 1805, before would 
required wash out. 

The test specimens were prepared compacting the soil approximately 
the same density, and about the same moisture content, the material 


foot meters per 
per square 
second foot 
(2) (3) (4) (5) (6) (7) 
min 
ave 
j 
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placed the model. The specimens were in. square section, and in. 
long. The soil was placed and compacted about five layers. 

The test consisted supporting the specimen its ends, and applying 
uniformly increasing load the third points until rupture occurred. The load 
was appliednormal the directionof compaction, effects lam- 
inations the test specimen. 

The test results showed modulus rupture between 3000 psf and 4000 psf 
for the material placed. conclusive results were obtained the effects 
saturation. However, the tensile strength was apparently significantly re- 
duced due prolonged submergence. 

Measuring Devices. 

Embedded Washout Rate Indicators:—To get informationon the rate which 
the interior portions the fuse plug were being eroded, series washout 
rate indicators were installed during the course construction. These were 
simply electrical switches which were held closed the earth pressure with- 
the embankment. the erosion progressed the earth pressure would 
relieved, the circuit broken, and would out the observer’s station. 

The indicators were installed vertical plane along the center the 
crest three vertical lines shown Fig. 17. general, they were in- 
stalled 2-ft increments elevation. The wires from each indicator went 
vertically down the clay foundation, upstream the concrete sill and then 
out over the Idaho abutment. This routing prevented damage the wiring be- 
fore operation the respective indicators. 

Grid Lines Downstream Face:—A 5-ft grid was paintedon the downstream 
face the fuse plug model, shown Fig. 16. Three independent observers 
recorded the elapsed time between the start the test and the removal the 
grid intersections, aid establishing the washout rates for the various 
materials. 

For the final determination washout rates, careful study the motion 
picture, taken the time the large scale model washout, was made. The 
grid lines were the dimensional control, and repeated stop watch measure- 
ments established the time. 

Description Washout. 

Provisions for Starting Test:—A 10-ft long section the crest was left 
lower elevation than the remainder form the pilot channel, through which 
the water would flow and start the breaching action. Close control the res- 
ervoir water surface elevation was required, but difficult achieve. order 
get the proper head for the test, without danger premature washout, sand- 
bags were employed the pilot channel. 

Four layers sandbags were used close the pilot channel the model 
crest El. 1748.8. The sandbags were arranged wire rope sling at- 
tached After the reservoir stabilized the desired eleva- 
tion, the sandbags could lifted quickly start the test. 

Final Test:—In the morning July 22, 1959, all preparations for the final 
test had been completed. The reservoir water surface was brought the 
proper level, one foot above the invert the pilot channel. sec after 
11:11 a.m. the initial breach began with the lifting the sandbags clear the 
pilot channel (time consumed clearing sand bags, sec), shown Fig. 
16. 

soonas were lifted, the water gushed through the pilot chan- 
nel and, flowing down the slope, commenced the erosion downstream rock- 
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fill starting the top the slope. Within sec after the start the test the 
erosion had progressed tosuch extent that the upper portion the core was 
acting weir, with the nappe clear the material below, Fig. (time 
The first break the core occurred min and sec after the start 
the test. Fig. (time 11:13+) shows the moment which the second break 
the core had just occurred. The time, calibrated from the motion picture 
record min and sec after the start the test. The breach action was 
considered completed this point further change occurred surface 
profile the water rushing through the gap the color the water. 


Existing 
Cofferdam 
EMBEDDED WASH-OUT INDICATOR EL. 
Line Line Stee! Sheet 


Piling Abutment 


BASE EL 1735.3 


DOWNSTREAM ELEVATION MODEL 
OBSERVED WASHOUT TIME 


The numbers on the signs along the crest designate 10 foot 
stations from the Sheet Pile Abutment. 

Time is in minutes and seconds from stort of wash-out. 

wash-out grid intersections. 

wash-out embedded wash-out indicators 


Ist break of core 107 
2nd breok of core 
(completion of breach) 


Approximate limit for 
concrete aggregate section 


Erosion rate- 4 ft/min Transition Erosion rate 5.6f1./min Breaching 
(Time 1Om 30s. avg) Time 2m 3is. Time 3m. 338. Time 
im. 276. 


DOWNSTREAM ELEVATION MODEL 
FIG, 17,—RECORD WASHOUT RATE 


The lateral erosion the model began during the initial breach period 
the rockfill material along the side slidinto the water rushing through the gap 
and the core was left overhanging along the side. the uniformly-graded 
rockfill and filter slid out from underneath, the core was left cantilevered and 
sections would break off intervals. The surface the uniformly-graded 
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FIG, MOMENT SECOND BREAK CORE, 
BREACH ACTION COMPLETE 


FIG, 18.—OVERHANGING FIRST BREAK CORE IMMINENT 
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material sliding intothe water assumed slope slightly the angle 
repose the material. This phase the washout min. 
Fig. (time 11:17) shows the characteristic slope the uniformly-graded 
rockfill and filter under full face lateral erosion. The change fromthe uniform 
the well-graded material occurs Station and the vertical break can 
seen the downstream face above El. 1795, the upper white line. 

The washout the well-graded material characteristically different from 
that the uniform material. Due apparent cohesion and small amount 
fine material the well-graded riverbed sand and gravel, the end the par- 
tially washedout fuse plug model assumed vertical attitude. The ero- 
sion this portion the fuse plug model proceeded the churning water 
undermined the vertical face. intervals the undermined portions the down- 
stream rockfill broke off and fell into the water and the material was carried 
away. the core became exposed and was left unsupported, sections would 
break off. Fig. (time 11:22) shows the vertical face the well-graded ma- 
terial and the churning, undermining action the water. 

soon the test was over, the spillway gates Brownlee were closed, 
and the powerhouse release was orderedcut down much possible that 
the examined early possible. 

12:30 the powerhouse release was reduced 6000 cfs, and kept 
there until 2:30 p.m. The Oxbow reservoir receded rapidly. 3:15 p.m. the 
water surface downstream the fuse plug was the same elevation the 
Oxbow reservoir, and the water started toflow into the reservoir. this time 
the reservoir water surface, well the tailwater pool, was El. 1738.5, 
below the temporary crest the Oregon spillway. 

Analysis Test Results.—The large scale field model test has clearly 
demonstrated the details the washout action. Information which had pre- 
viously been determined from the laboratory and flume model tests was con- 
firmed, and some new information was revealed which could not found from 
the other model tests due their limitation scale and size. 

The continuous motion picture been found most valuable, because 
enabled the engineer study the washout action closely.and repeatedly. 

After the Oxbow reservoir and the pond downstream the cofferdam had 
receded, the model foundation was examined. 

All information obtained, including the river gage records Oxbow and Hells 
Canyon, have been carefully analyzed, and the results are follows: 


Breach Action:—The breach action was entirely expected. showed 
close similarity tothe small scale model. The core broke the 
flume model, first, small break, then large break, and then the 
breach action all the way the base. There was significant seepage into 
the rockfill the downstream edge the pilot channel before the initial breach. 

Lateral Washout:—The uniform material washout process was exactly 
estimated. The kinetic angle repose was about 45° for the uniform material. 
For the material with apparent cohesion and small amount fines it, the 
washout face was vertical, and the washout process was three-step action; 
undermining the base the vertical face, vertical slide the filter 
material, and the breakof the impervious core. The embedded indicator record 
also indicated vertical face erosion under water. 

Rate direct observations (three observers) and the 
motion picture record, the washout time function shown Fig. 13, and the 
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FIG, FACE WELL-GRADED MATERIAL 
(COMPARE WITH FIG, 


20.—SLOPE UNIFORM MATERIAL (COMPARE WITH FIG. 


results summarized follows: 


Breach (Complete) min and sec 
Rate washout, uniform material 5.6 per min 
Rate washout, well-graded material per min 


Hold-up time the core was insignificant. 
Order and Completeness:—The entire washout proceededin anorderly man- 


ner. Breach was completed before lateral erosion. There was evidence 
any portion the plug remaining beneath the water during the washout and 
Fig. shows the completeness the washout after the foundation was ex- 
posed next morning. The reservoir water surface the start test was El. 
1747.80. The reservoir water surface the end test was El. 1747.98. The 


maximum water surface during test was El. 1747.98. 
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FIG, 22,—LOOKING TOWARD LEFT ABUTMENT STEEL SHEET PILING 


Headwater and Tailwater:—From reservoir water gage records, was found 
that during the entire test, the headwater varied only 0.18 ft. The breach head 
was ft. The slope the ground downstream the model was approximately 
10%. There was settling the material the escape channel. The dike 
downstream the model was constructed study its effect the washout, 
but noticeable effect was found. 

Summary.—The large scale field model test has conclusively demonstrated 
the following: 


The structure showed the ability withstand water load with reser- 
voir surface above the pilot channel invert for 4.5 hr, and also rapid 
drawdown actually more severe than expected for the prototype. 

The breaching action was initiated with water head one foot above the 
pilot channel, and completed the base the model less than two minutes. 


WASHOUT 


uniform, cohesionless material essential for the downstream filter and 
rockfill zones the pilot channel achieve rapid and successful breach. 

The lateral washout process for thetwo different types materials was 
clearly demonstrated. The rate washout for the uniform, cohesionless ma- 
terial was found 5.6 per minute, and for the well-graded material with 
apparent cohesion 1.4 per min. Through proper proportioning the lengths 
involving the two kinds materials, suitable value for the total washout time 
the prototype fuse plug can achieved. 

Theimpervious core design satisfactory, because the “hold-up” action 
the core was found insignificant. 

The model fuse plug was completely washed out the concrete founda- 
tion sill orderly fashion, without unforeseen developments. 


CONCLUSIONS 


The small scale washout tests the laboratory have established definite 
pattern similarity between 1:20 and 1:40 scale model, composed the 
same materials manner. The 1/2-scale field model 
test has also shown similar washout mechanics. 

The ratio washout rates from the small scale model tests, for identical 
materials and construction, and scale ratio 1:2, checks very closely with 
the theoretically derived value. Therefore, can logically concluded that 
the actual washout the prototype can closely approximated apply- 
ing this ratio washout rates the results fromthe 1/2-scale field model test. 

The design assumptions the fuse plug have been verified, all respects, 
the actual large scale field test. Based upon the performance this test, 
the prototype fuse plug design, including the selection materials, will en- 
tirely satisfactory effecting successful breach, and complete washout, 
orderly manner, planned. 

using two types filter materials, considered definitely 
feasible and practical construct fuse plug for suitable washout rates along 
different portions and appropriate total duration the washout. 
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THROUGH LAYERED ANISOTROPIC POROUS 


David ASCE, and Jacob Bear2 


SYNOPSIS 


The paper summarizes investigation seepage from leveed rivers into 
low-lying adjoining agricultural lands. Only the flow rate and distribution 
seepage function the surface and subsurface boundary conditions were 
studied. Flow nets idealized cross sections were determined electric 
analogy model. Boundaries represented channel bed and bank, levee base, 
lower impermeable stratum, and the water table adjoining drained agricul- 
tural land. Variables studied included two subsurface layers different per- 
meabilities and anisotropies, layer arrangements and thicknesses, channel depth 
and width, levee base width, and water table slope. 


INTRODUCTION 


Seepage from rivers and natural channels serious problem many lo- 
cations. Significant contributions the literature this and other types 
seepage have been made several engineers (1), (2), (4), (5), (6), (11), 
(15), (17), (21), (26).3 When channel exceed the adjacent ground sur- 
face elevation, water moves through and under the confining levees into adja- 
cent lands. drainage facilities are inadequate. the soil becomes saturated, 


Note.—Discussion open until October 1961. extend the closing date one month, 
written request must filed with the Executive Secretary, ASCE. This paper part 
the copyrighted Journal the Hydraulics Division, Proceedings the American So- 
ciety Society Civil Engineers, Vol. 87, No. May, 1961. 

Presented the March 1960 ASCE Convention New Orleans, 

Assoc, Prof. Civ. Engrg., Univ. California, Berkeley, Calif. 

Formerly Graduate Research Engr., Univ. California, Berkeley, Calif.; present- 
Lecturer, Technion, Israel Inst. Tech., Haifa, Israel. 

Numerals parentheses refer corresponding items the Appendix Bibliography. 
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water ponds the surface, and crops are damaged. The Sacramento Valley 
California typical example area which suffers extensive seepage 
damages. Losses the period 1937-1953 (at 1953 prices) were reported 
$1,834,000 per year, $102 per acre affected (20). The problem seepage 
this area prompted the present investigation. 

The purpose the investigation was study the physical picture flow 
water from stream channels through subsurface strata the adjacent ground 
surface. Natural conditions governing these flows include geology and soils, 
topography, ground water, rainfall, runoff, and channel stage. The range 
conditions for each these variables, together with their possible combina- 
tions, makes apparent that not all situations could studied. Instead, key 
variables were selected for study that results could applied not only 
the Sacramento Valley, but also other areas with similar conditions. The 
present work deals, therefore, with general aspects seepage from rivers. 

Most seepage studies the Sacramento Valley have been concerned with 
field conditions the affected areas, periods and durations flooding, and 
economic losses. Comprehensive surveys are summarized reports the 
Joint Committee Water Problems the California Legislature (9) and 
Plumb, ASCE, and others (20). Because the present study was lab- 
oratory investigation seepage flow nets for steady flows, these reports and 
others (10), (16) furnished valuable information relating construction 
idealized cross sections. Results from electric analogy model were ana- 
lyzed and presented graphical form facilitate interpretation and applica- 
tion various field conditions. view the number and complexity cross 
sections investigated for quantities and distribution seepage, model 
investigations rather than analytic studies were necessary (application ana- 
lytic techniques two simplified special cases shown Appendix I). Hy- 
draulic models have often servedas valuable tool supplement analytic ap- 
proaches (24, pp. 307-321). The electrical analogy was most suitable for the 
present investigation. Analysis thistype applicationto seep- 
age studies have been presented Pavlovsky (19), Vreenden- 
burgh and Stevens (22), (27), (28) and others (8), (12). recent sum- 
mary the theory the analogy has been prepared Malavard (14). 

most cases, subsurface strata through which seepage flows occur are 
nonhomogeneous and anisotropic. main emphasis the present 
investigation was seepage through layered anisotropic soils. 


SEEPAGE CROSS SECTION 


Conditions Sacramento River servedas for selecting cross 
sections toinvestigate. The variety, wide diversity layered 
alluvial deposits along the river precluded the possiblity determining seep- 
age through actual geologic sections without having extensive field data avail- 
able. Instead, about fifty representative hypothetical cross sections were con- 
structed for model tests (25). 
Almost all seepage occurs the zone above the highest more less con- 
tinuous clay layer, which acts lower impermeable boundary. the per- 
meable zone various layers deposits can found. Based types geo- 
logic deposits, layers could reduced those having relative permeabilities 
50, and These layers are essentially horizontal, but may vary inthick- 
ness and their relative vertical positions. most seepage occurs within 
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less than mile the river, the model was limited maximum distance 
3,000 inland fromthe levee. maximum depthof 100 from the water table 
the impermeable clay layer was also selected. symmetry assumed 
about the center the river, section only one side the river need 
investigated. 

Other assumptions made for the model investigation were that the river 
channel rectangular with depth few cases other depths were 
also studied) and width 600 ft. Levees were assumed situated ad- 
joining the river channel, impermeable, and have base width 100 
the elevation the inland water table. impermeable levee, course, 
implies that all seepage passes under the levee. Depending upon the mode 
construction and materials employed, seepage through the levee may may 
not significant item. 

final variable necessary evaluating seepage that the anisotropy 
permeable zones. Most alluvial strata are deposited horizontally with pro- 
nounced tendency for agreater uniformity grain size and depositional struc- 
ture horizontally than vertically. Commonly this results condition ani- 
sotropic permeability where the horizontal permeability may greatly exceed 
that the vertical. Ratios vertical horizontal permeability 0.1, and 
0.01 were selected give range representative conditions. 

The position the water table river stage and ground surface 
determines relative seepage rates and the length time before seepage will 
appear ground surface. Seepage deep from short duration 
flood stage may not affect surface conditions, whereas location with high 
water table may have marked surface effect. 

most irrigated agricultural land the water table controlled within 
narrow limits subsurface drainage facilities. The level commonly main- 
tained below ground surface. For steady flow conditions studied 
the electric analogy model, the water table was assumed stationary. 
This might represent levels controlled either subsurface drains pond- 
ing ground surface. Both horizontal and sloping water tables were investi- 
gated. 

Field (9) and laboratory (23) investigations have confirmed that seepage 
from flood flows proportional duration-height factor. This study was 
limited the steady case with emphasis the distribution seepage inland 
from the river. Factors river stage magnitude and duration are variables 
which are necessary applying these results obtain quantitative seepage 
for given situation and location. 

Based these considerations the general cross section studied 
shown Fig. From Fig. the channel depth, has range values 
100 ft; the permeable zone depth, constant 100 ft; thickness 
upper layer, ranges from 100 ft; thickness lower layer, has 
range ft; length cross section from centerline river, has 
range 1,000 3,300 ft; range assigned for the river stage above 
levee base, the range channel half-width, 100 300 ft; the levee 
base width, constant 100 ft; distance landward from toe levee, 
ranges from 600 2,900 ft; ranges are assigned for the permeability 
the upper layer, Ky, the lower layer, the ranges for both the anisotropy 
the upper layer, and the lower layer, are 1/1 1/100; 
and the ranges for the ratio permeabilities the two layers, 

The channel depth and the channel half-width were constant for all but 
few model tests. Significant variables are the thicknesses the upper and 
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lower layers, and respectively. These are related 


IMPERMEABLE RIVER 


GROUND WATER TABLE LEVEE 
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FIG, 1.—DEFINITION SKETCH CROSS SECTION ASSUMED 
FOR MODEL ANALYSIS 


ous combinations the significant variables (25). 


THEORY THE ELECTRIC ANALOGY 


The electric analogy well-known device for studying the field flow 
fluid through porous media (8), (18). The analogy based the similarity 
between the differential equations which describe the flow fluid through por- 
ous media, and those which govern the flow electric current through conduc- 
ting materials (14). Methods for treating anisotropic layered media trans- 
formation into equivalent isotropic fields scale distortion have been des- 


cribed Maasland (13). 


The lengthof cross section was varied required obtain practical ver- 
tical scale define the seepage flow pattern each model test. Although the 
river stage determines the rate subsurface flow, values were assigned 
since has effect the distribution seepage. Similarly, the permea- 
bilities the layers, and were unassigned, but their relative values 
expressed their ratio form key variable the investigation. The 


anisotropies the individual layers are expressed the ratio their re- 
spective directional permeabilities, and Based the 
variables defined Fig. and study field conditions seepage areas 
the Sacramento Valley, series model tests was outlined covering vari- 
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Model scales and dimensions were based 


which the subscript denotes model-prototype ratio. selecting 
the common boundary between layers was made coincide after transforma- 


tion. illustrative example typical cross section showing prototype and 
model dimensions shown Fig. 


DESCRIPTION THE MODEL 


The electric analogy model for this study consisted several components, 
including the electrolytic tank, the electrical circuit, the conducting medium, 
the analogy boundaries, and the drawing table. photograph the assembly 


the conducting medium was electrolyte, large shallow lucite tank, 
with dimensions in. in. in., was constructed hold the liquid. 
The bottom the tank was plane and was levelled insure uniform depth 
electrolyte. all tests the maximum length the tank was utilized. was 
found that making test the largest possible size the model minimized 
inaccuracies and boundary errors. 

The electric circuit for the model shown Fig. 60-cycle cur- 
rent was connected the electrolytic tank. The potential divider accurately 
divided the voltage into 1,000 subdivisions that readings 0.1% were pos- 
sible. Visual identification the the probe position the tank was 
obtained with cathode ray oscilloscope connected the probe and the poten- 
tial divider (Helipot). Whenthere was potential difference between the probe 
and the desired potential set the potentiometer, current flowed through 
the oscilloscope. 

The resistance box provideda means for establishing wide variety fixed 
potentials electrodes the tank. Most tests required only the limiting po- 
tentials 100%; but represent sloping water tables, intermediate 
potentials with specified differences were necessary. The box contained forty 
precision resistors 150 milliohms each, connected series. order 
avoid adjustments when intermediate potentials were being connected the 
electrodes the tank, resistances the resistance box were made relatively 
small compared the resistance the electrolyte the tank. this way 
only very small amounts current were consumed the intermediate points 
the main line. The deviation variable current from linear distribution 
potentials along the line was negligible. 
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Various concentrations copper sulfate solution well tap water 
were employed the conducting medium for these tests. Liquid conductors 
are advantageous when large number tests with different boundaries are 
studied. Solutions could poured into the model readily, and after each 
test could drained into storage containers. Also, where ratio 


(a) PROTOTYPE 


(b) MODEL 


FIG. PROTOTYPE AND MODEL DIMENSIONS FOR 
ONE CROSS SECTION 


two layers 50:1 was under investigation, this difference was obtained 
selecting solutions witha corresponding conductivity ratio. Liquid depths over 
the entire cross section were kept equal about 2in. Withalternating current, 
copper electrodes, and copper salt the tank, problems polarization 
were encountered the electrode surfaces. 


FIG, 3,—VIEW ELECTRIC ANALOGY MODEL 


SALINITY METER 
OSCILLOSCOPE 


STANOBY 
SWITCH 


INSULATING 
BOUNDARY 


POTENTIAL 
RESISTANCE BOX OIVIDER 


150mn 


FIG, CIRCUIT FOR ELECTRIC ANALOGY MODEL 


LAMP 
FUSE a SWITCH 
PROBE 
e 20p 
| 


May, 1961 


Three kinds boundaries were required the model: permeable and im- 
permeable outer boundaries, and inner boundaries between layers different 
permeability. Permeable boundaries were made copper strips in. wide 
bent fit the configuration the boundary and carrying given potential 
means connection with the resistance box. For boundaries representing 
sloping water tables, variable potential was required. Because the pressure 
constant (atmospheric) water table, the potential everywhere pro- 
portion the elevation. Therefore, such boundary was reproduced se- 
ries closely-spaced copper electrodes (Fig. 3), each connected its cor- 
responding potential the resistance box. The potentials along the boundary 
were expressed percentages the total head difference betweenthe channel 
water level and the lowest elevation. Impermeable boundaries were 
nonconducting strips lucite, in. wide, the lucite wall the tank. 

The inner boundaries had provide waterproof seal between the two liq- 
uids, and had conduct electricity freely from one layer the other, but each 
point boundary had maintain unique potential value. The boundary was 
constructed modifying procedure first suggested Vreedenburgh and 
Stevens (22), (28). Lucite strips were notched and perforated 3/16-in. in- 
short length 1-mm diameter copper wire was inserted ineach 
hole and bent sothat the ends crossed the notch above. All boundaries were 
held place modeling clay the bottom. 

Outer boundaries were interchanged during each test. Electrodes per- 
meable boundaries enabled equipotential lines electricity and flow 
sketched, while placing electrodes impermeable boundaries allowed stream- 
lines located. Variable potential boundaries could interchanged only 
after computing the redistribution potentials along such boundaries. 

Movement the probe was recorded drawing paper 1:1 scale 
means pantograph connecting the probe and pencil. connecting points 
the same potential, each equipotential line was obtained. Sketching time for 
one set lines cross section amounted about min. 

general, with proper instrumentation, any degree accuracy can ob- 
tained electric analogy model. the present work was concerned with 
ground-water flow, basic data such boundaries, permeabilities, and aniso- 
tropy ratios were only approximately known. highdegree accuracy, there- 
fore, was not required the idealized model cross sections. Results were 
estimated accurate within 5%. 


THE FLOW NET 


The curves constant (equipotential lines) and constant 
(streamlines) form two orthogonal families curves which together form the 
flow net. Each two neighboring streamlines form flow channel, which the 
amount flow equal the difference between the two the cor-- 
responding streamlines. flow net constructed withan integer num- 
ber flow channels, each transmitting the same amount flow. Similarly, 
the number equipotential lines chosen such that the drop head between 
adjacent ones will the same throughout. the total loss head, 
denotes the number equipotential drops 4H, and defines the number 
stream channels with flow rate per unit width each, then Darcy’s 
law 
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which the streamline spacing and the equipotential spacing. The 
ratio a/b remains constant throughout any homogeneous isotropic region. 

For reasons simplicity was decided draw the flow nets for the pres- 
ent work with 20. This gave potential drop betweentwo equipoten- 
tial lines equal the total available head and flow flow chan- 
nel equal the total flow from the river. The total flow per unit dis- 
tance along the river channel was given 


and average a/b from several measured “squares” was found sufficient 
insure reliable flow value. When two layers are present, the ratio a/b 
different each zone, that 


draw equipotential lines, the river bank and ABCD Fig. 
were connected tothe 100% electrical potential, and the ground-water leveland 
the end the model EFGH Fig. 2(b)] were connected the electric 
potential. The impervious base the levee impervious bottom (HI), 
and the right edge the model which line symmetry (IA) were made 
non-conductive lucite. From theoretical standpoint the line should 
have been infinitely distant from the channel, but the error introduced as- 
suming about 1,000 2,000 from the channel negligible. For 
these boundaries and given ratio, the equipotentials were drawn. 

For drawing streamlines, impermeable and fixed potential boundaries were 
interchanged. the sametimethe conductivities layers were inter- 
changed. 

The flow net drawn directly from the isotropic model orthogonal one. 
The prototype itself anisotropic; therefore, the actual flow net not or- 
thogonal. When the drawing reduced distorted scales, that is, the hori- 
zontal reduction differs from the vertical, the resulting drawing also non- 
orthogonal flow net (see Figs. and 6). 

Almost all tests had upper boundary formed fixed horizontal water 
table, the elevation which was taken the datum level 0). This ap- 
proximated conditions existing fields adjacent river where drains con- 
trol the water table. this way the problem unknown phreatic surface 
was avoided; phreatic surface can studied electric analogy only 
trial and error. However, some fieldcases interest the controlled steady 
water table not horizontal, but slopes away from the river witha gradient 
from about 1:500 1:1500. Two tests were performed with sloping water 
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tables study these situations. The lowest point the water table (at dis- 
tance 3,000 from the center the river) was taken the datum level 
anda slope 1:500 was assumed for the water table. example 
giving prototype dimensions shown Fig. 

The equipotential lines for thesetests were drawn the 
while the streamlines were drawn constructing orthogonal family cur- 
ves the equipotential lines according the rules graphical flow net con- 
struction (3). 


RESULTS 


Results the model tests were summarized inthe form plates and tables 
(25). Two plates, shown Figs. and illustrate cross-sectional flow nets. 
Seepage quantities and distribution were summarized tabular form. For 


948' 


FIG, CROSS SECTION WITH SLOPING WATER TABLE 


each cross section, seepage factor was obtained from measured values 
the upper layer, and from Eq. 


Therefore, the total seepage from one side the channel for given cross 
section obtained multiplying the factor a/b appropriate values 
and The fraction total seepage emerging through the vertical side 
depth was determined wellas the distribution seepage ground surface 
function distance inland from the levee. 

dimensionless scaling factor was computed for each cross section 
that the test data can extended other anisotropies the upper and lower 
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layers. The factor related the permeabilities and scales 
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which ry, and denote the ratio the figure length 
the corresponding prototype distance and other terms are previously de- 
fined. 


ANALYSIS RESULTS 


Seepage flow nets for the cross sections studied depend upon the variables: 
(1) channel depth, (2) channel width, (3) levee base width, (4) anisotropy, (5) 
arrangement and thickness layers, and (6) water table slope. From the test 
results effort was made evaluate the effect each factor seepage. 

Effect Channel Depth.—Discharges 
ratio d/D Fig. general, the discharge increases with channel depth, 
the rate rise being dependent uponthe anisotropy ratio. For increase 
depth ratio 0.25 1.00, the seepage increases for anisotropy ratio 
1:100; increase approximately 45% occurs for increase depth ratio 
from 0.25. For anisotropy ratio 1:10 the increase discharge 
for increase depth ratio from 0.25 1.00 and the same for in- 
crease depth ratio from 0.25. 

Analytic evaluations channel depths and are discussed Ap- 
pendix 

Effect Channel Width.—Most the tests were performed with half- 
channel width equal 300 ft, that direct observation this effect for 
the various anisotropy ratios and thicknesses layers could not made. 
However, was apparent that relatively small portion the total seepage 
occurs through the channel bottom when the upper layer pervious one not 
extending appreciably below the channel bottom. 

tests which the bottom the upper pervious layer extends well be- 
low the channel, the percentage seepage through the bottom increases 
approximately 50%. The percentage depends upon the geometry the cross 
section and the anisotropy ratio. When the upper layer semipervious, most 
the seepage takes place through the bottom thickness this layer 
smaller than the depth the channel. 

The effect channel width for homogeneous anisotropic case was demon- 
strated tests for bottom widths 300 ft, 200 ft, and The resulting 


are plotted against the depth 


seepage factors were 1.08, 1.10, and 1.15, respectively. Similarly, 


small differences were obtained inthe analytictreatment zero 
nel I). 

Effect Levee Base Width.—Because all tests were performed with levee 
base 100 ft, conclusions could drawn the effect this fac- 
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ANISOTROPIC MEDIA 
tor. For the theoretical case channel with zero depth, described Ap- 
pendix the effect levee base width was determined (Fig. 14). For given 
B/D ratio, the discharge decreases with increasing levee base width b/D. This 
trend believed generally true for the investigated cross sections. 

Effectof Anisotropy.—Fig. shows discharge function the anisotropy 
ratio the more permeable layer. general, the seepage factor increases 
the anisotropy ratio becomes smaller. This relation holds both for the ho- 
mogeneous cases (Curve and for two-layer cases (Curves II, and IV). 

was apparent that flow the semipervious layer, when the upper 
layer, mainly vertical that its anisotropy ratio unimportant. Also, 
whenever the semipervious layer lies below the pervious one, the flow trans- 
mitted through very small sothat the effect its anisotropy again neg- 
ligible. 

The distribution seepage through the side the channel depends 
anisotropy ratio. For both the homogeneous and the layered cases with the 
more permeable layer top, the percentage seepage through the banks in- 
creases the anisotropy ratio decreases. When the less permeable layer 
top, the fraction seepage through the banks increases with increase 
the anisotropy ratio. Trends are shown Fig. 10; some cases the effect 
negligible. 

The distribution seepage with landward distance from the levee shown 
Fig. for three tests with homogeneous cross sections and for three tests 
with isotropic semipervious layers top. Although data are limited model 
dimensions required for transformations, these demonstrate that seepage 
distributed further inland the anisotropy ratio decreases. 

Effect Arrangement and Thickness Layers.—In comparing the effect 
the position the relatively permeable layer that the relatively imper- 
meable layer, was noted that the discharge greater when the more per- 
meable layer top. For the same situation, greater amount seepage 
appears near the levee. 

The effect the relative thicknesses the two layers shown Fig. 
for selected tests. These show, for different anisotropies, that the discharge 
decreases thickness upper permeable layer decreases| Fig. 12(a)], 
but the converse true the upper layer the relatively impermeable one 
Fig. 12(b)]. 

Effect Water Table Slope.—The effect was studied two tests. both 
cases the water table slope away from the levee was 1:500, but the total drop 
(in feet) represented different percentages the total head direct com- 
parison between the two tests could not made; however, individual compar- 
isons were possible between these tests and test having similar cross sec- 
tion with horizontal water table. Compared discharge 1.08 for 
the horizontal water table, discharges 1.3 (for water table drop 0.3 
and 3.88 (for water table drop 0.6 were obtained. Slope the 
water table also affects seepage distribution landward from the levee. Seepage 
distributed farther inland with sloping away from the levee than 
for horizontal water tables. The larger the portion the total available head 
above the lowest point the slope, the greater the inland extent seepage. 


CONCLUSIONS 


Model tests demonstrated that seepage from leveed river into adjoining 
low-lying lands affected minor extent variations channel depth 
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and width for the conditions studied. The greater the base width imper- 
meable levee, the seepage. was found that witha smaller anisotropy 
ratio (greater horizontal permeability) the seepage factor was greater 
tended farther inland. Only rough estimates anisotropy are necessary 
define seepage flow pattern. two-layered strata, seepage greater the 
more permeable stratum the upper one and increases with the thickness 
the more permeable stratum. With water table sloping away from the river, 
the seepage quantity and inland movement increases with slope. 

The electric analogy model was found provide versatile means 
for evaluating seepage under variety boundary conditions. The model 
general application, being capable representing infinite number ge- 
ologic conditions. Only lack subsurface data precludes analysis seepage 
from given cross sections. 

Results defined seepage for variety geologic conditions, particularly 
those where two strata govern the flow. means dimensionless scaling 
factor, good estimates can obtained seepage quantity and distribution for 
many other cross sections which only approximate test conditions. 


APPENDIX I.—ANALYTIC EVALUATION EFFECT CHANNEL DEPTH 


Analytic determination seepage practicable for only very few sim- 
plified cross sections. Solutions are presented here for two isotropic cases, 
one with channel depth equal zero and second with fully penetra- 
ing channel D). These supplement the model results and their interpre- 
tation. 

Case Channel with Zero Depth showsthe cross sec- 
tion for this case. With regard tothe resulting flow net, Figs. 13(b) and 13(c) 
are equivalent Fig. 13(a). Fig. 13(c) then used starting point for 
set Schwartz-Christofel transformations. The z-plane shown Fig. 13(d) 
transformed into the -plane Fig. 13(e) means the transformation 


sin and and Fig. 13(e) are defined 
and 


The -plane shown Fig. 13(f) equivalent the shown Fig. 
13(e). this transformation the entire width (2B) the channel consid- 
ered; accordingly, the resulting discharge will correspond the whole chan- 
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nel. The -plane (Fig. 13(g) obtained from the -plane rotating 


The -plane[ Fig. 13(h)] obtained -plane reducing the scale 
afactor Finally, the upper half the -plane transformed into the 
inside the rectangle the -plane Fig. 13(i)] the transformation 


(12a) 


and 


1-s2 


which and and are shown Fig. 13(i). 


The total flow from the channel (2Q) equal the flow through the rec- 
tangle Fig. 13(i) which 


which the hydraulic conductivity the porous media and the head 
difference between the water the channeland the water table. The quantities 
and are shown Fig. 13(i). 

comparision can made between the theoretical result and model tests 
for the same cross section. the above considerations apply isotropic 
media, any given cross section has first transformed into isotropic 


one. The discharge determined the ratio that the scale can re- 
duced arbitrarily. Consider the cross section where in., and 
Then 


D 2 D 
0.471 


and from Eq. 


1.880 


The result the model test was 1.06 For another cross section the 


gives seepage discharges for various ratios and 


The distribution seepage with distance from the levee can also com- 
puted. The complex potential the rectangle Fig. 13(i)] 


and 
(14b) 


Let denote the distance from the origin the horizontal direction Fig. 


(a)]. Inthe transformed z-plane this distance represented along the 
imaginary axis. Along this axis 


e -e 
or 
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The integral must evaluated three separate parts: 


For this case 


channel). 


FIG, 14,—SEEPAGE FROM CHANNEL ZERO DEPTH 
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remains constant and only the value changes. This follows from the fact 
that the integral has real part, which the same that part (1), and 
imaginary part which gives the change potential along the base the 
general, values for and obtained from however, only the values 
are considered here. 


where both integrands (and also the integrals) are real their intervals. 
valuation tables yields 


which 
and 


The following example shows how the preceding analysis can applied for 
determining the distribution seepage. 
Example.—At 100 from the levee the same cross section described 
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(1.775)2 


49.6° 


From Eq. 


F(0.727; 49.6°) 0.92 


proximately 53% obtained the model. 

Case Fully Penetrating Channel sketchof this cross section 
appears Fig. 15. After series transformations, similar those des- 
cribed for Case the following results can written (25) 


and 


The total discharge (for one side) 


The seepage distribution can obtained similar way tothat described for 
Case 


Fig. shows the relationship between the discharge (expressed 


Cases and lead the following conclusions: 


and the width the levee expressed 


Increasing the widthof the channel increases the discharge very rapidly 
first, but has almost effect for 
For fixed ratio the discharge increases asthe ratio decreases. 


fully penetrating channel, its width has effect the discharge. 
The discharge from the channel increases the lengthof the base the levee 
decreases. 


Although these conclusions were obtained from the theoretical cases 
channel with depth and fully penetrating channel, believed that these 
trends remain true also for the cross sections investigated the model. The 
analytic study was based aquifer infinite extent, whereas the model 
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sults thetwo methods reveal that differences were negligible, thus justifying 
the model work. 

general, the theoretical approach described previously applicable only 
limited number simplified cross sections homogeneous aquifer. 
Even for these cases, the determination the complete flow net requires many 
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hours claculation for each case. From practical point view, therefore, 
the electric analogy model was essential investigate the seepage through 
cross sections approximating field conditions. 


The investigation presented this paper covers one portion continuing 
research program hydrology sponsored the Water Resources 
Center the University California. The study was conducted the Hy- 
draulic Laboratory the College Engineering Berkeley, California. In- 
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formation field seepage conditions was supplied the California Department 
Water Resources and the Geological Survey. Josselin Jong 
furnished advice the theoretical computations. 
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FISH PASSAGE THROUGH HYDRAULIC 


Glenn Von Gunten,! ASCE 


SYNOPSIS 


One the major problems connection with water resource development 
the Pacific northwest the United States has been safe fish passage 
dams. Particular concern has been felt regardto fish safety dams 100 
more high. This concept has been greatest factor establish- 
ment the basic plan for development the Lower Columbia and Snake 
Rivers. Also, may that the apparent safety accorded dams under 100 
high has been adverse factor migratory fish protection. There now 
promising evidence that successful passage can achieved through hydraulic 
heads several hundred feet, and that the hydraulic characteristics the 
water passageways are far more significant than heights dams. 

This paper reviews basic concepts fish passage through dams and their 
significance with particular reference given passage Test 
procedures and summarized results fingerling survival passing through 
high speed model turbines and prototype turbine under 450 hydraulic 
head are furnished. The numerous data presented demonstrating correlation 
fish survival with various turbine characteristics show promise leading 
the achievement successful passage through high head turbines. 


INTRODUCTION 


Anadromous fish have had and will probably continue have great 
influence water resources development the Pacific northwest. Hydro- 


open until October 1961. extend the closing date one month, 
written request must filed with the Executive Secretary, ASCE, This paper part 
the copyrighted Journal the Hydraulics Division, Proceedings the American So- 
ciety Civil Engineers, Vol. 87, No. May, 

Presented the August 1960 ASCE Conf, Seattle, Wash. 

Chf., Planning and Reports Branch, U.S, Army Engr, Corps Engrs., Walla 
Walla, Wash. 
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electric power developments, primarily coastal streams and the middle 
Snake River area, have been extensively curtailed due the concern main- 
taining commercial and sport fisheries dependent those streams 
spawning grounds and migration passageways and from the ocean. Curtail- 
ment such development has been exceedingly controversial issue, prima- 
rily because the great economic potential for power development and the 
lack evaluation the effects such developments would have the migratory 
fish population. 

The various species salmon spawn fresh-water streams. The young 
fingerlings move from their hatching grounds the sea where they live until 
maturity. After maturing the ocean, they return the the streams 
spawn and die, thus completing their life cycle. Steelhead trout follow 
similar life cycle but less definable pattern and not necessarily die 
when spawning. 

Dams without adequate fish passage facilities obstruct the natural passage- 
ways the migrants. structure developing hydraulic head only few 
feet may create barrier upstream adult migrants, and reservoir im- 
poundments often submerge their spawning grounds. Downstream migrants 
may exposed severe physical stresses passing through hydraulic 
structures. 


UPSTREAM MIGRANTS 


the past much stress has been placed the question the ability 
upstream migrant adults negotiate fish ladders. However, during recent 
years the question upstream passage has been resolved the extent that 
any one several methods appear feasible. This advance has been 
largely due the Fisheries Engineering Research Program financed and 
sponsored the Corps Engineers, and the activities other agencies. 


appears that the primary problem associated with upstream passage the 
ability intercept attract the migrants entrances providing 
proper locations and hydraulic conditions, rather than the ability the mi- 


grants negotiate ladder once found and entered. 


DOWNSTREAM MIGRANTS 


Young fingerling migrants move the sea with the current. dams they 
are involuntarily swept through spillways, turbines, outlet works. the 
Columbia River, fingerlings vary from in. to5 in. long, except for steelhead 
trout which vary from in. Under natural conditions they migrate 
without being exposed major pressure changes. being swept through 
hydraulic structures dam, they are atthe mercy the prevailing forces. 
The extreme turbulence and pressure changes spillways and turbines in- 
flict severe stresses fish passing through them. Some conditions unques- 
tionably cause injury and death, but for the most part, specific correlation 
these conditions distress has not yet been defined. 

There considerable evidence derived from various research programs 
which indicates substantiates the following statements: 


Static pressures associated with hydraulic heads several hundred 
feet are not fatal. Although under certain conditions distress evident, 
lasting injury apparent. 
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FISH PASSAGE 
The apparent mortality passing over through spillways dams 
relatively low head minor. High free-fall spillways discharging into 
plunge pool also appear provide relatively safe passage. 

Negative pressures cause distress and may may not cause permanent 
injury and death. turbine passage and laboratories negative pressures 
apparently have caused distress, and some cases, death. Most observations 
indicate that, under conditions prevailing dams, lasting ill effects are not 
apparent except areas extreme hydraulic turbulence areas pressure 
near the point vaporization. The length and nature exposure 
both temperature and pressure changes and the degree vacuum, may all 
factors significance. 

Exposure area vaporization and cavitation causes death. 

Mechanical injury may experienced fish passing through turbines 
with small clearances between blades relation size fish. The physical 
injury apparently varies with the relative dimensions and speed. Under proto- 
type conditions, indications are that blade clearances are greater than the 
length fish, the mechanical injury not significant normal turbine speeds. 


DOWNSTREAM MIGRANT PASSAGE THROUGH TURBINES 


The mortality tests which have been performed McNary, Big Cliff, 
Baker, Elwha, and Bonneville dams were primarily designed evaluate the 
extent of, and not explore the cause of, death. Intensive review those 
test results and the research data, failed reveal factual evidence corre- 
lation between mortality and either high static hydraulic pressures negative 
pressures normally experienced throughout the major portion hydraulic 
passageways power plant. Because this, the Walla Walla District, 
Corps Engineers, initiated investigations toward evaluating mortality 
relationship physical conditions. The purpose was find means 
safely passing downstream migrants through turbines. The possibility 
studying mortality passing fish through small-scale turbine performance 
models was investigated, and findings were promising. 


MODEL TURBINE TESTS 


June 1959, exploratory tests were made the Morgan Smith (now 
Allis-Chalmers) Hydraulic Turbine Test Laboratory York, Pennsylvania, 
determine the feasibility investigating mortality model turbines. 

low-head turbine performance test stand (Fig. was used conduct 
the tests. The closed hydraulic system was capable providing turbine 
operating heads ft. Hydraulic control and performance calibration 
the model excellent, this test stand has been used extensively for 
turbine design research and acceptance tests for many years. The McNary 
12-in diameter Kaplan model runner (Fig. 2), intake, scroll case, wicket 
gates and draft tube were used for the tests. 

4-in diameter pipe in. long, with flap valve the lower end and 
screw cap and air line fitting the upper end, served fish lock 
move fish into the model. 4-in pipe, extending downward from the lock, 
led the fish release point approximately midway thevertical opening 
the center intake opening. Nylon bobbinet material was used construct 
net capture the fish after they passed through the turbine (Fig. 3). 
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FISH PASSAGE 


nearly all tests, two species fish were used, generally hundred 
fish each species. similar group each species, utilized control, 
was subjected the same water, transportation, and handling conditions, but 
did not pass through the turbines. Primary species used were largemouth 
bass (Micropterus salmoides) and fathead minnows (Pimephalus notetus), 
both ranging from 1.5 2.4 in. long. 

Fish and various services were obtained from the Pennsylvania Fish 
Commission. The fish were held the Commission’s Pleasant Gap hatchery 
and delivered each testing day the hydraulic laboratory York, Pa. 
the finish each day’s testing, both test and control fish were returned 
the hatchery and placed hatching troughs for holding and recording 
delayed reactions. pathologist the Pennsylvania Fish Commission ex- 
amined test fish the hydraulic laboratory record types injuries sus- 
tained the test fish under different test conditions. Injuries detected 
turbine-passed test fish seemed fall into two categories: strictly mechan- 
ical injuries caused turbine blades, and internal injuries caused rapid 
pressure changes. Mechanical injuries ranged from slight nick some place 
the body, complete decapitation severance body. Pressure-caused 
injuries were characterized most often internal hemorrhages, protruding 
eyeballs, and deflated air bladders. 

general, with the runner setting tailwater level and normal operational 
conditions free from cavitation, the mortality increased indirect relationship 
speed. The mortality experienced was largely due physical injury. 
Artificial braking, operation abnormally low speed for given net head, 
resulted less total mortality than the normal high speed. The mortality, 
however, was greater than that experienced the same speed under normal 
design Similar tests the same net heads and speeds, but with 
the turbine runner set high level reference tailwater, resulted 
high mortality. The high negative draft tube pressures and associated ex- 
treme cavitation conditions caused extensive pressure injuries fatal 
nearly all survivors mechanical injury. summary the test results 
given Table and shown graphically Fig. September 1959 the 
Martin power plant Francis-type runner (Fig. was substituted for the 
McNary runner. significant differences pressure death due 
mechanical injury were experienced normal speed. However, abnormal low 
speeds resulted greater increase mortality than was experienced with 
the Kaplan runner, Apparently, the Kaplan blades simulated conditions 
closer proximity normal speed than could simulated with the Francis- 
type runner. summary the test results using the Francis runner given 
Table 1(b) and shown graphically Fig. 

The tests were highly successful demonstrating that the research 
fingerling mortality passing through turbines practicable. addition, 
the tests appeared quite conclusive respect the following: 


With small model turbine, indiameter, having blade clearances 
less than the length test fish specimens, wide variations survival could 
achieved, depending speed and tailwater elevation. 

Mechanical injury could distinguished from pressure injury, and 
was apparently directly proportional speed, that is, exposure physical 
contact with blades. 

Sufficient survival from mechanical injury was achieved allow anal- 
ysis pressure injury. 
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FISH PASSAGE 


TABLE MODEL TESTS 


Headwater Elevation Revolutions Total 
Elevations Centerline per Survival 
Runner Minute (Test 
Fish) 
(1) (2) (6) 


Experiment No. 


C1(K) OCT 500(n) 
D2(K) JUN 500(n) 
A2(K) OCT 95(r) 
B3(K) OCT 1400(n) 
C1(K) JUN 1400(n) 
B2(K) OCT 500(r) 
B1(K) OCT 1400(n) 
C2(K) OCT 1400(n) 
D1(K) JUN 1400(n) 
D2(K) OCT 1200(n) 
E2(K) OCT 500(r) 
E3(K) OCT 1400(n) 
E4(K) OCT 1300(n) 


Al(F) OCT 
C2(F) OCT 
A2(F) OCT 
C3(F) OCT 
A3(F) OCT 
B3(F) OCT 
B1(F) OCT 
OCT 
D2(F) OCT 
B2(F) OCT 
OCT 
F1(F) OCT 
F2(F) OCT 
E3(F) OCT 
OCT 
D3(F) OCT 
E2(F) OCT 


(McNary Model, in, Adjustable Blades, Size Relationship Model 
(Martin Dam Model, in, Diameter, fixed Size Relationship Model 
reduced speed 


(a) Kaplan 
76 
76 
(b) Francis Runner 
5,41 0.16 5.25 313(n) 
+31.45 1,00 30.45 768(n) 


Extremely high mortality prevailed where adverse hydraulic conditions 
resulted cavitation. 
The marked increase mortality associated with increases negative 
pressures low heads, together with lack any evidence that the positive 
head was influential mortality, encouraged extending the research 
high-head prototype installation. 


HIGH-HEAD TURBINE TEST 


The Fishery agency members the Corps Engineers’ Fisheries 
Engineering Research Technical Advisory Committee expressed interest 
the results the model tests and gave support continued study, which 


NORMAL SPEED HIGH TAILWATER 
REDUCED SPEED HIGH TAILWATER 
NORMAL SPEED LOW TAILWATER 
REDUCED SPEED LOW TAILWATER 


SMALL SYMBOL, INDIVIDUAL TEST 
SYMBOL, TEST AVERAGE 


SURVIVAL 


HUNDREDS 
FIG. 5.—KAPLAN RUNNER 
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led prototype experiments high head power plant. After compre- 
hensive study power installations the Pacific northwest, was decided 
that the city Tacoma’s (Washington) Cushman hydroelectric plant 

was the most suitable project which investigate the effects tailwater, 
partial loadings, and the possibility successful passage fingerlings 
through high-head turbines. Following this investigation project under the 
Corps Engineers, North Pacific Division, Fisheries Engineering Research 
Program established and conducted the Walla Walla District. Fig. 
view the Cushman penstock, looking down toward the power plant 
from point just below the surge tanks. 

The Cushman Project has hydraulic head approximately 460 ft. Effi- 
ciency the units within the range normal practice and cavitation 
characteristics are considered normal possibly subnormal. The turbine 

runner diameter in, developing 37,500 300 rpm with 440-ft head. 
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The runner has fifteen blades. The blade clearance varies from approxi- 
mately in. the inlet side in. 1/2 in. the inner discharge side. 
The draft tube vertical, discharges chamber, and pro- 
vided with discharge The turbine used the tests served 
10.5-ft diameter penstock anet hydraulic head approximately 450 ft. 
The turbine discharge approximately 800 cfs design head and rated output. 

series marks (tattoos) for the fish used the tests were set 
provide certain mark for combinations tailwater eleva- 
tion and wicket gate opening that, regardless when fish was captured, 
the conditions under which passed through the turbine would known. 

Fingerling test fish were inserted penstock El.-402 msl under ap- 
proximately static head through pressure lock (Figs. and 9). The 
test specimens were collected large net the draft tube exit, which 


NORMAL SPEED HIGH TAILWATER 
REDUCED SPEED HIGH TAILWATER 
NORMAL SPEED LOW TAILWATER 
REDUCED SPEED LOW TAILWATER 


SMALL SYMBOL, INDIVIDUAL TEST 
LARGE SYMBOL, TEST AVERAGE 


SURVIVAL 


FIG, RUNNER 


screened the entire discharge exclusive minor leakage around the net 
anchorage frame (Fig. 10). The average recovery all test specimens placed 
the penstocks was 90%, The average recovery all control fish released 
the center the net entrance was 82%. Pertinent power plant data are 
given Table and test results are summarized Table Although test 
procedures preclude refined evaluation, variations results duplicate tests 
appear relatively small. Where variations individual test results ap- 
pear significant, attributed the fact that catches from one test some- 
times overlapped those similarly marked fish from identical test. The 
fish were furnished, marked, handled and personnel the Wash- 
ington Department Fisheries and the U.S. Fish and Wildlife Service. 

Under the more favorable tailwater and hydraulic conditions normal 
turbine operation best efficiency, and average survival more than 75% 
was experienced, adjusted mortalities experienced with control 
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TABLE 2,—PHYSICAL PLANT DATA CUSHMAN NO, 
HYDROELECTRIC PLANT UNIT NO, 


Plant Location 
North Fork the Skokomish River near Potlatch, Wash. 


Manufacturer ------------------------------------ 


Normal Head ------------------------------------ 450 
Maximum Head ---------------------------------- 485 
Discharge Normal Head and Rated Output--------- 800 cfs 
Maximum Tailwater msl 


Minimum Tailwater Elev, -------------------------- 


Turbine Runner Inlet Diameter -------------------- 
Turbine Runner Discharge Diameter --------------- in, 
Maximum Guide Vane Opening --------------------- 1/4 in, 
Minimum Clearance Between Runner 

Vanes Varied from --------------------------- 


Length Total ----------------------------------- 1,310 
Length Lock Scroll --------------------------- 


TABLE 3.—TEST FISH SURVIVALS CUSHMAN NO, HYDROELECTRIC PLANT 
UNIT NO, 


Low Tailwater 
(4) 


Medium Tailwater 
(3) 


High Tailwater 


Gate Opening 
(2) 


(1) 


Average all tests. 


SURVIVAL 


NOTE: 
CURVES REPRESENT AVERAGES 


OF A rs 


TAILWATER RISE FEET 


FIG, 11.—1960 CUSHMAN NO. TESTS, UNIT NO, 


Turbine: g 
Draft Tube with Hydraucone 
in. 
Penstock 
from 10'6"' to 8'0"' 
65% 77.3% 70.9% 65.5% 
80% 73.7% 55.1% 
100% 73.5% 69.1% 


FISH PASSAGE 


Under the most adverse conditions partial turbine loading and low tail- 
water, the survival was reduced approximately 45%. Figs. 11, 12, and 
show relationship survival turbine characteristics and tailwater eleva- 
tions. The only apparent significant variables test conditions were turbine 
output and elevation tailwater. The tailwater designations (Table are 

averages generalized categories associated with high, low, and medium 

tides. Tailwater elevations were further influenced the discharge other 

units; consequently, detailed correlation survival each the variables 


PER CENT 


INDIVIDUAL POINTS ARE AVERAGES 
TESTS FOR EACH GATE POSITION TESTED 
GATE OPENING, PER CENT 


FIG, 12,—1960 CUSHMAN NO, TESTS, UNIT NO, 


MEDIUM 
LOW 


NOTE: 
CURVES REPRESENT AVERAGES 


UNIT EFFICIENCY, PER CENT 


FIG. 13.—1960 CUSHMAN NO, TESTS, UNIT NO, 


cannot made. Variations survival 10% more were experienced 
with each the variables individually. 


SUMMARY 


The exploratory tests survival rates fish passing through model and 
prototype turbines have definitely demonstrated that turbine performance 
characteristics are important, and possibly have greater influence sur- 
vival than all other factors. all test programs, high mortalities have been 
observed passing fish through turbines operating low efficiencies and 
through cavitating conditions. Since the extent the water 
passages can highly influenced design and operational considerations, 
the prospects improving survival are promising. 
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The exploratory research date (1961) has indicated definite avenues 
future study, and also gives much promise that passage through turbines may 
made relatively safe. The investigations the future should aimed 
three problem areas: 


Structural design turbines reduce mechanical injury. 

Physical setting the runner relationship tailwater conditions 
reduce negative pressures. 

Scheduling turbine operation obtain optimum hydraulic conditions 
throughout the water passages. 


Walla Walla District plans conduct additional turbine model tests with 
fish October 1960. Additional tests are also being planned for the spring 
1961 Cushman No, explore further both mechanical and pressure in- 
juries and, possibly, achieve direct correlation between model and proto- 
type results related both pressure and mechanical The primary 
purpose these tests will significance fish size, turbine 
blade clearance, and turbine speed related fish survival. 
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FOREWORD 


The purpose the Task Force Aerated Flow Open Channels 
investigate the status present (1961) knowledge the phenomena aerated 
flow open channels and prepare report the subject. The personnel 
the Task Force and brief record its activities are presented under 
the heading “Acknowledgments.” 

The reason for the interest the phenomenaof aerated flow mechan- 
ical mixture involved with the proper design steep channels chutes. 
air entrained the flow, the mixture increases volume “bulks.” 
The bulking aerated flow requires higher side walls than non-aerated flow. 
Thus, the phenomenon has important economic aspects the design chute 
spillways for earth dams, chutes for wasteways irrigation canals, and high 
overflow spillways. 

The process air entrainment chutes inextricahly involved the 
mechanism the generation turbulence. The actual presence large 
quantities air the form bubbles, doubt, affects the decay turbu- 
lence. The complete problem relatedtothe problem suspended sediment, 
except that solid particles tend settle toward the bottom where turbulence 
generated and air bubbles rise the surface. 

Although much has been learned the mechanism the phenomenon 
experimentation, rational solution has not been accomplished, the knowl- 
edge Task Force members. The objective this Report include 
empirical relationships which should helpful engineers preparing 
economical design channel carry aerated flow. Reliable empirical 


open until October 1961. extend the closing date one month, 
written request must filed with the Executive Secretary, ASCE, This paper part 
the copyrighted Journal the Hydraulics Division, Proceedings the American So- 
ciety Civil Engineers, Vol, 87, No, May, 1961. 
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design criteria are not known have appeared the technical literature 
published the United States. 


HISTORICAL RESUME 


full account the history development the knowledge aerated 
flow not within the scope this Report, nor would serve any useful 
purpose except for historians. brief historical resume, however, 
value afford the engineer aperspective problem, well recog- 
nize outstanding contributions the advance general information the 
subject. short list references appears Appendix followed more 
extended bibliography Appendix 

Three significant events have attentionof the Task Force mem- 
bers. The first important event concerns pioneer researchin Germany, about 
1925. The next two events involve the initiation important research the 
United States and published symposium entrainment air flowing 
water. The occurrence the last two events were stimulated small 
measure the interest the original Committee Hydraulic Research 
which preceded the organization the Hydraulics Division ASCE. 

The Committee Research the Hydraulics Division inherited the 
functions the old committee, upon organization the Division. 1957, 
the name the Committee Research was changed the Committee 
Hydromechanics and new Committee Research was formed. indeed 
appropriate, therefore, that the sponsoring committee this task force should 
maintain continuing interest the difficult problem bulking caused 
aerated flow. 

Instruments measure the amount air mixed with the water form 
important aspect the advance the technical knowledge. Ehrenberg- 
described simple. yet effective, device for measuring the concentra- 
tion air entrained the flow chutes constructed the laboratory 
Versuchsanstalt fur Wasserbau, Vienna, Austria. His apparatus was simply 
tube, pointing into the flow, and communicating through elbow vertical 
glass tube much larger diameter. The Pitot-like tube was rated unaerated 
flow with known velocity measuring the time fill known volume the 
vertical tube. The operation was repeated where the air bubbles 
rose the surface the vertical tube. the increased time fill the tube 
aerated flow, Ehrenberger calculated the air concentrationin the air-water 
mixture. 

The Ehrenberger apparatus produced results which have been questioned 
because the lack similitude into the Pitot-like tube aerated and 
unaerated flow. The apparatus was definitely crude compared with the 
modern electronic method, but ingenious, nonetheless. 

Ehrenberger published empirical relationships defining air concentration 
follows: 


-0.05 


for sin 0.476 


-0.05 -0.74 


Numerals parentheses refer corresponding references Appendix 


} 
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AERATED FLOW 
which the percentage water the toal air-water mixture, denotes 
the hydraulic radius, and the angle the channel slope with the horizontal. 
Ehrenberger considers his formulas valid for 0.30 (0.98 ft). The width 
his wood chute was 0.25 (0.82 ft). Nodoubt substantial portion the air was 
entrained turbulence generated the side walls. Therefore, the formula- 
tion his results would not applicable chute different width. 

Experimental work was initiated the United States about 1939. During 
that year Stevens, ASCE, was chairman the original Committee 
Hydraulic Research. Upon the recommendation ASCE, the Engineering 
Foundation made small grant Lorenz Straub, ASCE, Director, St. 
Anthony Falls Hydraulic Laboratory. ASCE had intended that this would 
“seed” money, stimulate other sponsorship, which indeed proved be. 
Later the funds were supplemented support from the University Minne- 
sota, the United States Navy, and other sponsors. 

The early adjustable metal flume St. Anthony Falls was wide. 
effort achieve true two-dimensional flow, with sidewall effect the air 
entrainment, the flume was rebuilt width in. 1950. 

The most important advance the science aerated flow research 
St. Anthony Falls Hydraulic Laboratory was the development better methods 
measuring the concentration entrained air. The early method measur- 
ing the entrained air volumetrically was give way electronic method 
1950. This apparatus was described paper Lamb and 
Killen (2). Essentially involved pair electrodes set the line the 
flow direction. The electrical resistance the laboratory water without 
entrained air was measured, effect. The concentration air aerated 
flow could then calculated measuring the increase resistance across 
the electrodes. 

Two series experiments were conducted with the wide flume and the 
electrical resistance apparatus. Results series with the smooth steel 
wall were unpublished report Straub and Anderson (3), for 
the United States Navy 1955. The results the test series walls 
roughened with granular particles were reported Straub and Anderson (4) 
1958. Significantly, Straub was awarded Research Prize ASCE 
1958 for his work aerated flow. 

Symposium the Entrainment Air Flowing Water was published 
the Transactions ASCE 1943. The paper Standish Hall (5) re- 
ported field tests three chutes operated the Pacific Gas and Electric 
Company. also included field test data the Kittitas chute the Bureau 
Reclamation, United States Department the Interior (USBR). For reasons 
large width and mimimum amount vertical and horizontal curvature, 
the Kittitas results proved the most desirable for later analysis. 

The tests the Kittitas chute were made Thomas and reported 
the USBR 1938 (6). Thomas measured the water-surface profiles 
well the mean velocity the salt-velocity method. With known water 
discharge, the quantity entrained air could estimated. 


PRESENT STATUS KNOWLEDGE 


Chutes.—For purposes this Report, narrow chute defined 
one which has width less than five times the depth. Air entrainment can 
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ordinarily considered relatively simple two-dimensional problem 
for wide chutes. Three boundaries, the bottom andtwo sides, contribute the 
generation turbulence which causes air entrainment narrow chutes. 
probale that the effect from one side wall reaches nearly the opposite 
wall, thus compounding the air-entrainment effect narrow chutes. Fig. 
demonstrates development the side-wall effect the Uncompaghre Chute 
very clearly. Much the data reported Hall (5) fall the category 
narrow chutes. The discussions this paper embraced attempts generalize 
the results for the calculation air concentration and the resultant bulking 
the flow. 


FIG, 1.—AIR ENTRAINMENT FROM SIDE- 
GENERATED TURBULENCE 


FIG, 2,—AIR ENTRAINMENT CAUSED SIDE-WALL 
EFFECT NARROW CHUTES 
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later years, other attempts generalization the air-entrainment 
phenomena were based empirical conclusions 
were roughly applicable narro they leadto over-design spillway 
chute side walls when applied chutes. Fig. shows the side-wall effect 
the USBR Kittitas Chute. 

Wide Chutes.—Maintaining the definition for narrow chute expressed 
previously, wide chute considered one which has width greater 
than five times the depth. example wide chute may seen Fig. 
which shows the Ft. Peck Close scrutiny the water surface during 
operation led the conclusion that there was negligible amount bulking 
from air concentration. Such observation was impressive when contrasted 
with high air concentration indicated previous empirical criteria based 
narrow chute 


FIG, 3,—LOW AIR ENTRAINMENT FROM BOTTOM-GENERATED TURBULENCE 
WIDE SPILLWAY CHUTE 


Hydraulic design engineers have been concerned that the use design 
criteria based unselected field data for narrow chutes had influenced the 
tendency toward costly over-design side walls. recent years, engineers 
the United States Army Waterways Experiment Station, Vicksburg, Miss., 
(WES), analyzed the smooth-channel data obtained St. Anthony Falls (3) 
combination with selected Kittitas field observations. The Kittitas data were 
selected the basis the observations less than one-fifth the 
width. 
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Straub and Anderson had found the air concentration, two-dimensional 
smooth channel function S/q2/3, which the sine the slope 


and the unit discharge. 
The WES issued hydraulic design criteria chart June, 1957, based 


the foregoing analysis the Straub and Anderson smooth-channel results. 


This chart presented much more economical basis for the design wide 
chute spillway walls than did previous criteria from unselected narrow chute 
results. 

mentioned previously, Straub and Anderson (4), presented the results 
tests rough channels December, 1958. They concluded that for 
artificially roughened channel, air concentration was function 
This function slope and unit discharge was based empirical rela- 
tionship between and which the shear velocity isa so- 
called transition depth defined therein. The new Straub and Anderson data 
rough channels were also analyzed comparison with selected Kittitas 
field results. The large unit discharge the prototype chute Kittitas pro- 
duced small values and consequently small concentrations air. 
The curve best fit for the combined laboratory and field data was deter- 
mined use digital computer andis shown Fig. The standard error 
for air concentration was 0.061. The Task Force believes that this rela- 
tionship offers realistic basis for the design wide chute spillways within 
the limitations the present status knowledge. 

Overflow Spillways.—The relationship represented the curve Fig. 
believed applicable wide ungated overflow spillways, although prototype 
observations are yet available for verification. laboratory study 
Bauer (7) made the State University Iowa, Ames, covered 
investigation the development the turbulent boundary layer over- 
flow spillway. His observations were not carried beyond the point the pro- 
file where the turbulent boundary layer reached the surface. this part 
where the phenomenon air entrainment begins. 

The Task Force not aware any research which deals with the effect 
spillway crest piers air entrainment. Examination the operation 
overflow spillways the field indicate that the turbulence generated the 
pier noses can cause substantial air entrainment. Indeed, believed that 
for modern spillways with large design head compared with the width the 
gate bay, the pier disturbance has greater effect air entrainment and its 
bulking, than does the bottom turbulence alone, The effect piers 
air entrainment the Chief Joseph Dam the Columbia River may seen 
Fig. 


NEED FOR FUTURE RESEARCH 


This Report has presented empirical criterion guide the judg- 
ment the engineers the design wide steep channels. Nevertheless, 
much remains learned the basic mechanics the generation turbu- 
lence, its energy the water surface approached, and the actual 
phenomenon air entrainment. 


The curve best fit was 
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STRAUB AND ANDERSON RANGE, 74 LABORATORY 
OBSERVATIONS ON A ROUGH CHANNEL 


KITTITAS RANGE 
18 SELECTED 
OBSERVATIONS 
ON A CONCRETE 


MEAN CONCENTRATION 
UNIT 


a= ANGLE OF CHANNEL 
BOTTOM WITH HORIZONTAL 


= STANDARD ERROR. 


CURVE BEST FIT FOR ALL DATA: 


MEAN 


FIG, 4,—AIR ENTRAINMENT BASED TESTS ROUGH 
LABORATORY CHANNEL AND FIELD TESTS 
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Theory.—The problem related the more fundamental one the genera- 
tion turbulence the boundary, its diffusion upward through the liquid, and 
the action the vertical components the water surface. More 
needs known the reflection dissipation the vertical components 
the surface with various degrees turbulence and the actual piercing the 
surface turbulence cells. undoubtedly the latter phenomenon which gives 
rise air entrainment. 

The present (1961) knowledge turbulence based primarily the great 
mass observations with the flow air. The work aeronautics not 


FIG, SPILLWAY PIERS AIR ENTRAINMENT 


concerned with the phenomenon turbulence the vicinity liquid with 
free water surface. Unfortunately, the hot wire anemometer, which the 
principal measuring device wind tunnel turbulence studies, not readily 
adaptable use with air-water mixtures. Nevertheless, theoretical studies 
turbulence open channels should start with study the theory and test 
results aerodynamics. 

Laboratory Equipment.—Although the combined effect both side walls 
and bottom are important certain structures, the more simple two- 
dimensional problem wide channel should approached first. dis- 
advantage studying wide channels with steep slopes that large discharges 
are needed obtain substantial depth. The fundamental problem turbulence 


open channels can studied with fairly light slopes. Practically uniform 
flow can obtained with reasonable channel length, the bottom artifi- 
cially roughened. fairly wide channel with rough bottom and smooth side 
walls would then have very little side-wall effect. 

More laboratory work perfect instruments which will measure, 
least, the longitudinal and vertical components turbulence velocities 
water. The Hubbard-Ling hot-film anemometer (8)(9) one approach and 
several laboratories have experimented with pressure transducers mounted 
stream-lined supports. 

The three-dimensional problem could well investigated first studying 
the side-wall effect, alone. wide laboratory flume again indicated 
minimize reflection from the opposite wall. The experimenter would probably 
choose roughen only one wall and maintain the bottom and opposite wall 
smooth. 

The behavior the turbulence from discrete wall protuberances, whether 
knobs vertical strips, the disturbance reaches the surface needs 
studied. Little the effect joint offsets channel walls. The effect 
piers overflow spillways belongs this class problems. The Task 
Force believes that the fundamental research could best performed the 
laboratory, where the geometrical variables may changed with little 
expense. 

Field Tests.—More observations existing prototype structures should 
made and publicized, only profession the important design 
problems involved. The problem holding measuring instruments high 
velocity flows the prototype severe itself. Provisions for measure- 
ment such rudimentary variable velocity are badly needed the pro- 
totype. Velocity observations would yield information for friction-loss studies 
well for basic information onhigh velocity turbulent flow. The provision 
for high-velocity measurement must normally done during construction, 
embedded metal which support the bottom span for attach- 
ing the support from Embedded metal has been installed for the attach- 
ment instrument piers the bottom the Ft. Randall chute spillway. Ob- 
servations this installation must await the availability test flows, which 
common problem with all field tests. 

Air entrainment caused the disturbances crest piers overflow 
spillways may critical condition for the design side walls this type 
spillway. General field observations the pier effect may well serve 
useful guide the design future spillways, until more basic work can 
accomplished the laboratory. 
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The original task force membership was: 
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account the urgency other professional business, DeLapp 
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FORMS BED ROUGHNESS ALLUVIAL CHANNELS* 


SYNOPSIS 


Laboratory experiments large recirculating flume, addition field 
studies, have established that resistance flow and sediment transport 
alluvial channels are related the form bed roughness. The forms bed 
roughness that have been observed the laboratory flume and natural 
alluvial channels can divided into two regimes flow, the basis their 
shape, resistance flow, and sediment transport. These divisions are the 
lower flow regime and the upper flow regime. Between the two regimes there 
transition zone. The bed forms are ripples and dunes the lower flow 
regime and plane bed, standing waves, the upper flow regime. 
The bed form the transition zone variable, ranging from that which 
typical the lower flow regime that for the upper regime. 

Resistance. flow the lower flow regime relatively large and sedi- 
ment transport relatively small. Conversely, the upper flow regime the 
resistance flow relatively small and sediment transport large. the 
transition zone both resistance flow and sediment transport vary appreci- 
ably. 

The variables that determine the form bed roughness are the slope 
the energy grade line, the depth flow, fall velocity effective fall diameter 
the bed material, and the shape the channel. The fall velocity func- 
tion the size bed material, the shape factor the sediment, the mass 
density the sediment, the mass density and the viscosity the water- 
sediment mixture. 
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Laboratory and field studies indicate that the change resistance flow 
that occurs when the form the bed roughness changes from dunes plane 
bed standing waves, accounts for the discontinuity observed stage dis- 
charge relationships those streams that flow both regimes during 
runoff event. 


INTRODUCTION 


Flow alluvial channels complicated the fact that the form the 
roughness elements and, hence, the resistance flow and the sediment trans- 
port are functions the characteristics the fluid, flow, bed material, and 
channel, Because the limited scope most studies and the interdependency 
and the number variables involved there are many conflicting statements 
the literature. This was brought out Norman Brooks,3 A.M. ASCE and 
those who discussed his paper. However, clearer insight into the alluvial 
channel problem results when the various forms bed roughness that can 
occur alluvial channel, the affect the form bed roughness, 
and the variation resistance flow and sediment transport with forms 
bed roughness are known. 

Herein described the bed configurations that have been observed ina 
large recirculating laboratory flume with sand bed and alluvial channels 
the field; the magnitude resistance flow and sediment transport as- 
sociated with these bed forms; the variables that influence the bed 
tion and, where known, how these variables influence the bed form, possible 
relationship that may prove useful predicting the configuration the bed 
introduced and one field situation for which the change bed form 
prime importance described. 

This analysis outgrowth continuing study alluvial channels 
the United States Geological Survey, Dept. Interior (USGS) Colorado 
State University, Fort Collins, Colo. The investigation conducted primarily 
recirculating flume 150 long, deep, and wide. Slope can 
varied from per 0.01 per and discharge from cfs cfs. 
supplement the large flume recirculating flume wide, deep, and 

The depth sand bed the flumes was approximately 0.6 ft. With shal- 
lower depths sand the dunes could not fully develop and portions the flume 
floor were exposed. The 0.45 data quoted this report were previously 


Streams with Movable Beds Fine Sand,” Norman Brooks, 
Transactions, ASCE, 123, 1958, 
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published4 and the 0.28 data are includedherein (Table 1). more detailed 
description the facilities and methods collecting the data are 


REGIMES FLOW 


commonly recognized that four regimes flow occur open channels 
with rigid boundaries.6 Two these regimes frequently occur natural 
alluvial channels. These are the tranquil-turbulent flow and rapid-turbulent 
flow regimes. Tranquil and rapid flow are adequately defined either the 
specific energy diagram the Froude number, When the flow 
tranquil and when the flow rapid. The Froude number defined 
which the velocity flow, feet per second; denotes the gravita- 
tional acceleration, feet per second squared; and refers the depth 
feet. 

Normally the mean velocity and average depth are used compute 
the Froude number. However, important recognize that with the ex- 
treme variability flow conditions that occur the cross section 
alluvial channel may advantageous consider local values velocity 
and depth and the corresponding Froude number help explain observed 
phenomenon. not uncommon confronted with situation where the 
Froude number, based average and less than unity and yet the 
same cross section, the local Froude number may exceed unity. 

The absolute magnitude the Froude number either the tranquil the 
rapid flow regime depends onthe scale the system and only quantitatively 
significant for the system under consideration. For instance, the large re- 
circulating flume which wide, the dunes only occurred when 0.3< F<0.6, 
whereas large deep river, dunes can occur whenF 0.3. Also the 
large flume, the beginning motion occurs 0.15, whereas, 
very small flume using the same bed material and width-depth ratio, the be- 
ginning motion may occur 1.0. 

the water accelerates over the artificial natural humps the 
stream decelerates over the depressions This illustrated 
Figs. 1(a) and 1(b). Fig. 1(a) also illustrates separation zone and the 
existence strong circulation the trough area dune bed and the fact 
that boils appear the water surface just downstream the crests the 
sand waves. The separation and the turbulence generated dissipates consid- 
erable energy which increases the resistance flow. 

When the water decelerates over the humps and accelerates the 
troughs illustrated Figs. 1(c) and Fig. 1(c) illustrates symmetrical 
sand and water waves sinusoidal form that have been described mathemati- 
cally with limited accuracy the flumes these sand and water 


Proceedings, ASCE, Vol. 86, No. May, 1960, 
“Studies Flow Alluvial Channels; Flume Studies Using Medium Sand,” 
Simons, Richardson, and Albertson, Water Supply Paper 1498A, 
Geol, Survey, 1961, press, 
“On the Four Regimes Open-Channel Flow,” Robertson and Hunter 
Rouse, Civil Engineering, March, 1941, 
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waves are commonly observed the rapid flow regime. There little sepa- 
ration and circulation the rapid flow regime when plane bed and standing 
waves exist. these cases, the dissipation energy reflected resistance 
flow primarily the result shear the bed and the formation waves. 


TABLE 1.—BASIC VARIABLES FOR RUNS WITH 0.28 SAND THE FLUME 


0.045 10.73 1.00 Ripple 
13B 0,062 1.00 1,19 Ripple 

0.069 0.86 1,25 Ripple 
0.073 0.59 1,23 Ripple 

0.108 7.19 0.57 1,20 Ripple 
0.116 8.61 1,21 Dune 

0.131 20.39 1.07 1,19 Dune 
9.90 0.65 1,23 Dune 
0.134 1,02 1,21 Dune 
0.136 0.65 1,24 Dune 

0.136 0.88 1,22 Dune 
0.141 11.96 0.61 2.45 1,24 Trans, 
0.150 5.50 0.44 1,56 Dune 
0.153 14,92 0.60 1,31 Plane 
0.172 21.98 0.82 3.35 Trans, 
0.199 21.85 0.72 3.79 Plane 
0.229 15.72 0.55 3.57 1,23 Plane 
0.278 15.70 0.52 3.77 1,22 Plane 
0.470 21.76 0.58 4.69 Antidune 
0.533 0.43 4.50 1,23 Antidune 
0.593 21.34 0.56 4,76 Antidune 
0.815 21.33 0.54 4.93 1,38 Antidune 
8.34 0.30 3.48 1,35 Antidune 

0.930 15.26 0.40 4.77 1,37 Antidune 

1,007 0.57 4.69 1,36 Antidune 

Slope poorly defined. 


cfs sec per sec ture ppm 
(1) (2) (3) (4) (5) (6) (7) (8) (9) 

0.023 7.76 1.01 1,38 10.9 3.3| Ripple 
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BED ROUGHNESS 


With antidunes the resistance flow related the shear the bed, the 
formation waves, and the energy dissipated the breaking waves. 


FORMS BED ROUGHNESS 


The forms bed roughness observed and alluvial streams 
are illustrated Fig. These forms bed roughness were described and 
divided into the tranquil and rapid flow regimes.4 Additional laboratory 
studies the writers and field investigation, indicate that more logi- 
cal classify forms bed roughness into lower and upper regime con- 
nected transition zone than work with the tranquil and rapid flow re- 
gimes that are defined the Froude number. 

The lower and upper flow regime classification upon magnitude 
resistance flow, bed, and water surface configuration, and mode and magni- 
tude sediment transport. the lower regime the bed roughness ripples 
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FIG, 1.—RELATION BETWEEN WATER SURFACE AND BED CONFIGURATION 
THE TRANQUIL FLOW AND RAPID FLOW REGIMES 


dunes both, the resistance flow large and the bed material load 
relatively small. The majority the bed material load moves close the 
bed over the backs the ripples dunes and the crest these roughness 
elements some the bed load avalanches down the faces the ripples 
dunes where temporarily stored causing the ripples and dunes advance 
downstream, The other part the bed material load carried onward, some 
the surface the boil areas. 

When the bed plane when standing waves and antidunes develop the 
resistance flow small, the bed material transport relatively large 
and nearly continuous motion. With ripples and dunes bed load closely 
related the velocity the roughness elements. With plane bed, standing 


gan . 
Tranquil Flow, Alluvial Channel 
(a) 
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waves, antidunes the bed load closely related the velocity the in- 
dividual grains moving contact with the bed. the transition region which 
connects lower and upper flow regimes there discontinuity the re- 
sistance flow and sediment transport relations. 

There are reports the literature aplane bed with bed material move- 
ment which develops before ripples the lower regime flow. Based upon 
experimentation, this seems small flume phenomenon. the large 
flume whenever general movement bed material occurred (d< 0.5 mm) 
ripples dunes ultimately formed. 

Ripples are usually spaced the order 0.5 1.5 and have ampli- 
tudes 0.01 0.1 ft. The spacing the dunes usually greater than ft, 
and their amplitudes range from 0.10 feet depending the depth 


(b) Dunes with ripples superposed 


Boil Boil incipient 
and moving upstream 
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Breaking antidune wave 


FIG. 2.—FORMS BED ROUGHNESS ALLUVIAL CHANNELS 


flow and the characteristics the bed material. For example, small al- 


ft, whereas, the Mississippi River sand waves dune form have been 
recorded Carey, ASCE andM. A.M. ASCE which have 
lengths several hundred feet and amplitudes large ft. Dunes 


defined here are frequently referred bars. 

the transition zone for the same boundary shear the bed configuration 
may range from plane bed fully developed dunes depending primarily 
antecedent conditions. Generally, the bed consisted long dunes small 
amplitude that decreased amplitude the boundary shear was increased. 
However, when the bed form consisted dunes, slope depth could in- 


“Systematic Changes the Beds Alluvial Rivers,” Carey and 
Keller, Proceedings, ASCE, Vol. 83, No. August, 1957. 
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creased relatively large values before the bed form changed the plane 
bed standing waves the upper flow regime. Conversely the bed was 
plane, the slope depth could decreased significantly before dunes de- 
veloped. 

plane bed developed the flume experiments following the transition 
zone with the 0.28 sand. With the 0.45 sand plane bed did not de- 
velop following the transition. Instead, standing waves developed. The standing 


and water waves that occurred with the 0.45 sand and that did not 


occur with the 0.28 sand were sinusoidal form, phase, and did not 
move either upstream downstream with time any great extent. The 
amplitude and spacing these waves depended upon the size the flow 
system and the characteristics the fluid and the bed material, The ampli- 
tude the water waves increased with depth and were 1.5 times the 
amplitude the sand waves. 

The antidunes are very similar standing waves except some the sand 
and water waves continued grow amplitude until they broke. Antidunes 
have been observed natural streams which the median diameter the 
bed materials have ranged size from fine sand coarse gravel. the 
flume studies the breaking the antidune waves usually occured when the 
amplitude the water waves were about twice the amplitude the corre- 
sponding sand waves. this time, the water surface the trough was 
approximately the same elevation the crests adjacent sand waves. These 
sand and water waves moved upstream prior breaking. Either one two 
waves broke several waves broke more less simul- 
taneously. After the antidune waves broke, anewtrain waves developed and 
the antidune cycle repeated itself the waves subsided without developing 
their break point and then reform break abate again. indicated 
ASCE the Froude number, where antidunes first form, 
increases with increasing size bed material. 

the large flume the steeper slopes, the antidune activity was the 
form chutes and pools. Flow was rapid and accelerated the chutes and 
flow was tranquil the pools, Fig. 2(h). 

When antidunes break they resemble the hydraulic jump and can analyzed 
with some success such, particularly with two-dimensional flow. Large 
quantities sediment are carried into suspension the breaking antidune 
wave. 

the breaking wave and immediately downstream, the velocity rela- 
tively small. The storage water the short reaches where the waves are 
breaking causes the discharge vary with time and this action tends set 
slug flow. Mendano Creek, which located the San Luis Valley 
Colorado, antidunes set slugs water that had depth about twice average 
depth and traveled down the channel spaced intervals about 350 
nearly fps. the time this field observation, the average discharge 
was approximately 120 cfs, the median diameter the bed material was 0.3 
mm, and the channel slope was 1.67%. 


VARIATION RESISTANCE FLOW WITH BED ROUGHNESS 


The variation resistance flow, measured using the Darcy-Weisbach 
with the form bed roughness and the influence size bed material 


“Hydraulic Criteria for Sand Waves,” Langbein, Transactions, Amer. 
Geophysical Union, 1942, 615-618, 
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and slope the form the bed roughness, resistance flow, and sediment 
transport are illustrated Fig. and Table 

The resistance flow relatively small for plane bed prior beginning 
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FIG. 3.—VARIATION AND CONCENTRATION WITH CONSTANT DEPTH 


plane bed standing wave case where there motion bed material. This 
understandable because, the first case, resistance results from grain 
roughness that fixed; whereas, the second case the grain roughness 
not fixed but moving and form rolling friction exists. 
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general motion the bed material was initiated increasing the 
boundary shear, ripples formed that caused large increase resistance 
flow. With sufficient increase slope, and depth constant, dunes formed, 
and with further increase slope transition conditions developed. The ripples 
formed from the 0.28 bed material had slightly larger resistance 
flow than the ripples formed 0.45 bed material. Conversely, the dunes 
formed from 0.45 bed material had greater resistance flow than the 
dunes formed 0.28 bed material (Table 2). 

the transition zone, because the bed form was unstable, the resistance 
flow was extremely variable. general, the magnitude the resistance 
flow was between that for fully developed dunes and plane bed. The vari- 
ation resistance flow the transition zone was larger for the 0.45 
bed material than for the 0.28 bed material. Also note that the transition 
zone occurred much flatter slope when using the 0.28 sand. 

With still further increase boundary shear the upper flow regime was 
reached. this regime flow the bed roughness was plane bed, standing 
waves, antidunes some combination these. Resistance flow was 
minimum for given bed material throughout the plane bed standing wave 
range antidunes develop, the resistance flow increased slightly 
and continued increase with further increase boundary shear. The rate 
increase resistance flow was much more rapid for the finer sand and 
was reflected directly increased antidune activity. 


VARIATION BED MATERIAL LOAD FUNCTION BED FORM 


The large flume studies indicated that the bed material transport varied 
with form bed roughness and size sand, shown Fig. and Table 
The transport for the ripples reached maximum 150 ppm for the 0.28 
bed material and 100 ppm for the 0.45 bed material. With dunes 
the form bed roughness, the maximum concentration bed material trans- 
port was about 800 ppm and 1200 ppm, respectively, for the 0.28 sand 
and the 0.45 sand. 

For the 0.28 sand, the bed material load the upper flow regime 
ranged from 1500 ppm for the plane bed 42,000 ppm with antidunes. The 
range concentration bed material discharge for the 0.45 sand was 
from 4,000 ppm with standing waves 15,000 ppm for the antidunes. 

The bed material load the transition zone was unpredictable but varied, 
for both sands, from slightly greater than the maximum concentration for 
the dunes concentration that was characteristic that for the plane bed 
the standing waves. 


FACTORS INFLUENCING FORM BED ROUGHNESS 


The variables that fix the form bed roughness alluvial channels are 

which 


and refers the depth, denotes the slope the energy grade line, 
defines the median diameter the bed material, measure the size 
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distribution, refers the mass density the water-sediment mixture 
(normally taken that for water unless considerable fine sediment pres- 
ent), the mass density the bed material, denotes the acceleration 
gravity, the fall velocity the bed material, describes the seepage 
force, shape factor the cross section, shape factor the 
particle, defines the viscosity the water-sediment mixture, and refers 
the effective median fall diameter the bed material. addition tem- 
perature, also depends the concentration fine sediment. 

The full significance all these variables has not been revealed (as 
1961) because the investigations have not been broad enough include all 
conditions. However, the effect some the variables are least quali- 
tatively known and the following paragraphs they are summarized. Re- 
ferring Eq. velocity bed material load can substituted for bed 
roughness. 

Slope S.—The effect slope was illustrated Fig. With depth constant 
change slope increases the boundary shear. With sufficient increase 
shear the form the bed roughness changes. develop given form bed 
roughness alluvial channels that has the samebed material, the slope must 
large the discharge small the discharge large. 

Depth D.—The effect depth not well defined because the range 
depth which has been studied small. However, the experiments conducted 
the large flume the variation depth has been adequate indicate that 
change depth can change the bed form. natural alluvial channels which 
slope essentially constant bed form changes have been reported that have 
resulted from changes depth.10,11,12 However, there combi- 
nations slopes and sizes material where the bed form will not change 
with large change depth. 

addition the changes that can occur with changes depth, 
there are the uncertain effects changes depth bed material load and 
resistance flow when the bed form does not change. With ripples there 
relative roughness effect. The resistance flow decreases with increase 
depth. Whereas with dunes, the flume data indicate that resistance flow 
increases with increase depth. However, has not been verified that re- 
sistance flow with dune bed form increases continuously over large 
range depths. 

Bed material load also affected changes depth flow. Based upon 
recent study available field data, has shown that because 
velocity increases with depth and concentration decreases with depth, 
peculiar but logical relationship exists between sediment load, depth flow, 
and velocity. Within certain range velocities, velocity held constant, 
the bed material load increases with increasing depth, within another range 
velocities, bed material load decreases with increasing depth. 


“Discontinuous Rating Curves for Pigeon Roost and Cuffawa Creeks Northern 
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86, June, 1960, 
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The effect the combination depth and slope the shear term, DS, 
ASCE, and Ning Chien, mention few. All have indicated 
that shear important parameter, whereas stated “that for the 
laboratory flume was found that neither the velocity nor the sediment 
discharge concentration could expressed function the 
bed shear stress, any combination depth and apparent 
Fig. that what Brooks states true for the transition zone. However, 
the upper and lower flow regimes shear closely related velocity and 
sediment concentration. 

Size Bed Material.—The effects the physical size bed material 
the bed form are (1) its effect grain roughness, and (2) its influence 
the fall velocity the effective fall diameter, which measure the mo- 
bility the particle. The mobility the particle, when related the fluid 
forces, the utmost importance determining whether the bed form 
ripples, dunes, plane bed, standing wave, antidunes. 

the second effect the physical size the bed material, measured 
fall sieve diameter, important only related the 
other variables Eq. However, fall diameter has the advantage elimi- 
nating the shape factor and density the particle variables. Knowing the 
fall diameter, the fall velocity effective fall diameter the particle, 
any fluid any temperature, can computed the fluid’s density and vis- 
cosity are 

The importance size bed material the form bed roughness was 
indicated the experiments using the 0.45 and the 0.28 sands. The 
length the dunes (crest crest trough trough) observed while ex- 
perimenting with the 0.28 sand were much longer than the dune lengths 
for the 0.45 sand. Although the amplitude the dunes for the two sands 
were the same, their shapes were different. The dunes formed 0.45 
sand were more angular. The large spacing betweenthe dunes and the smaller 
separation zone and turbulence the troughs accounts for the fact that some 
ripple runs exhibited larger value than the dune runs where.the 0.28 
sand was used. The flow over the relatively long plane backs these dunes 
was quite efficient, reducing the over-all resistance flow value less 
than that computed for the ripple case. The dunes, using the 0.45 sand, 
had much larger values than the ripples the same sand the dunes 
the 0.28 sand because the closer spacing the dunes, the large sepa- 
ration zone the dune troughs, and the greater magnitude turbulence. The 
increase resistance flow that occurred with increase sand sizes 
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ment,” Shields, 1936, Trans, from German Poh and van Vehelin, 
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“The Flow Cohesionless Grains Fluids,” Bagnold, Trans., 
Royal Soc, London, Vol. 249, No. 964, 1956, pp. 235-297, 

Discussion Einstein and Chien “Mechanics Streams with Movable 
Beds Fine Sand,” Brooks, Transactions, ASCE, Vol. 123, 1958, pp. 553-562, 
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BED ROUGHNESS 


different from the trend reported Leopold, ASCE, and 
Maddock, 

The change from lower upper flow regime occurred much smaller 
shear value, and the rate change was faster for the finer Also, using 
the fine sand the form bed roughness changed from dunes plane bed 
when 0.4, and the bed remained plane the boundary shear was in- 
creased until then antidunes developed. Using the coarser sand, the 
dunes decreased amplitude with increasing boundary shear when 0.6 
and transition conditions then standing waves developed. 

That plane bed developed with tranquil flow with the 0.28 sand and 
not with the 0.45 sand attributed the greater mobility the finer 
sand and the limitations the experimental equipment. The boundary shear 
required eliminate the dunes was sufficiently small for the more mobile 
finer sand that the change could occur depth and slope such that the 
Whereas, with the coarser sand the discharge capacity the flume 
system was not sufficient eliminate the dunes depth and slope where 
the change would occur with Field observations indicate that streams 
with 0.45 coarser bed material, but having large depths and flat slopes 
may have plane bed form bed roughness when 

Fall Velocity and Effective Fall Diameter.—As indicated the preceding 
section, the fall velocity the effective fall diameter the bed material 
major importance determining the bed form, resistance flow, and 
sediment transport. Because the size concept entrenched the litera- 
ture, effective fall diameter used preference fall velocity. Effective 
fall diameter sediment particle defined the diameter sphere 
having specific gravity 2.65 and fall velocity quiescent distilled 
water infinite extent temperature 24°C equal the fall velocity 
the particle falling alone any quiescent stream fluid stream temperature. 
This definition neglects the unknown effect turbulence fall velocity. 

The effect the various factors given Eq. fall velocity are well 
known although they are often overlooked. For instance, known that vis- 
cosity and mass density fluid varies with temperature and affects the fall 
velocity, but the fact that fine sediment affects the viscosity18,20,21 and the 
mass density the fluid are often overlooked. increase concentration 
suspended fine sediment (wash load) adecrease temperature increases 
the viscosity and specific density the liquid and decreases the fall velocity 
and effective fall diameter the bed material. Particles different physical 
diameter have the same effective fall diameter their fall velocities are the 
same. With the same effective fall diameter the forms bed roughness, re- 
sistance flow, and bed material load should approximately the same 
spite the difference physical size. 

With respect temperature effect there have been conflicting statements 
the literature. Vito Vanoni, ASCE, and Brooks22 reported that with 


“The Hydraulic Geometry Stream Channels and Some Physiographic Implica- 
tions,” Leopold and Maddock, Jr., Paper 252, Geol. Survey. 

“Fluidity and Plasticity,” Bingham, Graw-Hill Book Co, Inc., New York 

Properties Disperse Systems,” Hermans, ed., North-Holland Pub- 
lishing Amsterdam, The Netherlands, 
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4 


100 May, 1961 


Hubbell23 reported that with increase temperature there 
increase the resistance flow. These two apparently contradicting state- 
ments, that show the effect temperature variation, are easily explained 
considering the forms bed roughness which existed the alluvial channels 
question, and the effect temperature change the fall velocity and, 
hence, effective fall diameter the bed material particles. 

decrease temperature increases the viscosity the water and de- 
creases the fall velocity the sand. Thus, decrease temperature de- 
creases the effective fall diameter given sand. Consequently, sand 
bed covered with ripples and the temperature the water decreased, 
the mobility the particles increased due the decrease effective fall 
diameter the sand, larger ripples form, and resistance flow increases. 
the other hand, the form bed roughness near the transition 
zone and there reduction temperature, the reduction effective fall 
diameter causes the dunes shear out greater extent perhaps causes 
the bed become plane, which (ineither event) accompanied decrease 
resistance flow. Both these phenomena are reversible. These two 
variables (temperature and fine sediment) have been studied the USGS. 

Seepage Force.—In alluvial channel there usually either inflow 
outflow from the channel through the bank and bed material that causes seep- 
age forces because there loss head. Although the effect seepage 
forces bed form has not been studied known that there inflow, 
the seepage force acts reduce the effective weight the sand. there 
outflow from the channel, the seepage force acts the direction flow and 
increases the effective weight the sand. With inflow the extreme case, 
the seepage forces can large enough set quick-sand condition and 
the weight the material equilibrium with the seepage forces. The 
seepage force reduces increases the effective size the bed material 
contact with the bed changing the effective weight the sand the bed. 
Once particle sand loses contact with the bed, the seepage forces longer 
have direct effect it. 

Gradation the Bed Material.—The comments concerning the effect the 
size the material, fall velocity, and seepage force indicate that the bed 
material can adequately described fictitious diameter. This obviously 
not the case. Although, extensive studies into the effect gradation bed 
material form bed roughness, and resistance flow have not been made, 
would appear that alluvial stream with bed material composed only 
one size would not react fluid forces the same the bed material had 
the same median diameter but large range sizes. 

Vanoni and Brooks22 investigated the effect gradation bed material 
bed material load and concluded that for given median diameter bed 
material the load increases the range sizes gradation increases. 
Einstein24 also recognized the importance the gradation bed material. 

Shape Factor the Cross Section.—The effect channel shape rather 
nebulous. However, wide shallow channels have been simulated and under 
certain conditions two more forms roughness have developed simul- 
taneously, for example, small dunes part the cross section and plane 


“Investigations some Sedimentation Characteristics Sand Bed Stream,” 
Hubbel al., Open-File Report, Geol. Survey, Lincoln, Nebr. 

“The Bed Load Function for Sediment Transportation Open Channel Flows,” 
Einstein, Bulletin 1026, Dept. Agric., 1950, pp. 
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bed, standing waves antidunes the other. The resistance flow vari- 
able and difficult evaluate whenthere multiple roughness. The resistance 
flow dependent the percent the bed covered with each bed form. 
Multiple roughness and the variable resistance flow associated with this 
condition has been observed alluvial streams the writers. 

Current (1961) studies resistance flow and sediment transportation 
indicates that the forms bed roughness for given alluvial bed material may 
quite different flumes different size, other conditions being essen- 
tially the same. 

Using the large flume, the forms bed roughness observed the two 
regimes flow were Fig. These forms bed roughness 
seem agree with field conditions. Using the 0.45 bed material the 
small flume the observed bed configurations were different. When bed move- 
ment began, ripples did not form the small flume. fact, ripples never 
developed this flume with this size bed material. When bed shear the 
small flume was increased sufficiently, the plane bed was replaced dunes. 
These dunes were more two dimensional than those which occur the large 
flume and these dunes never, any time, had ripple roughness superposed 
them was the case the large flume. 

the rapid-flow regime, conditions were quite similar both flumes. 
However, the resistance flow the smallflume was much greater because 

flow was more nearly two dimensional which caused the total width the 

flume occupied the waves and antidunes that developed. 

The results this comparison indicate the many problems that the ex- 
perimentor faces who tries effectively utilize the data collected various 
investigators from their flumes different size and design. 


PREDICTION FORM BED ROUGHNESS 


Thus far, completely adequate method predicting form bed rough- 

ness has been developed. Albertson, ASCE, and others presented25 
method suitable for predicting bed forms for the laboratory case, but 
was not completely satisfactory for the field case. Similarly pre- 
sented figure relating Shield’s parameter 7d, the Froude number, 
and the form bed roughness. This figure was not satisfactory for the 
laboratory case but was superior for the field case. another possibility 
Fig. was developed which modification the figure Hsin-Kuan Liu, 
ASCE for ripple formation.27 Fig. separates the bed forms for the vari- 
ous sand sizes for both streams and flumes, fairly well. Because the fall 
velocity used, also takes into account all viscous and gravitational forces 

that appear affect the bed form. However, Fig. may not separate the bed 
forms when larger depths are studied. Until additional data are available that 

cover adequate range depth, will difficult develop the criterion 

necessary predict form bed roughness alluvial channels with confi- 

dence. 


Sediment-Ripple Formation,” Hsin-Kuan Liu, Proceedings, ASCE, Vol. 84, No. 
February, 
“Total Sediment Transport Alluvial R.J. Gardi, thesis presented 

Colorado State Univ., Fort Collins, Colo., 1960, partial fulfillment the re- 

quirements for the degree Doctor Philosophy. 

“Mechanics Sediment Formation,” Hsin-Kuan Siu, Proceedings, ASCE, Vol. 
86, No. May, 1957. 
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the meantime, the possibility predicting the bed configuration based 
upon the mechanics surface waves and the appearance the water surface 
should considered. This method not applicable design but ex- 
tremely useful for analysis. However, some training required the part 
the individual before can apply effectively. 

Whether the flow rapid tranquil can determined observing the 
direction artificially induced surface wave using the wave 
celarity concept. the disturbance generated oppose flow moves upstream 
with respect fixed point observation, flow tranquil, and 
swept downstream, flow rapid. 

also usually possible determine form bed roughness and whether 
the lower upper flow regime exists alluvial channel observing the 


FIG, 4.—FORMS BED ROUGHNESS ALLUVIAL CHANNELS 


condition the water surface. the upper flow regime, the water surface 
likely appear relatively smooth and fast there will some evidence 

standing waves and antidunes [as illustrated Figs. 2(e), 2(f), and 
which immediately fix the general form bed roughness that exists. 

the lower flow regime the velocity will relatively slow and the major 
rises water surface (usually boils) will out phase with the dunes. For 

rippled bed, the water surface will quite plane and placid except for very 
small depths. this case small waves willbe generated the water surface 
the sand ripples. With dune bed configuration, there will turbulence 
generated the water surface the form boils downstream the dunes. 
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Usually the color the water different the boil area due the large 
concentration suspended sediment carried the water surface within the 
boil. The strength the boil activity dependent the magnitude, spacing 
and shape the dunes, and the depth flow. the transition zone the water 
surface will rather plane with minor boils appearing the water surface. 

The classification stream into the proper regime flow and form 
bed roughness observing the water surface combined with some knowledge 
the bed material would give some information about the transport sedi- 
ment and resistance flow. This information could then applied the 
interpretation stage-discharge relations and other flow phenomena. 


FIELD PROBLEMS 


The conditions observed the field and the large flume agree quite 
well, However, the slope energy grade line through given reach 


Transition 
Ripples Bed, Standing 
Dunes Waves, Antidunes 


Ripples and Dunes, 


Falling Stage 


Stage, ft. 


- 
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Rising Stage 


(a) (b) 


FIG, 5.—TYPICAL QUALITATIVE STAGE—DISCHARGE CURVES 
FOR ALLUVIAL CHANNELS 


channel varies within narrow limits and the depth varies greatly the 
discharge changes; but the flume studies the slope can varied widely 
and depth variations are small due the physical limitations the pumping 
plant and the There ample evidence that change depth can shift 
the form bed roughness from dunes transition, plane bed, and antidunes. 
The reduction resistance flow which occurs dunes are eliminated 
accounts for the discontinuity stage-discharge relations which have been 
observed the laboratory and for many natural streams.10,11,12 Typical 
qualitative stage-discharge curves which illustrate the effect change bed 
roughness stage are presented Fig. Einstein and their 
discussion Brook’s recognized that the depth-discharge curves 
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illustrated Fig. could occur. However, present evidence, they 
thought these curves would only rarely ever occur nature. 

Fig. 5(a) the break which shown caused change bed rough- 
ness from dunes plane bed standing waves. The break occurs the 
rising stage larger discharge than the falling stage. There usually 
considerable scatter the stage discharge relation with ripples and dunes 
the bed due the wide variation resistance flow. Under certain cir- 
cumstances, this variation can explained. For example, Fig. 5(b) the 
magnitude resistance flow lags the actual discharge, that is, the change 
bed roughness lags the change discharge. This results smaller re- 
sistance flow and smaller depth than would normally occur for equilibrium 
flow. The reverse occurs the falling stage. this instance, the peak 
discharge, assume the form bed roughness large dunes. discharge 
decreases, the large dunes are not reduced inamplitude and spacing rate 
keeping with the reduction discharge, and the resistance and depth are 
larger than for equilibrium flow the falling leg the hydrograph. This 
results loop type stage-discharge curve which resembles hysteresis 
curve. 

The loop curve Fig. 5(b) may reversed. Then the falling stage points 
would plot below the rising stage points the two curves may coincide with 
the falling stage points plotting the same curve the rising stage points, 
the rising and falling stage curves may cross. The form the curve which 
results depends upon the rate change discharge with time, the magnitude 
the discharge, the characteristics the bed material, the effect the fine 
sediment load, and the characteristics the channel. Under certain circum- 
stances, temperature could play the same role the fine sediment. 


SUMMARY 


The forms bed roughness that exist alluvial channel because 
similarities form, resistance flow, and sediment transport are divided 
into categories lower flow regime, transition zone, and upper flow regime. 
Based upon laboratory and limited field investigations, the major forms 
bed roughness the normal order occurrence with increasing slope with 
depth, bed material, temperature and concentration fine material constant 
are follows: 


Lower flow regime 


Ripples, 
Dunes with ripples superposed, and 
Dunes. 


Transition (Bed roughness ranges from dunes plane bed standing 
waves) 

Upper flow regime 

Plane bed, 

Standing waves, and 


Resistance flow and sediment transportation vary greatly form bed 
roughness changes. Consider the two sands with median diameters 0.45 
and 0.28 which have been thoroughly investigated large recircu- 
lating flume. For dune bed configuration fvaries from 0.0489 0.1490, and 
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the corresponding bed material load ranges from 100 ppm 1,200 ppm. For 
the antidune condition, varies from 0.0247 0.0672 and bed material load 
ranges from 5,000 ppm 42,000 

The form bed roughness depends primarily the slope the energy 
grade line, depth, fall velocity effective fall diameter the bed material 
and shape the channel. the transition zone the bed form also depends 
the antecedent form the bed. There are several other variables whose 
effect the bed form may only minor importance. They are the seepage 
force, the gradation the bed material, and the cohesive properties any 
fine sediment the bed, and the physical size the sand grain rough- 
ness. 

has been determined that fall velocity effective fall diameter has 
significant effect the bed form. Sands with large difference physical 
size have been found have the same bed form, resistance flow, and 
sediment transport when their effective fall diameters were the same. The 
factors affecting fall velocity have important effecton the bed form. These 
are the size the bed material, viscosity, and mass density the water- 
sediment mixture, mass density the bed material, and the shape the sand 

Knowledge the forms bed roughness that develop alluvial channels 
and the variables which influence bed roughness essential the theory 
resistance flow alluvial channels, the design stable channels and the 
refinement sediment transport theories. 

further improve existing concepts flow alluvial channels 
important proceed with the following studies: 


The effect depth flow; 

the effect size bed material using larger sand sizes than 0.5 mm; 

the effect gradation bed material; 

the mechanics related the meandering and braiding streams; 

the effect width-depth ratio; 

the mechanics flow around bends alluvial channels establish 
stability criterion and magnitude the energy losses; and 

the mechanics diffusion, dispersion, and turbulence alluvial 

channels, 


The development new instruments essential part future work. 
Developing methods measuring turbulence, instantaneous velocity profiles 
and boundary shear open channels necessary for the continued develop- 
ment the theory flow open channels. 
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ESTIMATING POTENTIAL 


SYNOPSIS 


Methods computing potential evapotranspiration analytical procedures 
have been based the application the turbulent-transport and energy- 
balance concepts. Empirical formulas, correlating some temperature function 
and adjusting for day-time hours, have proved valuable practical utilization. 
This latter approach has been used formulate simple computational pro- 
cedure whereby average daily potential evapotranspiration represented 
proportional the product day-time hours squared and the saturated water 
vapor concentration (absolute humidity) the mean temperature. The day- 
time factor was determined from consideration the disparity between net 
radiation and temperature, latitudinally, and the fact that transpiration re- 
stricted during darkness since the leaf stomata are closed. 

Computed values potential evapotranspiration obtained the new pro- 
cedure are compared with those obtained the more complex Thornthwaite 
method (24)2 and other methods currently employed. General applicability 
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seems justified from comparisons between observed and computed values 
potential evapotranspiration, both yearly and seasonal basis. 


INTRODUCTION 


The determination evapotranspiration professional civil 
engineers, hydrologists, meteorologists, and agriculturists. these and 
others, the problems water and the conservationof water supplies—in view 
the “population explosion” and the increasing water resources— 
more than academic. 

Since evapotranspiration represents the transport water from the earth’s 
surface back the atmosphere—the reverse precipitation—it assumes 
important role the hydrologic cycle and given careful consideration 
assessing the water resources particular region. Given estimates 
evapotranspiration then, through appropriate accounting system soil 
moisture, becomes possible indicate the water balance area, 

Evapotranspiration controlled meteorological factors-radiation, wind, 
temperature, and humidity; soil moisture tension; soil moisture conditions; 
and plant physiological factors. Its direct measurement difficult and climatic 
data have been used number relations (2,7,20,24) obtain useful esti- 
mates. These procedures rely either the correlation between water evapo- 
rated from open water surface vegetative area empir- 
ical relationship with climatic parameters and day-time hours. When soil 
moisture not limiting, evapotranspiration almost entirely dependent upon 
the amount solar energy the surface and the resulting tempera- 
tures (27). Such condition optimum soil moisture results “potential” 
evapotranspiration. 

The concept potential evapotranspiration was introduced 
Thornthwaite (24) and understood represent the rate which water 
would removed from extended area covered actively growing 
green crop, completely shading the soil, and with non-limiting supply 
water. The rate evapotranspiration not influenced the advection 
air that has been heated over drier areas, and the ratio energy utilized 
evaporation that heating the air remains constant. Kohler, ASCE 
(13) has suggested that the evaporation from large water surface must 
equivalent potential evapotranspiration. Inthis respect, the several methods 
adaptive computing evaporation from open water surface seem applicable 
the estimation potential evapotranspiration, provided the necessary data 
may obtained. the other hand, soil and plant factors may significantly 
alter the water lost transpiration that simpler methods based tem- 
perature and day-time hours may effective. The purpose this paper 
formulate such functional relationship. 


ANALYTICAL METHODS 


The turbulent-transfer approach the determination reservoir evapora- 
tion has received considerable attention. This method views evaporation 
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equal the flow water vapor through the surface layer air and deter- 
mined the product the vertical gradient vapor pressure and the rate 
mixing the air. The rate mixing, expressed coefficient diffusion 
austauch coefficient, depends upon the rate change with height wind 
speed. Such equations derived (22), (19), (25) could used for the deter- 
mination evaporation over any kind surface where the difference sur- 
face properties expressed through roughness constant (22). The practical 
application this approach severely restricted due the accuracy re- 
quirements measurement wind and humidity two levels determine 
the rate change wind and humidity with height (26). 

the extensive Lake Hefner study (29) lake evaporation, simple 
empirical formula 


relating evaporation, the difference vapor pressure-saturated the 
surface temperature, and the atmospheric vapor pressure, (as Dalton 
formulated 1802); adjusted f(u), and proportionately con- 
stant, gave results comparable those obtained the complex mass- 
formulas. Surface temperature required for this simplified ap- 
proach that little use since pan evaporation usually measured 
whenever the water temperature observed. The Lake Hefner study also 
confirmed that the energy-budget approach, when appropriate data are avail- 
able, may used successfully estimating lake evaporation. 

The problem the energy balance approach estimating evaporation 
measure assess all other sources and sinks energy and leave evap- 
oration the only unknown. The complete energy balance crop volume 
has been expressed Tanner (23) but may ordinarily considered 
consisting net radiation, R,; heat flux the soil conduction, heat 
flux the air convection, and the heat vaporization water, ex- 
pressed 


Usually the soil heat flux neglected but may measured directly such 
the net radiation while the sensible and latent heat terms can separated 
using the Bowen ratio method (23). 

There has been considerable discussion the validity equating the 
constants the heat and vapor transport equations the Bowen ratio and 
apparently such assumption only valid for conditions over water sur- 
face. using the energy budget under conditions restricted soil moisture, 
forced estimating evapotranspiration, Ray Linsley, ASCE (15) 
states that the Bowen ratio longer useful means eliminating the 
sensible heat exchange from consideration. 

determining the evapotranspiration from irrigated plot, Maurice 
Halstead (8) has illustrated that combination both the turbulent-transfer 
and energy approach required. evaluate the “oasis-effect,” Halstead 
recast turbulent-transfer equations that they employed temperature and 
humidity differences between the surface and one height, and introduced 
another aerodynamic equation indicate the height which the influence 
the oasis extends (8). Assuming conditions net radiation, vertical wind 
gradient, roughness, temperature, and humidity such that the evaporation 
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from 6-ft tank, when surrounded completely dry area, was 0.45 
per hr; then, increasing the size the “oasis” 1-mile field, the com- 
puted evaporation rate was reduced less than1/3 that the 6-ft tank. This 
extreme example demonstrates that actual evaporation can exceed that re- 
sulting from net radiation by, considerable amount, depending upon the size 
the field and the difference between moisture conditions the field and its 
surroundings, Even for field cotton extending miles upwind, observa- 
tions confirm that energy was gained from some drier area upwind (14). 

Penman (20) has successfully combined the energy budget and aero- 
dynamic approaches compute open water evaporation. For the energy bal- 
ance estimate, Penman neglected net heat flux the water and considered 
the Bowen ratio remain constant: and assuming identity the heat 
transfer and vapor transfer coefficients the aerodynamic relations, the 
need for knowing the surface temperature was eliminated. Conversion fac- 
tors have been applied these hypothetical open water evaporation values 
estimate what termed potential evapotranspiration, obtained com- 
paring tanks with turf and water; and for 2-yr period observation they 
were found range from 0.6 during the winter 0.8 the summer, 
Rothamsted, England. The dominant factor determining the seasonal cycle 
the conversion factor the range day-time hours which leads bigger 
summer values and smaller winter values more northerly latitudes but 
only small annual range low latitudes (21). Such coefficients only convert 
pan evaporation equivalent area well watered vegetation and, therefore, 
cannot represent potential evapotranspiration due the “oasis” effect that 
increases with the aridity locality. Besides, the Bowen ratio takes 
larger values arid climates (23) and the empirical aerodynamic relation 
varies with climate. Even so, important generalizations about potential evapo- 
transpiration still hold, stated Penman (21): 


"For complete crop covers different plants having about the 
same colour, i.e., the same reflection coefficient, the potential trans- 
piration rate the same, irrespective plant soil type.” (Physi- 
ological responses the plant, particularly high temperatures, may 
alter this 

potential transpiration rate determined prevailing 
weather.” 

transpiration short-green cover cannot exceed the 
evaporation from open water surface exposed the same weather.” 


EMPIRICAL RELATIONS 


number empirical methods have been advanced estimate actual and 
potential evapotranspiration, all employing some temperature function. Early 
work Hedke (11) produced equation whereinthe seasonal estimated con- 
sumptive use (evapotranspiration) was considered “effec- 
tive” heat expressed day degrees and taken the difference between the 
monthly temperature and germinating temperature for Using yearly 
data the same fashion, Robert Lowry and Arthur Johnson, ASCE 
(17) found high correlation between consumptive use, under where 
soil moisture was not limiting, and accumulated degree days with base 
32°F. 
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Blaney, ASCE, and Criddle, ASCE (2) correlated actual 
measurement consumptive use with the product monthly mean tempera- 
ture and the monthly day-time yearly day-time hours, 
obtain applicable coefficients for different irrigated cropping conditions 
the West. The average coefficient varies from around 0.6 for field crops 
1.2 for rice over the growth eriod under arid and semi-arid conditions. The 
coefficients, however, show marked seasonal variation and for grapefruit 
the vicinity Phoenix, Ariz., they range from 0.55 January 0.75 
August (2). 

widely used method for estimating potential evapotranspiration (P.E.) 
that devised Thornthwaite (24) through correlating tempera- 
ture with evapotranspiration determined from the water balance for valleys 
where sufficient moisture was available maintain active transpiration. The 
equation, when obtaining 30-day month values centimeters, takes the form 


which represents day-time hours expressed units hr; the 

mean monthly temperature degree centigrade; heat index obtained 


the summation the twelve monthly indices (i), which 


and the value obtained from cubic function the heat index (J), 


commenting the formula, Thornthwaite (24) stated that its work-ability 
requires the use nomograms and that the chief obstacle developing 
rational equation the lack understanding the reason why potential 
evapotranspiration corresponding given temperature not the same 
everywhere. 


DEVELOPMENT METHOD 


seeking any new relationship represent average potential evapo- 
transpiration function temperature, due consideration must given 
the geographical variation between temperature and radiation. Average 
daily temparatures July over the eastern United States vary from 80-85° 
the southern states 65-70° over the northern tier states. contrast, the 
average incoming solar radiation exhibits only small variation (6) and for 
clear skies the incoming radiation constant quantity the summer sol- 
stice for latitudes 25° 50° over North America (10). 

Since net radiation comprises the heat source area, aside from ad- 
vected heat certain localities, relationship between temperature and net 
radiation desired. Hamon (9), upon examination limited net radiation ob- 
servations and adjusting observed incoming radiation computing total long- 
wave nocturnal radiation, found that for monthly data net radiation could 
represented linear function the product day-time hours and tempera- 
ture, 
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which the average daily net radiation langleys (calories per square 
centimeter); isthe average duration units hr; and 
the average daily mean temperature degrees fahrenheit. 

The need for considering wind average estimates potential evapo- 
transpiration may reckoned from the effect wind open water evapora- 
tion and the response transpiring plants. For daily wind movement 
100 miles, the percent error evaporation with percent change wind 
the order 0.25% dry regions (16). Penman states (21), reference his 
formula (developed humid region), that possible have 100% over- 
estimate evaporation the aerodynamic term leading only 10% over- 
estimate open water evaporation. Perhaps more conclusive substantiating 
the elimination wind variable under natural conditions the response 
plants. Yamaoka (31), laboratory study transpiration from plants 
exposed constant temperature and relative humidity, found that the trans- 
piration rate remains constant under particular radiation intensity for wind 
speeds above per sec (about miles per hr). 

then appears justifiable consider wind constant factor any 
climatological estimate average potential evapotranspiration. For condi- 
tion, however, where “oasis” situation exists any estimate actual evapo- 
transpiration would highly dependent advected heat and, therefore, the 
wind. 

Adhering the concept potential evapotranspiration and accepting 
constant Bowen ratio, may relative humidity exists 
some height near the surface which would result proportionaility 
between the saturated vapor density the meansurface temperature, and the 
average difference between the saturated vapor density the surface and the 
vapor density the reference height. Therefore, other things being equal, and 
referring the turbulent transport equations, average potential evapotran- 
spiration may considered proportional the saturated water vapor 
density the air temperature near the surface, represented 


Considering that transpiration phenomenon due the function- 
ing the leaf stomata (open only indaylight) and that turbulence occurs prin- 
cipally during the daytime, adjustment for day-time hours required. 
Also, the adjustment enters again account for the disparity between net 
radiation and mean temperature, latitudinally. 

From the preceeding considerations, seems feasible formulate the 
final simplified expression for potential evapotranspiration 


(8) 


which represents the average potential evapotranspiration inches 
per day; the possible hours sunshine units hr; the satu- 
rated water vapor density (absolute humidity) the daily mean temperature 
grams per cubic meter; times 10-2; 0.55, chosen give 
yearly values potential evapotranspiration indicated observations re- 
ported the literature (24), (17). 


EVALUATION 


Yearly values potential evapotranspiration obtained from the new 
relation (Eq. for which 0.55, are practically identical those obtained 
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the highly complicated Thornthwaite method (Fig. with consistent de- 
parture about noted only for latitudes southern Florida. further 
evaluation yearly data has been made comparing average computed 
values potential evapotranspiration consumptive use data for valleys 
covering wide variation climate obtainedby Lowry and Johnson (17) and 
Williams, ASCE (30), listed Table and shown Fig. 


VALUES THORNTHWAITE METHOD, INCHES 


FIG, 1.—COMPARISON COMPUTED AVERAGE ANNUAL 
POTENTIAL EVAPOTRANSPIRATION 


important test the reliability the new relation (Eq. whether 
such simplified procedure adequate representing the seasonal 
tions. Several comparisons are made. 

Average monthly values water loss from irrigated valley, Mesilla 
Vally New Mexico, are study Lowry and Johnson (17). The 
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LINE EQUAL VALUES 


West River, Vt. 


Swift River, 


WILLIAMS 


CONSUMPTIVE USE, INCHES/YR 


Valley and State 


AVERAGE POTENTIAL EVAPOTRANSPIRATION, INCHES/ 


FIG, USE FOR BASINS COMPARED WITH 
COMPUTED POTENTIAL EVAPOTRANSPIRATION 


TABLE 1.—DATA OBSERVED WATER LOSS AND COMPUTED AVERAGE 
POTENTIAL EVAPOTRANSPIRATION FOR FIG, 


Mean Estimated 


(1) 


Mesilla, Mex, 

Pecas, Mex, 

Sangamon River, 
Green River, Ky. 
Tallapoosa River, 
Mad River, Ohio 

Skunk River, Iowa 

White River, Mo. 
Platte, 
Black River, Wisc. 
Cypress Creek, Tex. 

Wagon Wheel Gap, Colo. 
River, Colo, 


Lake Cochituate, 


114 May, 1961 
+ 
Area, in Length of 
years inches inches 
(2) (3) (4) (5) 
109,000 33.4 
5,000,000 31.4 31.4 
307,000 25.8 27.3 
1,890,000 27.0 27.6 
462,000 23.8 24,1 
222 15.6 15.0 
43,000 18,0 17.2 
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correspondence these observed values actual water loss computed 
average potential evapotranspiration satisfactory seen Fig. 
Estimates the average water loss for periods adequate water supply 


has been obtained Fox, ASCE (5) for basins Wisconsin the 

COMPUTED AVERAGE 
TRANSPIRATION 


(INCHES/ MONTH) 


MONTHS 
USE AND POTENTIAL EVAPOTRANSPIRATION 
22 (—) BIG EAU PLEINE RIVER 
(STRATFORD, WISCONSIN) 
TREMPEALEAU RIVER 
COMPUTED NORMAL POTENTIAL 
.12 EVAPOTRANSPIRATION 
WISCONSIN) 
uw 
. > 
am vee. MAR. APR. ome AUG. SEPT. am DEC. 


FIG, EVAPOTRANSPIRATION FROM WATER BALANCE FOR BASINS 


water-balance method, His determined average values water loss and the 


computed curve average potential evapotranspiration compare tavorably 
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Both Penman and Blaney (4) have computed estimates evapotranspira- 
tion their respective methods for Asheville, C., and for the climate 
this location should approximate potential evapotranspiration. addition, 
computations the Thornthwaite method and Eq. have been listed 
Table Results from the new procedure are fully satisfactory comparison 
with the other widely used methods and succeeds correcting short-coming 
the Thornthwaite method wherein winter-time estimates are considered too 
low and mid-summer values excessive (18). 

extensive investigation has been made the yearly water loss and 
lake evaporation New England and New York (12). The difference between 
the yearly lake evaporation this study, obtained adaptation the 
Penman equation (16), and values average potential evapotranspiration, 
obtained Eq. show small departures seen Fig. This close corre- 


TABLE 2,—COMPUTED AVERAGE MONTHLY EVAPOTRANSPIRATION FOR 
ASHEVILLE, C., SEVERAL METHODS 


Month 


January 
February 
March 
April 

May 

June 

July 
August 
September 
October 
November 
December 


. 
Ue 


(Totals) 


Methods: Blaney (a), Penman (b), Thornthwaite (c), and (d). 


obtained Eq. show small Fig. This close corre- 
spondence only expected inhumid regions “oasis” effect small. 

Since Eq. has shown every indicationof representing unique and simple 
procedure for obtaining useful estimates average potential evapotranspira- 
tion, computations have been made using monthly mean data for the region 
east the Rockies. map generalized average yearly values presented 
(Fig. 6). 

The conditions required for potential evapotranspiration eliminates the 
advection air heated over surrounding drier areas and allows for equilibrium 
conditions where evaporation proportional the energy Evapora- 
tion losses from pans soil tanks are always excess potential evapotran- 
spiration—ranging from differences humid regions excesses 
hot-arid regions. For large irrigated areas valleys this discrepancy 
practically eliminated. Under non-irrigated conditions, the actual evapo- 
transpiration depends upon the available soil moisture andthe potential evapo- 
transpiration may adjusted for its estimation. 

The decrease evapotranspiration relation soil moisture viewed 
differently various investigators. Thornthwaite (18) has adopted system 
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that decreases evapotranspiration direct function the remaining soil 
moisture the basis field showed the fraction available energy 
used evaporation proportional soil moisture. From data 
Wadleigh (28), the suppression plant growth andthe soil moisture above the 
permanent wilting point are inthe same fashion with soil moisture 
This indicates that the adjustment potential evapotranspiration 
obtain actual evapotranspiration should made proportional the plant 
available moisture under conditions where soil moisture limiting. Whether 
not soil moisture “limiting” depends upon the soil and soil aeration, 
the plant, and meteorological factors. 


CONC LUSIONS 


The equation developed estimate average potential evapotranspiration 
(Eq. adheres the desirable features requiring only readily available 
mean temperature data-converted saturation water vapor density and day- 
time hours, expressed easily computable fashion. General applicability 
seems justified view the close correspondence between observed and 
computed values average potential evapotranspiration, both ona yearly and 
seasonal basis for widely scattered localities. 
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SYNOPSIS 


and transverse spacings. 


the roughness elements. 


The experimental results are compared with those some other investi- 
gators. addition, some observations are made onthe von Karman turbulence 
coefficient and the problem depth determination rough channels. 


INTRODUCTION 


Notation.—The letter symbols adopted for use this paper are defined 
where they first appear, the illustrations the text, and are arranged 


alphabetically, for convenience reference, the Appendix. 
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series experiments determine the effect roughness spacing 
open channel flow reported. The experiments were conducted 8-ft 
wide adjustable-slope flume. Roughness elements consisting 
sheet metal baffles were placed the bedof the flume various longitudinal 


The data are analyzed terms the von Karman-Prandtl concepts 
turbulent flow near boundary. Logarithmic flowformulas are developed 
which the roughness factor size and spacing characteristics 
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Many important phenomena involved the flow water open channels 
are only partially understood. One the most basic these the effect 
boundary roughness the flow. 

Over the years, number empirical roughness formulas have been 
developed. Almost all are based the Chezy formula 


which the mean velocity flow, denotes the Chezy resistance func- 
tion, defines the hydraulic radius the channel cross section, and the 
slope energy gradient. Formulas developed Manning, Ganguillet and 
Kutter, Bazin, and others relate roughness. These formulas 
have little rational foundation. They give fair results when applied over 
the fairly narrow range flow conditions which they were based, but fre- 
quently lead results that are seriously inerror when applied outside this 
range. Experiments Sayre3, ASCE, and Sayre and 
show conclusively that the roughness coefficients formulas 
not necessarily remain constant for particular boundary roughness. 

the last few decades, the development the boundary layer theory has 
opened the way toward rational understanding some the flow phenomena 
involving boundary roughness. 

Using the von Karman-Prandtl concepts turbulent flow near rough 
boundary addition data from his own experiments with pipes having sand- 
coated interiors, Nikuradse® developed artificial roughness standard for 
pipes, that based the diameter the sand grains constituting the 
roughness, 

The application and others the Nikuradse roughness 
standard open channels has been quite successful describing grain-type 
roughness wide, open channels. However, has been found inadequate for 
describing certain other types roughness which the relative spacing 
addition the relative size the roughness elements important bound- 
ary characteristic. Dune and ripple formations onthe beds alluvial channels 
are considered typical the latter type roughness. 

Using square strips extending across the bed the channel roughness, 
Ralph ASCE, studied the effect the longitudinal spacing 
the strips. and Robinson and reported experiments 
which the size geometrically similar roughness baffles was varied, but the 
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ratios longitudinal and transverse spacings baffle height were held con- 
stant. For particular roughness pattern they demonstrated that the Chezy 
resistance function depends only the relative roughness (ratio flow depth 
baffle height), assuming rough boundary conditions. 

These and other studies have contributed toward the establishment 
reproducible artificial roughness standard that will applicable open 
channels characterized “dune”-type roughness. Itis hoped that this study 
will stimulate further interest toward this end. 


ANALYSIS 


Theoretical Analysis.—The von Karman-Prandtl equation for velocity dis- 
tribution near rough boundaries the form 


2.30 


which the velocity distance from the boundary, defines the von 
Karman turbulence coefficient, the distance above the bed the channel, 
denotes the measure roughness size, defines constant integration, 
denotes the shear velocity theboundary. Its derivation appears 
many standard fluid mechanics texts. has been verified experiment 
over wide range conditions pipes, fluids moving past flat plate, and 
open channel flow provided that secondary circulation not important 
consideration. 

expression for the mean velocity can obtained integrating Eq. 
over the depth flow the following manner 


V Yn (3) 


that results the expression 


(4) 

K log ee eevee ee 


which the normal depth (the depth flow occurring when slopes 
energy gradient, water surface, and bed are equal), and defines experi- 
mental constant dependent the roughness spacing. From the Chezy and 
Darcy-Weisbach equations 


which the Chezy coefficient, denotes the hydraulic radius, defines 
the slope the energy gradient and the Parcy-Weisbach resistance co- 
efficient. Eq. can thus written the form genéral resistance equa- 
tion such 


2.30 
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0,81 


Equations proposed Keulegan, Powell, and Robinson and Albertson are 
the form Eq. 6(a). The von Karman-Prandtl equation for rough pipes 
the form Eq. 6(b). 

conclude that must factor dependent upon the arrangement and shape 
the roughness alternate form Eq. 6(a) would 


which appears roughness parameter the size, shape 
and spacing the roughness elements. other words, should completely 
describe the boundary roughness. 

Two fundamental assumptions are the foregoing analysis. The 
first that the boundary hydrodynamically rough (roughness protruding 
above the laminar sublayer) and that consequently viscosity has negligible 
effect. The second that the channel sufficient width, that the bed 
roughness great relative the side-wall roughness both, elimi- 
nate any appreciable sidewall effect. 

Dimensional variables used the dimensional analysis 
were chosen represent the experimental conditions. They can 
grouped under three headings. 

Boundary Characteristics.—The boundary characteristics are follows: 


width channel; 

height roughness baffle; 

width roughness baffle; 

longitudinal spacing roughness baffles; 

transverse spacing between roughness baffles; 

thickness roughness baffle; 

shape factor roughness elements; 

roughness parameter describing size, spacing and shape rough- 
ness elements, and having the dimension length; and 

shape channel. 


Flow Characteristics.—The flow characteristics are: 


mean velocity flow; 
hydraulic radius; and 
boundary shear. 


Fluid Characteristics.—The following comprise the fluid characteristics; 


mass density fluid; 
dynamic viscosity fluid; and 
difference specific weight water and air the fluid interface. 


completely defines the roughness the theoretical analy- 
sis, the relationship existing among the variables describing the roughness 
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can stated 


The general relationship among the remaining variables 


Choosing repeating variables, the Buckingham Pi-theorem 
yields -terms among which the following relationship exists: 


y/p 


R’ R’ ’ 
wide channel, canbe substitutedfor and there are side-wall 
effects B/R can eliminated. Because the shape the channel remains 


rectangular and wide throughout the experiments, can also eliminated. 
Eq. can now restated 


which and are the Reynolds and Froude numbers respectively. 
use 


Eq. can simplified and rewritten 


the height the roughness elements large relative the thickness 
the laminar sub-layer, viscous effects are negligible, and can eliminated. 
The Froude number assumes importance only surface waves 
disturbances are present. Consequently, boundary hydrodynam- 
ically rough and whenthere are surface eaves, Eq. reduces the simple 


relationship 


Turning now the problem roughness geometry, Eq. stated that 


(7) 


The ratio b/a and the shape factor the roughness were held con- 
stant and, therefore, canbe Becausetis very small (approximately 


16 
7 
> 
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0.05 relative all other roughness dimensions, t/a can considered 
effect equal zero, and hence negligible importance. Because the dimen- 
sions the individual roughness baffles remain constant the ratio X/a 
independent baffle size, andcanbe replacedby parameter that appears 
dependent only the baffle spacing. Eq. can now written sim- 
pler form 


reasoning, the roughness size and spacing dimensions can combined 
into single dimensionless density parameter b/x(e+b). This parameter 
numerically equal the ratio between the all the roughness 
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FIG, 1.—FLUME AND CIRCULATION SYSTEM 


elements projected perpendicularly the direction flow the numerator 
and the total area the flume bed the denominator, 


x(e+b area flume bed 


Another expression for the spacing parameter then 


DESCRIPTION EXPERIMENTS 


The experiments were conducted 8-ft wide 2-ft deep rectangular 
flume which measured length. The flume and circulation system are 
shown schematic form Fig. 
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The bed and walls the flume were constructed 1/2-in. plywood rein- 
forced 2-in. ribs and cross members. The flume rested pair 
6-in. I-beams which were supported screw jacks 10-ft intervals. 
Changes the longitudinal slope the flume could accomplished ad- 
justing the jacks. 

the flume entrance flow was distributed and excessive local turbulence 

The water-surface level the flume could controlled adjusting the 
vertical-slat tailgate the downstream end the flume. The tailbox, 
single-stage propeller pump, and pipeline completed the major 
elements the circulation system. 

The discharge was controlled valve the line. Discharge was 
measured means orifice plate located the 14-in. pipe, 
downstream from the pump. 

Water-surface slope and flow depth were determined five stationary 
point gages which were spaced over length along the flume. The 
point gages were read with the aid engineer’s level. 11/32-in dia- 
meter Prandtl-type pitot tube with tilting manometer was used obtain 
velocity profiles the vertical. The water temperature was measured with 
ordinary laboratory mercury thermometer. 

Sheet metal baffles measuring 6-in wide and 1/2-in high were used 
roughness elements. shown Fig. the baffles were placed definite 
longitudinal and transverse intervals, thus forming symmetrical repeating 
patterns. Each baffle was centered the space between baffles the rows 
immediately upstream and downstream, hence alternate rows were the same. 
All patterns were symmetrical about the centerline the flume. Half-width 
baffles were placed adjacent the side walls order maintain the same 
geometry there along the centerline. The photograph Fig. shows one 
the roughness patterns. All the patterns are shown schematic form 
Fig. 

performing the experiments, the discharge, slope, and roughness pattern 
were the independent variables, and the normal flow depth was the primary 
dependent variable. Once predetermined combination independent variables 
was established, the procedure during any given run consisted essentially 
manipulating the tail-gate setting until the water-surface slope equaled the 
slope the flume bed. This could done quite accurately intentionally 
establishing and then interpolating between slight and backwater curves. 
The condition equal slopes the water surface and the bed defined normal 
flow. The depth normal flow was then measured, and velocity profiles were 
obtained least three verticals across the flume. 

Tests were conducted discharges varying from cfs cfs and 
slopes 0.001 and 0.003 for each the six roughness patterns. total 
twenty-seven runs was made. 


AND DISCUSSION RESULTS 


Data describing the overall flow characteristics and the roughness are listed 
Table The relationships developed the theoretical and dimensional 
analysis are largely substantiated the experimental results the form 


Engrg., Colorado State Univ, Research Foundation, October, 1958, pp. 
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dimensionless graphs that are presented and discussed the following sec- 
tions. addition, information pertaining the von Karman turbulence coef- 
ficient, and the problem depth determination rough channels also 


presented. 
Chezy Coefficient Function Relative Roughness.—Fig. plot 
the Chezy resistance function against relative roughness based the baffle 
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FIG, 2.—DEFINITION SKETCH ROUGHNESS BAFFLE INSTALLATION 


height. The data fall into family parallel lines when plotted semi- 
logarithmic paper. Each line represents data obtained from experiments with 
particular roughness pattern. Robinson’s datal0 which are superimposed 
the plot, represent experiments which the same roughness pattern, but two 
different baffle sizes, were used. The roughness baffles used Robinson 
were geometrically similar those used the writers. 
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Data obtained H.J. ASCE, for conditions stable flow 
are also plotted Fig. These tests were conducted 2-ft-wide flume 
with 3/16-in. brass cubes roughness elements. Itis interesting note that 
Koloseus’ data were obtained the super-critical (F>1) flow range, whereas 
Robinson’s and the writers’ data all represent flow conditions well within the 
sub-critical range. The data are all described the equation 


which similar Eq. 6(a) given the theoretical The constant 
for each line, that each roughness pattern, defined extrapolation 
the value which The writers’ and Koloseus’ data show 


3.—EXAMPLE ROUGHNESS PATTERN 


clearly that function the longitudinal and transverse roughness 
spacing, and Robinson’s data show that independent absolute roughness 
size. 

Fig. plot the Chezy resistance function against the relative rough- 
ness based the roughness parameterX Description the boundary rough- 
ness the single parameter groups all the data about single curve 


«The Effect Free-Surface Instability Channel Resistance,” Koloseus, 
ment the requirements for the degree Doctor Philosophy. 
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FIG, DIAGRAM THE SIX ROUGHNESS PATTERNS 
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FIG, RESISTANCE FUNCTION WITH RELATIVE ROUGHNESS 


Robinson 


ROUGHNESS SPACING 
described the equation 


— = 6.06 1 6 
(18) 


that similar Eq. 6(c) the theoretical analysis. The parameter has 
dimension length and particular boundary roughness. The 
baffle size and spacing. The properties and Eq. are combined 
according the relationship 


Cp = - 6.06 log (19) 


The ratio X/a was suggested spacing parameter the dimension- 
analysis. 

Eqs. and are essentially equivalent. However, definition the relative 
Grouping all the data about single curve the most obvious advantage. 
Description the roughness single term greatly facilitates the inclusion 
the roughness effect various dimensionless parameters. addition, 
Eq. possesses the same advantages simplicity application does the 
Manning formula which the single term ndescribes the boundary roughness. 

Should demonstrated, however, that the coefficient 6.06 varies appreci- 
ably with roughness form channel shape but not with roughness spacing. 
Eq. would seem more significant than Eq. 18. 

The empirically-determined coefficient 6.06 Eqs. 17and18 equivalent 
von turbulence coefficient 0.38. This compares favorably with 
the coefficient 0.40 determined Nikuradse® for turbulent flow smooth and 
rough pipes, and the coefficients determined number investigators 
open channels. 

Figs. and some the points representing observed flow conditions 
for the two sparsest the writers’ roughness patterns are seen fall below 
the curve established the other points. Inthese runs combination small 
depths and relatively high velocities occurring directly over the crests the 
roughness baffles caused the formation breaking standing waves immediately 
downstream from the baffles. effect, these baffle-generated waves con- 
stituted energy sinks not accounted for either Eq. Eq. 18. This ex- 
plains the low values obtained for these runs, and also agrees with 
Koloseus’ for unstable flow which clearly demonstrate that channel 
resistance increased systematic manner free-surface instability. 

The tendency Robinson’s data drift below the curve larger values 
relative roughness attributed increase sidewall effect accompany- 
ing the increase depth. Robinson was only in. wide, and 
the depth flow several runs approached the same magnitude. 

Comparison Chezy and Manning Formulas.—Fig. compares 
with the Manning formula. Basically the same plot Fig. except that 
extended over two additional log cycles relative The straight 
dashed lines are extrapolations the experimental data, each line represent- 
ing particular degree boundary roughness. The solid curved lines repre- 
sent constant values Manning’s thatwere plotted onthe versus 


coordinates bysetting equal The line compares Eq. 
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FIG, 6.—VARIATION RESISTANCE FUNCTION WITH RELATIVE ROUGHNESS 
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with the Manning formula roughness commonly encountered 
practice. 

Fig. demonstrates conclusively that application the Manning formula 
the range roughness included the experiments would lead serious 
error. seen that the boundary described the baffle spacing ratios 
x/a and e/a 3.9 wouldhave values Manning’s ranging from approxi- 
mately 0.060 0.100 over the relatively small two-fold range depth covered 
the experiments. However, the region where the Manning curves are 
approximately tangent lines having slope log 6.06, the Manning formula 
and Eq. should give results which are fairly close accord. This the 
range conditions within which the Manning formula commonly applied 
open channel flow. 

Influence Roughness Pattern.— 

Resistance Flow.—Fig. plotofthe roughness spacing 
against the roughness density ratio that was suggested the dimen- 
sional analysis. Data from several additional sources are superimposed the 
plot. Robinson’s and Koloseus’ data have been the preceding sec- 
tion. Schlichting’s datal3 were obtained closed rectangular conduit 
measuring 4-cm high 17-cm wide. Baffles measuring 0.3-cm high 0.8-cm 
wide were fixed the top plate. 1946 roughness consisted 
and 1/8-in square strips extending continuously across the bottom and the 
sides 8-in-wide flume. Rand’s datal4 were obtained 2-ft-wide 
flume which 1/4-in. square extending across the bottom were used 
roughness. Bazin’s experiments were conducted 2-m-wide flume 
which strips 1-cm thick and 2.7-cm wide extended across the bottom and 
the sides. the Waterways Experiment Station triangular 
ridges 0.04-ft high with vertical upstream face and the downstream face 
sloped served roughness. The ridges extended across the bottom 
and the sides the 2.08-ft-wide flume. The brushed-concrete surface upon 
which the ridges were superimposed was rougher thanthe base surface the 
other experiments. 

Considering the differences roughness and channel forms seems that the 
various curves should compared only qualitative sense. There 
general tendency Fig. 8for the spacing parameter X/a), that propor- 
tional the hydraulic roughness factor roughness height, increase 
with increasing roughness density. For the continuous strip roughness 


becomes smoother the roughness density increases. For the triangular 
ridges, maximum roughness occurs about 0.2. Evidently, for com- 
parable spacings, the ridges are somewhat rougher than the square strips, 
particularly close spacings. the more distant spacings the comparitively 
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rough brushed-concrete surface that the triangular ridges 
evidently came into play. 

For the continuous roughness, interesting note that the roughness 
density parameter reduces a/x, because the transverse spacing between 
roughness elements zero. seems likely that the curves representing the 


/ 


LEGEND 
Sayre— 
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Koloseus— cubes, 
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FIG, 8.—VARIATION SPACING PARAMETER WITH ROUGHNESS DENSITY 


isolated roughness elements (baffles and cubes) also should reach maximum 
values but greater roughness densities that extend beyond the range 
the data. 
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would expected, the rising limbs the curves Fig. are arranged 
order increasing relative width (b/a) the roughness elements the 
spacing parameter increases. 

plot somewhat similar Fig. which contains more complete data for 
continuous-strip roughness has been given Joe Johnson, ASCE, and 


FIG, VELOCITY PROFILES 


Roughness,” Rowell, Transactions, ASCE, Vol. 11, 1946, pp. 555-559. 
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Velocity Distribution.—Profiles velocity distribution the vertical are 
plotted dimensionless form Fig. Each point represents weighted 
average least three readings taken the same elevation above the bed, 
but different verticals across the width the flume. The weights were 
assigned according the position the vertical with respect the rough- 
ness elements. 

The nature the velocity profiles clearly seems have been affected 
roughness pattern, particularly the longitudinal spacing the baffles. 
Velocity profiles for the four denser roughnesspatterns are quite well repre- 
sented the equation 


for which the constants were determined empirically. 
For the sake comparison the von equation 


for velocity distribution near rough boundary plotted 
broken line. The Nikuradse sand-grain roughness related the non- 
linear relationship 

12,2 


0.053 


that corresponds approximately over the range conditions 
tested. 

The deviation the points representing velocity readings taken for the two 
sparsest roughness patterns attributed the non-uniformity the rough- 
ness and hence the turbulence the sparse-pattern experiments. The raw 
velocity-profile data are 

General Resistance Diagrams for Open and Fig. were 
plotted for the purpose determining the order and range hydrodynamic 
roughness covered the experiments. Some available open-channel data rep- 
resenting other regions the roughness spectrum are included these plots 
means testing the applicability the Prandtl-von 
type resistance diagrams, that were developed for flow pipes, flow 
wide, rigid-boundary open channels. Assumption wide-channel conditions, 
that justified for all data shown Fig. and Fig. with the 
exception the smooth-boundary data which case was substituted for yp. 

Fig. 10, which open-channel adaptation the well-known general 
resistance diagram for flow pipes, illustrates the effect the Reynolds 
number, that viscosity, the resistance function. 

The laminar-flow equation 


1/2 
— = 0.577R ee (23) 
ve 


derived directly from the following equations describing steady, uniform, 
laminar flow 


dx dy (24a) 
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The smooth-boundary formula 


direct adaptation Blasius’ empirical relationship for flow smooth 
pipes. The other smooth-boundary equation 


smooth channel infinite width. The coefficient changed 
6.06 for the sake conformance with Eq. and Eq. 18, and the constant 
was changed compensate for the alteration the coefficient. 

Experimental data for the six roughness patterns fall well within the 


rough-boundary region that begins this condition the ratio 


the roughness the thickness the laminar sublayer one (see Fig. 
11). The data are described afamily horizontal lines, each having con- 
stant value relative roughness Hence, the region covered the 
data, the resistance function dependent only relative roughness and 
therefore independent the Reynolds number viscous effects. 

The laminar-flow and smooth-boundary data shown the plot were obtained 
from variety Laminar-flow data were obtained Owen18 
rectangular, 1.50-ft-wide brass-plate flume, and the Waterways Ex- 
periment rectangular, 2.30-ft-wide flume smooth cement 
finish. Smooth-boundary data were obtained Tracy, ASCE, and 
the Waterways Experiment the previously described cement 
flume. The transition from laminar turbulent flowfor the smooth-boundary 
data occurs the range 1,200. 

Fig. open-channel adaptation the transition function for pipes. 


means the abscissa that proportional the ratio the data 
n/X 


are arranged order hydrodynamic roughness. 


actions, ASCE, Vol, 119, 1954, 1157. 

“Study River Bed Materials and their Movement with Special Reference the 
Lower Mississippi River,” Waterways Experiment Sta., Corps Engrs., Vicksburg, 
Tech. Memorandum No, 17, 1935, 161 pp. 

Lester, Administrative Report, Survey, 1959. 
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The smooth-boundary equation 


obtained merely subtracting 6.06 log from both sides Eq. 26. 
The transition equation 


(29) 


1.58 


was obtained combining the smooth and rough-boundary equations the 
manner proposed Colebrook and White2! for pipes, and asymptotic 
both the smooth-boundary and rough-boundary equations. The broken-line 
curves, indicating the extremes ofthe transition region, approximate the posi- 
tions the uniform and non-uniform sand roughness curves the trnasition 
plot for pipes. 

Data for the six roughness patterns appear the rough end the scale. 
Although most the data representing Koloseus’ two patterns also fall the 
rough-boundary region, few points representing the sparser the patterns 
lap over slightly into the transition region. The Waterways Experiment Station 
datal9 for beds consisting gravel, 0.525 sand, 0.310 sand 
and smooth cement, and the steel-flume data Tracy and are also 
plotted Fig. 11. The alluvial-bed data all represent flow conditions before 
the beginning movement, that inallcases the boundary effect rigid. 
plotting the smooth-boundary data, was assumed that 0.00004 for 
the smooth-cement finish, and that 0.00001 for the steel surface. Eq. 
was used computing from the median diameters the 0.525-mm 
sand 0.00019 ft) and the 0.310-mm sand 0.00011 ft). The 4-mm gravel 
data fall raostly the rough-boundary region, overlapping into the transition 
region. The data for the 0.525-mm sand fall the transition region, for the 
sand the transition and smooth regions, and for the smooth- 
cement ana steel-flume surfaces, along the smooth-boundary curve. 

Scattering the data the smooth-boundary and transition regions not 
6.06 log direct difference rather than percentile deviation index. 
For example, considering hypothetical point located 6.06 log 


smooth-boundary curve one unit 25%. However, based upon C/s/g value 
15, which about average for the data the smooth and transition regions, 
deviation one unit amounts less than 7%. 

The ratios and are plotted along the abscissa scale Fig. 11. 
The ratio average over the approximate range relative 


and 0.320, the deviation the ordinate scale from the 


“Turbulent Flow Pipes, with Particular Reference the Transition Region Be- 
tween the Smooth and Rough Pipe Laws,” Colebrook, Journal, Engrs., 
11, 1939, pp. 
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roughness 10,000 included the data plotted Fig. 11. The fact 
that the rough-boundary region begins approximately appears 
significant. This compares with indicated22 Hunter Rouse 
ASCE, the approximate demarcation between the transition and rough- 
boundary regions for pipes. 

The wide band within which data may fall the can 
explained the basis roughness spacing size distribution both. 
Where the spacing between roughness elements large the size rough- 
ness elements not uniform, isolated elements even may project 
above the laminar sublayer sufficient numbers upset the turbulence 
pattern generated the smooth boundary, but not sufficient numbers 
cause rough-boundary turbulence. However, when the roughness elements are 
densely-packed and uniform size, smooth-boundary flow may persist 
approaching one and the transition from smooth rough boundary con- 
ditions would much more abrupt. 

The von Karman Coefficient.—One the purposes the study was throw 
some additional light the question how the von Karman might vary, 
fact does vary, with boundary roughness suggested Rand.23 

The coefficient initially appeared the literature “universal con- 
stant” turbulent exchange von Karman’s hypothesis24 turbulent simili- 
tude. the basis Nikuradse’s work® with smooth and rough pipes, the 
value has commonly been accepted equal 0.40, although Vito 
Vanoni25 ASCE, analyzing Nikuradse’s ranging from 
0.32 0.42 with 0.37 average. Evidence has been mounting 
steadily show that expressed the logarithmic velocity distribution 
equations, not constant, but systematically influenced such factors 
suspended-sediment concentration, secondary circulation, and perhaps 
boundary roughness. 

the experiments reported herein two methods determining were 
used. the first method was determined from analysis the velocity- 
profile data through use the relationship. 


mic cycle depth. The portions the velocity profiles that followed most 
closely logarithmic distribution were used determining and velocity 
measurements either the baffle crests water surface were given 
only secondary consideration. 

The second method determing was based Eq. from which 
seen that the slope each the lines drawn Fig. defined the quantity 


“Elementary Mechanics Fluids,” Rouse, John Wiley and Sons, Inc., New 
York, 1946, 376 pp. 


Discussion Rand “Wind-Tunnel Studies Sediment Movement,” Proceed- 


ings, Fifth Conf., State Univ. Bulletin 34, 1953, pp. 132-133, 
Ahnlichkeit und Turbulenz,” Theodor von Karman, Proceedings, 
Third Internatl. Congress for Applied Mechanics, Stockholm, Sweden, 1930. 

“Some Effects Suspended Sediment Flow Characteristics,” Vanoni, 
Proceedings, Fifth Hydr. Conf., State Univ. Studies Engrg., Bulletin 34, 1953, 
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Best-fit lines were drawn through the data representing each rough- 
ness pattern Fig. and average integrated for each rough- 
ness pattern was obtained from the relationship 


which the change the resistance function over one logarithmic 
cycle relative roughness. 

determined both methods are shown this plot. The curves compare the 
average determined each method, each point the curve 
representing average value determinedfrom all the runs for particu- 
lar roughness pattern. 

Fig. indicates that the roughness pattern instrumental causing the 
apparent variation The principal significance Fig. 12, however, lies 
the comparison between the average -values determined the 
velocity-profile method and the integrated method. Agreement remarkably 
good for the three denser patterns, but for the two sparsest patterns there 
little semblance agreement. probable that the integrated method pro- 
vides more reliable index than does the velocity-profile method, and 
that the large apparent -values obtained the velocity-profile method for 
the sparse patterns are indicative large-scale local vortices set the 
isolated roughness baffles, rather than such large variations This 
expected, however, because the theory was originally developed for 
cases wake-interference flow pipes, where the turbulence, close the 
boundary, essentially random nature. Extension this theory the 
isolated-roughness case, where the flow three dimensional the rough- 
ness non-uniform, both, is, therefore, neither logical nor justified. 

This being the case, the portion Fig. drawn solid line appears 
justify the adoption average value for the range conditions tested. 
effect, selection the constant 6.06 the coefficient Eqs. and 
amounted adopting the average value 0.38. 

cannot said that the experimental results provided entirely satis- 
factory answer whether the von Karman fact constant for clear 
water open-channel flow, how and why might vary not con- 
stant. For the type roughness studied, the velocity distribution method 
determining proved deficient respects for the sparse patterns. 
Without taking large number velocity profiles verticals various 
positions relative the roughness baffles, was impossible ascertain 
whether not truly representative picture the velocity distribution had 
been obtained. Limitations the equipment and deviations 
from the logarithmic law close the baffles and near the water surface added 
the complications. Nevertheless, based onthe experimental results, indica- 
tions are that, provided the conditions essentially wake-interference flow 
and uniform -roughness are met, constant for steady, uniform flow 
clear water. 

Depth Determination Rough Channels.—Various investigators, such as, 
have suggested that when roughness elements are closely spaced, the 


Discussion Rand “Artificial Roughness Standards for Open Channels,” 
Transactions, Amer, Geophysical Union, Vol. 35, No. August, 1954, 649-650. 
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depth flow should measured from the water surface some reference 
level between the crests and the troughs the roughness elements. 

implied the section describing the experiments, the depth flow 
measured from the water surface the flume bed was used the analysis. 
However, quite extensive study was made the velocity profiles obtained 
for the three densest roughness patterns before was decided not apply 
any depth corrections. 

detailed description the depth-correction analysis beyond the scope 
this paper. The negative results this analysis appear justified for quite 
dense roughness patterns, however, when the roughness elements consist 
thin isolated baffles such experiments. This type roughness 
geometry permits flow through the transverse spaces between roughness ele- 
ments. addition, because the proportion the total bed area actually taken 

the roughness elements— 
extremely small, the turbulence down close the bed level the flume 
probably not greatly inhibited. The results turbulence flume studies 
reinforce these hypotheses. these studies staggered pattern 
individual roughness baffles proved extremely effective maintaining 
large sediment concentrations suspension without appreciable deposit, 
whereas extensive deposits occurred comparable concentrations when 
the roughness consisted baffles extending continuously across the width 
the flume. These studies have application here because the basic 
similarity between the mechanisms mass transfer andmomentum transfer. 


0.008 for the densest pattern—is 


CONCLUSIONS 


Analysis the experimental data from this and other studies gives rise 
the following conclusions with respect evaluating the roughness rigid- 
boundary, rough open channels: 


The experimental results give additional evidence that the variation 
the Chezy resistance function with the relative roughness logarithmic 
nature. 

The equations 


are found apply over the range roughness conditions encountered the 
six arrangements roughness baffles tested. these equations, the constant 
depends only relative spacing roughness elements height 
andthe general roughness parameter function both the relative size 
and relative spacing the baffles. 

Eq. shown considerably more accurate than the Manning 
formula over the range roughness and flow tested. 


“Turbulence Flume Measure Bed Load,” Albertson, Transactions, 
Amer, Geophysical Union, Vol. 32, No. December, 1951, pp. 
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The ratio which equivalent the combined area all rough- 


ness baffles projected perpendicularly the direction flow divided the 
total area the bed, found adequate definition the roughness 
density. 

Distribution velocity the vertical direction when the roughness 
pattern relatively dense found similar that described the von 
equations for rough boundaries. The equation 


adequately describes the velocity distribution such pattern con- 
sisting sparsely-spaced, isolated roughness elements, however, gives rise 
persistent horizontal and vertical vortices, and consequently considerable 
local variations the turbulence pattern which cause marked deviations from 
the von type distribution. 

The general resistance diagram, which the resistance function 
plotted against the Reynolds number, and the Colebrook-White type transition 
function are found applicable uniform flow wide, rigid- 
boundary open channels. 

Indications are that the von turbulence coefficient inde- 
pendent roughness pattern, provided that conditions wake-interference 
flow and uniform roughness prevail. the final analysis effect con- 
sidered constant, equal 0.38, for the range conditions tested this 
study. 

Negative depth corrections are found unnecessary for roughness 
systems consisting staggered individual baffles, least when the trans- 
verse spacing ratio e/b about and the volume small, 
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prepared committee the American Standards Association with Society 
representation, and approved the Association 1942: 


Definition Dimensions 
Height roughness baffle. 
Width flume. 
Width roughness baffle. 
Chezy discharge coefficient. 
Constant integration. 


Experimental constant dependent rough- 
ness spacing. 


Transverse spacing between roughness baffles. 
Froude number. 

Darcy-Weisbach resistance coefficient. 
Acceleration gravity. 


Measure roughness size. 


Dn 


Equivalent sand roughness (Nikuradse). 
Roughness coefficient (Manning). 

Flow discharge. 

Hydraulic radius. 

Reynolds number. 

Slope energy gradient. 

Thickness roughness baffle. 

Mean velocity flow. 

Velocity distance from boundary. 
Longitudinal spacing roughness baffles. 


Distance above bed channel. 


Normal depth, that is, depth flow occurring 
when slopes energy gradient, water surface, 
and bed are equal. 


involved (water and air) fluid interface. 


Thickness laminar sublayer. 
11.6v 


y70/P 


von Karman turbulence coefficient. 


Symbol 
a = 
e 
L1/6 
= 
L 
L/T 
L/T 
4 = 
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Channel shape factor. 

Dynamic viscosity fluid. 

Kinematic viscosity fluid. 

Fluid mass density. 

Shape factor roughness elements. 
Shear intensity. 

Shear intensity acting bed channel. 
Roughness parameter. 

Spacing parameter. 


FT/L2 
F/L2 
F/L2 


r = 

= 
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EDDY DIFFUSION RESERVOIRS AND PIPELINES 


SYNOPSIS 


Previous studies have shown that the ocean eddy diffusion coefficient 
function eddy size. Reservoir tests with radioactive tracers confirm this. 
Tests eddy diffusion pipelines are analyzed and compared with Taylor’s 
theoretical solution. The pipeline coefficients are greater, expected, be- 
cause obstructions and curved sections the pipe. 


INTRODUCTION 


Diffusion ponds and pipelines combination molecular and eddy 
Molecular diffusion the result the heat motion in- 
dividual molecules and extremely slow liquids and gases, because the 
mean free path molecules short. Eddy turbulent diffusion, while 
fundamentally the same process molecular diffusion, the mass transfer 


Note.—Discussion open until October extend the closing date one month, 
written request must filed with the Executive Secretary, This paper part 
the copyrighted Journal the Hydraulics Division, Proceedings the American So- 
ciety Civil Engineers, Vol, 87, No. May, 1961. 

Engrg. Leader, Oak Ridge Lab., Union Carbide Corp., Atomic Energy 
Comm,, Oak Ridge, Tenn, 

“Modern Developments Fluid Dynamics,” Goldstein, Clarendon Press, Ox- 
ford, England, Chapter 1938, 

“Turbulent Diffusion,” Batchelor, Lecture Series, No. Univ. Mary- 
land, College Park, 

«Introduction Some Topics Turbulence,” Frenkiel, Lecture Series, 
No. Univ. Maryland, College Park, 1960. 

“Statistical Theory Turbulence,” Goldstein, Lecture Series, No. Univ. 
Maryland, College Park, 1960, 
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molecules eddies various sizes. Essentially, molecular diffusion 
submicroscopic movement, while turbulent diffusion macroscopic. Eddy 
diffusion has been found vary from 1,000 1,000,000 times greater than 
molecular diffusion open bodies Therefore, molecular diffusion 
will ignored the computations. 

Turbulent diffusion may taking place numerous ways within the same 
area. Large-scale eddies and small-scale eddies within them induce various 
rates mixing. The complete randomness the eddy diffusion process makes 
difficult phenomena formulate, and only recently with the applica- 
tion statistical methods the problem that some ideas the fundamental 
nature the process are being uncovered. present, (1961), only pos- 
sible approximate the magnitude the diffusion rate. 


HISTORY DIFFUSION STUDIES 


The earliest work turbulent diffusion was that Osborne Reynolds who 
differentiated turbulent flow between the primary character- 
istics the motion. The separation turbulent velocity into two components 
accomplished assuming the mean velocity along the direction travel 
the primary component and superimposed turbulent velocity having 
temporal mean velocity zero the secondary component. This can car- 
ried extreme taking the time and length scales small enough that 
even laminar flow would appear turbulent. 

Prandtl advanced turbulence theory relating the intensity turbulent 
flow parameter distance. The “mixture length” Prandtl represents 
the average distance small mass fluid will travel before loses its 
increment momentum the region into which The work 
Prandtl, however, bore close relationship the physical nature turbu- 
lence. Taylor next applied the concept the transfer vorticity rather 
than momentum, but still retained the concept mixture length which repre- 
sents the transverse path the vorticity transfer. This formulation accorded 
better with experimental results. Theodor von Hon, ASCE, at- 
tempted improve these efforts with his thesis turbulent similitude based 
the assumptions that the mechanism turbulence independent visco- 
ity, except the immediate neighborhood the flow boundaries and that the 
pattern secondary flow statistically similar from point point, varying 
only time and length scales. These assumptions lead equations the 
same form Prandtl’s and Taylor’s but withthe mixture length function 
the velocity distribution the mean flow. difference between momentum 
transfer and vorticity transfer that momentum transfer the eddy coeffi- 
cient function mixing length and the transverse velocity the stream, 
while vorticity transfer the eddy coefficient independent the secondary 
velocity. 

These methods lead relatively simple and semi-empirical solutions 
which are widely used engineering applications. For more adequate theo- 
retical treatment and closer representation the physical reality; they are 


New York, 1942, 91. 

Mechanics for Hydraulic Engineers,” Rouse, McGraw-Hill Book Co., 
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not sufficiently accurate enable fundamental research into turbulent flow 
proceed. The fundamental problems turbulence are best approached 
through statistical analysis. Taylor produced the first study turbulence 
statistical methods. The notions introduced Taylor the investigation 
isotropic turbulence were the study the correlation between two fluctuating 
quantities the turbulent flow, the study the decay turbulence, and the 
study the spectrum turbulence (that is,the distribution energy among 
eddies different sizes). The foregoing devices are means solve the main 
question the joint probability distribution the displacement any “n” 
marked particles. 

The theoretical work the field continues ever-eccelerating rate 
such physicists and applied mathematicians Heisenberg, Chandrasekhar, 
Onsager, Von Weizacker, von Batchelor, Townsend, Dryden, Gold- 
stein, Feriet, Loitsiansky, Taylor, and Kolmogoroff. This research has 
been mainly concerned with laminar flow the boundary layer, the initial 
stages transition from laminar turbulent flow, and the theory non- 
viscous compressible and incompressible flow. Dryden wrote® 1952: 


“Progress knowledge fully developed turbulent flow processes 
has been more moderate rate, and our present methods solu- 
tion) are still semi-empirical character.” 


The most frequent formulations used engineers continue some form 
the mixing length equations. 


MATHEMATICAL THEORY DIFFUSION 


Semi-empirical formulations the mixing-length type will used this 
study. Recognizing that the diffusion rate the size the eddies 
possible, average size eddy andthe corresponding empirical obtained. 
should noted that the ocean reservoir there are many sizes 
eddies operating within one another increase the transfer area that 
constant concentration attained. Because the random nature the 
eddies, the concentration distribution the ocean reservoirs cannot 
follow exact mathematical law, but must treated the statistical 
sense. That is, enough tests are conducted the averages will tend 
stabilize. Therefore, single tests must approached warily. 

The diffusion equation used the form 


for the elementary case convection one dimension, which, the 
concentration, denotes the eddy diffusion coefficient, the length, and 
represents time. 

The coefficient the constant proportionality between the rate 
diffusion across plane area and the concentration gradient normal that 
area; analogous thermal conductivity heat flow. Eq. equivalent 
the general equation conduction heat flow. The vast amount litera- 


Mechanics Supports Tasks Civil Dryden, Civil En- 
gineering, 22, 1952, 329, 


* 
4 
te, 
2 
4 
4 


154 May, 1961 


ture heat transfer thus available for the more involved studies 
diffusion. The most important source heat-conduction problems 
has been the text Carslaw and Jeager.9 


PREVIOUS STUDIES EDDY DIFFUSION COEFFICIENTS OCEANS 


Among the early investigators, Sverdrup, Johnson, and 
Fleming have tabulated previous results for eddy diffusion 
Their values for horizontal diffusion are 103 105 per sec and 
for vertical diffusion 10-3 per sec. This would imply that vertical 
diffusion can ignored when compared tohorizontal diffusion. They note that 
the eddy diffusion coefficient function the size the eddy. 

Later, Sverdrup, assuming diffusion function average velocity 
and mixture length (that is, the dimension the eddy) and independent 
the depth, hazarded guess values diffusion coefficient.11 For eddies 
radius 103 ft, 104ft, ft, and 106 ft, the diffusion coef- 


TABLE EDDY DIFFUSION COEFFICIENTS 


Diffusion coefficient, Diffusion coefficient, per sec 


(a) Ignoring the size (b) Taking into account the size 
the the eddies 


The diffusion coefficient equals 


ficients were estimated the order magnitude 2.0, 102, 103, 
and 104 per sec, respectively. Sverdrup suspected that these values 
might too large. 

After this, Burke, series off the coast California 
where photographed dye-diffusion patterns from arrived the 
eddy diffusion values shown Table 

The values the diffusion coefficients obtained this test were used 
estimate the spread the contaminant for “Operation Crossroad” (the 1946 
Atomic Energy Commission tests the Pacific Ocean). 


“Conduction Heat H.S. Carslow and Clarendon Press, 
Oxford, England, 1947, 

Hall, New York, 1942, 484, 

“Speculations Horizontal Diffusion the Ocean,” Sverdrup, Scripps 

«Horizontal Diffusion Dyes inthe Ocean,” C.J. Burke, Scripps Inst. Oceano- 


0.116 0.045 0.090 
0.052 0.322 0.645 
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Munk, Ewing, and Revelle measured the concentration 
reef water various parts the Bikini lagoon and arrived the following 
The vertical diffusion coefficient due wind shear 0.2691 
per sec; the horizontal diffusion coefficient for eddy size 103 0.5 
per sec; and for eddy sizeof 104 per sec. They 
present primary constant the eddy diffusion coefficient divided the 
radius the eddy which they found equal 0.0065 0.0164 per sec with 
0.013 per sec being the average. These values were obtained assuming one- 
dimensional flow the analysis the test results. 

Measuring the values for the dispersionof radioactivity following the 
underwater atomic blast, they arrived mean value for surface 
spreading 1.6 102 per sec for equal 9.83 103 ft. This yields 
value a2/r 0.0164 per sec. The coefficient diffusion depth 
150 was about one-third the coefficient diffusion the surface. 

Stommel, photographing sheets paper scattered from airplane 
onto the surface the sea Woods Hole, Mass., series values 
ranging from 2.252 10-3 102 per sec. fairly straight 
line was obtained when the diffusion coefficient was plotted against the 4/3 
power distance separating the sheets paper. The equation used 


F(L) 


which denotes the distance sheets particle, and con- 
stant. This the Richardson Lawof “neighbor diffusivity” and F(L) 3.03 a2, 


FIELD TESTS FOR DIFFUSION COEFFICIENTS RESERVOIRS 


series tests were conducted the diffusion rates reser- 
voirs. Radioactive tracers, because their lowcost, ease detection, small 
gravimetric quantities required, and absence density problems, were used 
measure the eddy diffusion Various reaches Reservoir 
No. the Metropolitan District Commission Framingham, Mass. (Fig. 
were used determine the eddy diffusion coefficient quiescent bodies 
water. This reservoir completely isolated from the potable water supply, 
and fishing and boating are prohibited. 

Method.—The tests were run the following manner: (a) small cove 
basin section the reservoir was selected; (b) transit stations allowing 
full view the reservoir were established; (c) dye and radioisotope mixture 
were affixed the end fishpole and deposited predesignated spot 
the reservoir; (d) the dosing apparatus was removed, and test tube was 
attached the end the fishing pole; (e) the rowboat then circled the initial 
spot and water samples were taken aseries concentric circles about the 
initial dosing point. Because the action the wind and water currents, the 
rate dispersion the tracer seldom exhibited radial symmetry about the 
point dosing. The dye provided useful guide for the selection sampling 


Transactions, Amer, Geophysical Union, Vol. 30, No. 1949, 

“Horizontal Diffusion Due Oceanic Turbulence,” Stommel, Journal 
Marine Research, Sears Foundation, December 30, 1949, pp. 205-214, 

“Radioactive Tracers Hydrologic Studies,” Parker, Transactions, 
Amer, Geophysical Union, Vol. 39, June 1958, 434-439, 
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points following the period initial sampling. the transit stations angular 
measurements were taken locate each sampling point. intersection 
the angular measurements from two stations the exact points sampling 
were plotted the map. This system allowed greater flexibility than could 
have been attained using fixed set predetermined sampling points. With 
the 12-ft fishing pole was possible approach the dye the boat without 
inducing too much turbulence the area. The dye,sodium fluorescein, was 
chosen only after series tests showed have the maximum visibility 
low concentration all dyes tested (eosin, picric acid, Bordeau red, scarlet 


INLET DAM 
RESERVOIR NO. 


DEAD INLET 


UPPER INLET POND 


ROUTE 135 SCALE MILES 


FIG, RESERVOIR NO, 


conc., acid black conc., malachite green oxalate, ecio glavcine, fuchin 
crystal). also the dye chosen for use air-sea rescue kits. typical 
reservoir test shown Fig. and Table 

Analysis.—In the analysis the test, the following mathematical model 
was used: “slug” pure radioisotope added point lake (of 
infinite expanse and constant depth) and rapidly mixed from top bottom 
form instantaneous line source constant strength. The mixing 
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the water assumed isotropic. Thereafter, the isotope subject 
two actions: (a) radial diffusion away from the source and (b) nuclear disin- 


TABLE 2,—SAMPLING POINTS AUGUST 14, 1951 TEST 


Time After August 14, 1951 Warm, sunny, wind 


109 Extremity 
613 
1270 Picked bottom 
Picked bottom 
1844 Picked wind shift from westerly 
1159 Surface wind 
5159 Surface wind 
34110 Surface wind 
30690 Surface wind 
20120 Front path 
16360 Front path 
7118 Front path 
Back path 
296 Back bottom 
Back bottom 
303 Back bottom 
276 Back bottom 
10770 Front path 
11670 Densest path 
dye 
20170 Dense 
13470 Densest 
Integrated depth sample 
Integrated depth sample 
252 Integrated depth sample 
223 Integrated depth sample 
Integrated depth sample 
Integrated depth sample 
Integrated depth sample 
Integrated depth sample 
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Grab samples path 
Grab samples path 
End path 

path 

path 

path 

Dosing point 

End path bottom 
End path bottom 
End path bottom 
End path bottom 


tegration together with removal self-purification processes (precipitation 
and deposition, adsorption, biological uptake, etc.). 
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Consider (Fig. element volumerdr the plane diffusion, 
which the radius from the source, the angular width the element 
and the height. The rate which isotope enters the element diffusion 


The rate which leaves diffusion and die-away 


which denotes the eddy diffusion coefficient reflecting the degree 
mixing the effective decay constant the radioisotope. 
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Equating the difference between these two rates the rate accumulation 
isotope within the element, one obtains 
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The appropriate solution Eq. for the boundary conditions stipulated 


given 


4at 


which the concentration any radius and any time after the isotope 
added and represents the amount radioactive material initially added 
per foot depth. ideaofthe functional relation existing between andt 
given (Fig. 4). The concentration distribution takes the form sym- 


= SMALL— 


FIG. FROM INSTANTANEOUS LINE SOURCE 


metrical bell-shaped hill (Gaussian frequency distribution) with circular con- 
centration contours. With increasing time the hill spreads out and the concen- 
tration the center decreases; also time goes the volume the hill 
decreases because radioactive disintegration and self-purification. 

instructive consider the history one particular concentration 
contour. Let the radius the circle having concentration which 
might correspond critical toxic level tolerance concentration. From 
Eq. found that 
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Eq. indicates that the critical radius first increases rapidly, then passes 
through maximum, and with progressive disintegration zero. The 
recession rapid toward the end. any time the concentrations within the 
critical contour are greater than and outside they are less than 
especial significance the way the initial concentration, appears 
the formula. tenfold hundredfold increase has relatively small 
effect 

Results.—The results this type analysis the pond tests are shown 
Table These values a2/ are somewhat higher than those obtained 
Munk. This may explained, part, the circumstance that Framingham 
Reservoir No. relatively shallow, and uniform concentration distributions 
the vertical direction were quickly established. the ocean, however, 
steep concentration gradients persist for long periods following the addition 
the tracer; the resulting continued the tracer causes 
reduction the dispersion the horizontal directions. From the data obtained 
previously, eddy diffusion coefficients oceans and ponds versus the size 
the eddies are plotted Fig. Ignoring Sverdrup’s estimates, possible 


TABLE 3,—RESULTS FRAMINGHAM RESERVOIR TESTS 


plot line average eddy diffusion coefficients versus the radius the 
eddy. The equation the line 
(6) 
The results this investigation indicate that Richardson’s law neighbor 
diffusivity obtains more than Munk’s suggestion constant for the eddy 
diffusion coefficient divided the radius the eddy. The use power 
function, rather linear relationship, would accord more with the physical 
system where with increasing radius the diffusing surface would increase 
faster than linearly. 


DIFFUSION PIPELINES 


Sir Geoffrey Taylor showed theoretically and experimentally that contam- 
inants pipelines diffused the following manner:16,17 


“Dispersion Soluble Matter Solvent Flowing Slowly Through Tube,” Sir 
Geoffrey Taylor, Proceedings, Royal London, Sect. 219, No, 1137, August 
25, 1953, 186. 

“The Dispersion Matter Turbulent Flow Through Pipe,” Taylor, 
Proceedings, Royal Soc, London, No. 1155, May 20, 1954, 446. 
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Laminar Flow.— 


Dissolved matter spreads through the combined action molecular 
diffusion radial direction and the variation velocity over the cross- 


section. 
The distribution concentration dissolved matter centered about 


mean which moved with the mean speedof the flow and symmetrical about 


it. 
The dispersion along the tube accord with vertical coefficient 


diffusion equal which the radius the pipe; denotes the 


mean velocity the flow; and represents the coefficient molecular dif- 
fusion the injected material the fluid. 

Turbulent Flow.—Taylor predicates that the universal distribution 
velocity pipe obtains (that is, the ratio the difference between the veloc- 
ity the center the pipe and the velocity any point the pipe the 
shear velocity, which the square shear stress the wall divided 
the density the fluid, function the radial distance the point 
question divided the radius the pipe). Further, assumed that Rey- 
nold’s analogy valid (that is, the transfer matter, heat, and momentum 
turbulence are exactly analogous). 


Matter transferred longitudinally convection (the mean velocity) 
with virtual coefficient diffusion equal 10.16 times the radius the 
pipe times the shear velocity. The shear velocity function only the 
Reynolds number for smooth pipes (that is, the coefficient friction). 

Matter transferred longitudinally turbulent (secondary velocities) 
diffusion, with coefficient diffusion equal 0.052 times the radius times 
the shear velocity. (Taylor assumes that the coefficient longitudinal diffusion 
equal the coefficient lateral diffusion.) 

The combined value for longitudinal diffusion coefficient is, therefore, 
equal 10.1 times the radius times the shear velocity. the American sys- 
tem notation: 


which the longitudinal diffusion coefficient, denotes the diameter, 
and the Weisbach-Darcy friction coefficient. 
For open channels 


which the hydraulic radius, represents the gravitational acceleration, 
and, the slope the channel. 

Most the previous theory has held that the transfer due the superim- 
posed turbulent component velocity the Therefore, 
most theoretical formulas for eddy diffusion pipes have been error 
factor about 200. This because they have completely ignored the trans- 
fer which dependent the primary velocity. 


2 
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The concentration dissolved material pipe time and distance 

equal 

which the total mass (or activity) the dissolved material. any 
given time after dosing the tracer, the concentration distribution symmetri- 
cal and fixed station downstream the curve concentration 
versus time not symmetrical due the diffusion the material passes 
the measuring station. Ordinarily, skewness the concentration-time dis- 
tribution not great when the measuring point distance downstream 
greater than 100 diameters. 

Numerical values the virtual coefficient diffusion may readily 
determined taking advantage the fact that the Eq. the error function 
equation. Taylor computed from plot the data scaling off one-half 
the distance which the concentration exceeds the peak value. Solving 
the equation the following relation obtained: 


0.5 


equals 50% the peak value, and the maximum concentration the meas- 
uring station. 

may transposed the form 


which one-half the time that the concentration above 50% the 
maximum concentration, and the distance from the dosing point the 
measuring station. 

Experimental work indicated that the formlation validfor both smooth 
and rough pipes, since the universal distribution velocity holds for both. 

For curved pipe the universal distribution velocity does not hold and, 
hence, the coefficient diffusion should expected greater than that 
computed from Eq. 


OTHER METHODS SOLUTION 


Other treatments Eq. for the evaluation are possible. Eq. may 
written the form 

this form linear relation obtains between log and (x-u t)2/t. 


(12) 
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Another method linearization obtained writing Eq. the form 


which the ordinate; the abscissa; and the 


slope the equation. 
also possible evaluate the standard deviation from moments the 


observed concentration-time distribution. The first two moments the con- 
centration-time distribution, and mg, (Eq. are 


which the theoretical detention period, X/u, and represents the total 


length pipe. 
The second terms the right-hand side Eq. and are ordinaril 


small and may neglected. they are neglected the following equation for 
may obtained solving Eq. and for 


3 


final method the approximate formulation used Archibald, 


ASCE, the standard deviation the concentration-time distribution, 
assumed one-sixth the time between the first appearance and 


last trace the tracer the station. The imperfection Archibald’s 
criterion the determination the time occurrence the last trace, 
since the latter obscured fluctuation background activity. Therefore, 
this method determining the diffusion coefficient much less accurate 
than either the method moments the Taylor method. 


DIFFUSION COEFFICIENTS COMPUTED FROM EARLY SALT-VELOCITY 
TESTS PIPELINES 


The results pipe tests reported the literature have been analyzed, 
and the diffusion coefficients have been determined the methods described. 
The results are presented Table 


“Radioactive Tracers Flow Tests,” Archibald, Journal, Boston Soc, 
Civ. Engrs., Vol. 37, No. 1950, 77. 
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Linearization (Eq. 12) has the disadvantage that forcing the data into 
normal curve, the experimental points that into this pattern produce 
some scatter, and the determination the straight line entails strong ele- 
ment subjectivity. Moreover, the physical significance the data lost 
linearization, the Taylor method (Eq. 11) computing the diffusion coeffi- 
cient the physical sense retained; much simpler and more straight 
forward. 

Taylor’s method requires knowledge the mean time flow (X/u) and, 
therefore, not completely satisfactory. None these methods give satisfac- 
tory résults instantaneous injection not made obstruction occurs 
the line reduce the effective diameter the pipeline not 
straight. Taylor’s method (Eq. and 11) does not indicate this. Therefore, 
the Taylor method should carried out conjunction with the method 
moments. The method moments entails bias the mean time flow 
remains unknown. Then the proper conditions have obtained, the results 
the two methods should approximately equal. 

comparison the values Table 4of the diffusion coefficients obtained 
the various methods makes evident that Taylor’s formula (Eq. 11) holds 
for both oil and water and for pipes flowing full and partly full. 

Ordinarily the coefficient diffusion obtained experimentally may 
expected exceed the theoretical coefficients, since the Taylor formulation 
pertains straight pipe. Existing pipes usually have curved sections. few 
the values Col. Table are less than unity. The discrepancy 
probably due errors the determination the friction coefficient the 
pipe. Also, mentioned previously, the linearization method dependent 
the choice the line best fit. 


TESTS PIPELINES 


series tests using radio-tracers was made determine the virtual 
coefficient longitudinal diffusion pipelines. These results were then com- 
pared with similar coefficients obtained from previous and with the 
coefficients computed the Taylor formula. 

The tests run old water-supply pipeline made available the 
Metropolitan District Commission and the Conservation Department Mass- 
achusetts. The pipe tested was part the conduit system between Whitehall 
Reservoir and Hopkinton Reservoir. These are out-of-service water-supply 
reservoirs and have been assigned the Department Consevation. The 
pipe in. diameter and 6,400 long. The drop the hydraulic grade 
between the inlet and outlet the line ft. The method testing was 
allow the pipe flow full long enough establish steady state conditions and 
then “slug” the pipe with radioactive tracer and dye. the outlet 
the pipe was arranged the outlet was sampled every min. 
until appearance the dye and every 1/2 min. thereafter until the dye dis- 
appeared. The sampling was continued for five six detention periods. The 
concentration the dye served rough check for radioactive measure- 
ments. volume each sample collected was precipitated.The pre- 
cipitate was dried and counted with end window Geiger tube with suitable 
sealing circuits. The results were compared with diluted sample the 
initial dose. 
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The results two the tests are shown Figs and The method 
moving averages was used compute the tails the curves. 

Table much greater percentage iodine passes the downstream 
gauging station than does phosphorus. Only small amount phosphorus 
reaches the sampling station test due the high absorption rates 
phosphorus the bacterial slimes the pipe walls.19 Further credence 
this theory Coffinet al., who state20 “p32 uptake plants 
and sooplankton minutes and not days weeks.” will 
noted Table that the coefficient diffusion Taylor’s method (Eq. 11) 
and the method moments (Eq. 16) not agree too well. 

The apparent reason for the failure the method moments give re- 
sults that agree with other methods is: 


The method moments (Eq. 16) uses “tails” concentration-time data 
which are obscured variations background counting rate. 

Taylor’s theory not good for those portions the distribution-time 
curve large distances from mean. This evident from the shape the 
curves, 

The Taylor method (Eq. 11) not good for skewed curves. The curves 
for the pipe test were skewed because the air pockets which remained 
the pipe. the last test, however, the pipe was vented before the test began, 
and can assume that most ofthe air was removed. The method moments 
then gives approximately the same results other methods. 

Another series tests were run through sewer line New Bedford, 
Mass., co-operation with Frank Heaney, ASCE, Fay, Spofford and 
Thorndike, consulting engineers. The sewer tested was the main collecting 
sewer New Bedford, Mass. Its route Rodney French Boulevard and 
through Fort Rodman, and discharges into Buzzards Bay (Fig. 8). The sewer 
the usual horseshoe cross section combined sewers with maximum 
width in. and height in. The slope approximately 0.001. one 
test the activity was introduced Hudson Street. Samples were collected with 
diaphragm pumps downstream stations the drill ground (1,987.7 ft) and 
the tide chamber (3,429 ft). these pump outlets the stream was sampled 
every min until the dye 20-sec intervals until the 
dye disappeared which time sampling was resumed; typical curves 

another test the tracer was dosedatthe screen house, 1,060 upstream 
from Hudson Street. The same sampling procedure was used with additional 
sampling station Hudson Street. Due faulty injection the activity the 
test bled into the system from the screenhouse well slowly. This reflected 
the high values for these conditions should expected and 
found that the apparent coefficient diffusion much larger than the theo- 
retical value. This emphasizes the necessity for instantaneous injection into 
the turbulent flow. Injection into the laminar sublayer will obviously produce 
laminar diffusion and later turbulent diffusion bleeds into the turbulent 
flow. Attempts measure the eddy diffusion coefficient the outfall were 


and Assimilation P32 Bacterial Slimes,” Reid, San, 
Dept., John Hopkins Univ., Contract AT-(10-1)-498, March 20, 1950, 
“Exchange Material Lake Studies the Addition Radioactive Phos- 
phorus,” C.C, Coffin, Hayes, L.H. Jodrey, and Whiteway, Canadian Jour- 
nal Research 27, 1949, 221, 
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not successful due the small amount radioactive isotopes used and the 
high solids sea water. 

The tests Massachusetts bear out the conclusions drawn from previous 
tests. With proper injection methods and measuring procedures the Taylor 
formulation and the method moments give satisfactory results for circular 
and horseshoe pipes, full and partially flowing full pipes, and oil, water, and 
sewage carrying pipes. 


relatively easy determine diffusion coefficients reservoirs and 
pipelines use comparison with most other methods 
heretofore used. Simple dosing and sampling devices can used obtain 
good results. 


ACKNOWLEDGMENTS 


The work described herein was performed under contract with the Sanitary 
Engineering Section the Division Reactor Devlopment the AEC 
Harvard University, under the direction Harold Thomas, J., ASCE. 
The writer would like acknowledge with thanks the help Thomas guiding 
the project and the participation Stevens Kleinschmidt, Mrs. Klein- 
schmidt, Jane Maclachlan, Thomas Ostrom, ASCE, and Gerald Parker, 
the execution the field tests. 

This paper has been submitted part the requirements for the Ph. 
Degree Harvard University, Cambridge, Mass. 


4 
4 
4 
i 
i 


a 
: 
. 
. 
: 
‘ 


2827 May, 1961 


Journal the 


HYDRAULICS DIVISION 


Proceedings the American Society Civil Engineers 


DISCUSSION 


paper part the copyrighted Journal the Hydraulics Division, 
Proceedings the American Society Civil Engineers, 87, May, 


173 


a8 
th 
a 
ory 
7 
3 
of 
| 
( 
3 
K 


id 
4 
Pan i 


EARLY HISTORY HYDROMETRY THE UNITED 


Closure Steponas Kolupaila 


STEPONAS writer feels particularly honored the 
fact that three the pioneers presented the paper took part the discus- 
sion: Stevens, ASCE, Hoyt, ASCE, and Blanchard, ASCE. 
not enough thank all contributors for their valuable cooperation; let 
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FIG, 63.—HAMILTON SMITH FIG. 64.—LOAMMI BALDWIN 


suffice say that the writer inspired, the attention his colleagues, 
continue further study the history hydrometry. 

Johnson deserves credit for recalling the “Early Engineering Center 
California,” good addition the chapter “Water the West.” Hamilton 
Smith world renowned for his authoritative that now very rare. 
The classical experiments carriedout Smith are mentioned modern text- 
books. The portrait this pioneer hydrometry (Fig. 63) reproduced 


January 1960, Steponas Kolupaila (Proc. Paper 2335). 

Engrg., Univ. Notre Dame, Notre Dame, Ind. 

“Hydraulics, the Flow Water Through Orifices, Over Weirs, and Through Open 
Conduits and Pipes,” London, 1886, 
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here87 the courtesy Daugelaand Bedford. Smithdied while sail- 
ing down Oyster River. 

Johnson mentioned his the first the world’s long-distance 
high-voltage power transmission lines was built California. connection 
with this statement, Ellis, ASCE, insisted letter the editor that 
the Missouri River Power Company had built 70,000-v line extending from 
their Canyon Ferry plant the town Buthe, distance miles, the 
year 1901, thus claiming predate the California project. The writer would 
like point out, however, that the 60,000-v line running 140 miles from the 
Colgate plant the city Oakland, California, that was described Johnson, 
was built the year 1899, thus refuting Ellis’ 

Hoyt, one the leading pioneers the USGS, from 1907 1950, recom- 
mended the inclusion several hydrographers that institution. the 
feeling the writer that these persons could more appropriately included 
the modern history hydrometry. Hoyt emphasized the manifold indirect 
duties hydrometric engineers, disseminators public information and 
educators the necessity conserving our water resources. expressed 
his opinion, aconsultant Senate Committee Water Resources, con- 
cerning the future those problems intimately affiliated with the history 
hydrometry the United States. The annual expenditure the USGS for in- 
vestigationof water resources, about million dollars, demonstrates the high 
recognition this important work. 

Shaw Cole, ASCE, who skillfully continues the illustrious work his 
father, Cole, mentioned some new developments, that, the writer feels, 
belong more appropriately the modern history hydrometry. 

Blanchard, one the distinguished pioneers the Lake Survey, added 
more infcrmationon the measurements Sault Ste Marie, and work 
different rivers connecting the Great Lakes. 

Frazier, who inspired this study, contributed valuable improvements 
with which the writer agrees completely. The pioneer work Loammi Bald- 
win, (Fig. 64) who has been referred the father civil engineering 
America, should included the early history hydrometry. 

The writer grateful Frazier for disentangling the amazing coincidence 
the two Nicollets, both investigators one part our country, but 200 
apart. 

The first was Jean Nicolet, Canadian explorer, born Cherbourg, France, 
about 1598, who immigrated Canada 1618, travelled across the Indian ter- 
ritory, that time called New France, and was the first white man visit 
Wisconsin. didnot reach the Mississippi River, however, and later returned 
Quebec. drowned the year 1642. river, country, and town 
Canada bear his name. 

The other was Joseph-Nicolas Nicollett, 1786 1843, French mathematician 
and cartographer, who came tothe United States 1832 the upper 
Mississippi the years 1836, 1838, and 1839. His report, written French, 
was translated English Ducatel and Alexander, and published 
Congressional 

The improvements the material concerning Hall, State Engineer 
California, and the portrait furnished Frazier, are highly appreciated. While 


the Town Durham, New Hampshire,” 1913. 

Report intended map hydrographical basin upper Mississippi 
River made Nicollet. Senate Document 237, 26th Cong., 2nd Session, 1843, 
Serial No. House Document 52, Vol. 2nd Session, 1845, Serial No, 464, 
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true that several methods were initiated Hall, who published large 
volume hydrometric data collected California, and several annual re- 
ports irrigation from 1880 1888, his methods were mentioned papers 
his assistant Grunsky only, and therefore were attributed Grunsky, 
although Grunsky himself did not claim them personally. This information 
Hall and Grunsky was corrected the recently published bibliog- 
raphy. 

Stevens, the eminent pioneer hydrometry, was particularly favor- 
able the writer’s intentions. recalled the origin the widely publicized 
theory our forests, that was deliberately falsified 
order further conserve the forests. This was done connection 
comparison droughts the Ohio basin 1824 and from 1925 1930. 
was particularly interesting know the personal opinion one the authors 
that story. Careful hydrometric investigations carried out experimental 
forest basins different countries had thoroughly explained the true hydro- 
logic role forests and clearings. 

Moskvitinoff presented some conclusions from his hydrometric work 
Central Asia from 1914 1915. The extensive investigations the accu- 
racy flow measurements, the best investigations known the writer, were 
included the 5th volume the Report the Hydrometric Department for 
that, regrettably, was never published. Fragments this work ap- 
peared Bulletins the Hydrometric Department Turkestan 

Among the many letters received the author connection with his pa- 
per was humorous contribution received from George Henry Ellis, ASCE, 
who wrote: 


“Years ago, when was swinging the current meter, there was 
group spectators asking foolishquestions. The foolishanswer some- 
times gave them grew, over periodof about thirty years, into the story 
the counting the fish.” 


The enclosed ballad entitled “The counting the fish” 1946 contains verses 
like “fishes cannot live water near power plant, with kilowatts squeezed 
excerpt good example hydrometric humor. 


APPENDIX BIOGRAPHICAL DATA 


Loammi Baldwin, May 16, 1780, North Woburn, Mass., June 30, 1838, 
Charlestown, Mass.; 1800 graduated Harvard Univ., built harbors, 
1821 Union canal Penn., 1825 Boston, built Connecticut-Hoosac tunnel; 
called Father civil engineering America. 


«physical data and statistics California,” W.H. Hall, Tables and memoranda 
data, collected and compiled the State Engrg. Dept. Sacramento, 

“Bibliography Hydrometry,” Kolupaila, Notre Dame, Univ. Notre Dame 
Press, 1961, See references 17, 157, andH 12, 

“Bibliography Hydrometry,” Kolupaila, Notre Dame Univ., Notre Dame 
Ind,, Press, 1961, See reference 

“Bibliography Hydrometry, Kolupaila, Notre Dame Univ. Press, Notre 
Dame, 1961, See reference 294, 
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William Hammond Hall, 12, 1846, Hagerstown, Md., Oct. 16, 
1934, (only this information was available). 1879 State Engr. California. 


David Carlisle Humphreys, Oct. 14, 1855, Chatham Hill, Va., Jan. 
10, 1921, Lexington, Va.; 1878 graduated Washington Lee Univ., 1879 Mo. 
River investigations, 1885 1921 Prof., 1904 Dean, Washington Lee Univ.; 
1895 1908 Geological Survey; 1897, James River expedition. 


Elwood Mead, Jan. 16, 1858, near Patriot, Ind., Jan. 26, 1936, Wash- 
ington, D.C.; 1882 graduated Purdue Univ., 1886 Prof. Irrig. Engrg. Colo- 
rado Agric. College; 1899 Chf., Irrig. Investigations, Dept. Agric.; 1907 
1913 river studies Victoria, Australia, 1913 Prof., Univ. California, 1917 
Calif. State Bd. Irrig., 1923, 1927, Cons., Palestine, 1924 1936 Commis- 


FIG, 69.—ELWOOD MEAD FIG, HANNA 


sioner, Bur. Reclamation, promoter Hoover Dam; reservoir Colorado 
was named Lake Mead his honor. 


Edwin Nettleton, Oct. 22, 1831, Medina, Ohio, April 22, 1901, 
Denver, Colo.; Medina Academy and Oberlin College, Ohio, 1870 Colo- 
rado town planner, 1883 1886 State Engr. Colorado. Adopted initial, 
middle name. 


Joseph-Nicolas Nicollet, July 24, 1786, Cluses, Haute-Savoy, France, 
Sept. 11, 1843, Washington, C.; cartographer, 1817 Paris astronom- 
ical observatory librarian, Prof. Louis College; 1830 ruined stock market, 
1832 immigrated S., 1836 1839 Miss. River investigations. 
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FIG, TRAINING CAMP EMBUDO, M., 1888 
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FIG, 71.—EARLY WEIR MEASUREMENT, LEFFEL, 1881 
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Hamilton Smith, July 1840, near Louisville, Ky., July 1900, near 
Durham, H.; 1854 work mining and hydraulic engrg, 1869 1881 Calif- 
ornia, water supply for mining, 1881 Venezuela, London, 1892, South Africa, 
Alaska, mining, 1895 New York, cons. engr. 


Extending cordial thanks all contributors, the author has taken the oppor- 
tunity toadd several portraits located through the kind cooperationof Frazier. 
They are: Sexton, the first manufacturer current meter the United 
States; Nettleton, State engineer Colorado; Rafter, State engi- 
neer New York; Humphreys, investigator the James River (cour- 
tesy the Washingtonand Lee University Library); Mead and Hanna, 
pioneers the Bureau Reclamation (courtesy the Seattle Public 
Library). The interesting sketchof the early weir measurement93 also pre- 
sented. (courtesy Purdue University, Gorr Library). The historical group 
the members the first hydrometric training camp, near Embudo, New 
Mexico, 1888, must also included the early history hydrometry. 


“Construction Mill Dams,” Leffel and Co., 1881, 
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GENERALIZED DISTRIBUTION NETWORK HEAD LOSS CHARACTERISTICS* 


guments, clarifications, and original supplementary information presented 
the thirteen engineers who contributed the nine discussions have placed the 
subject much clearer perspective. 

Little would gained reviewing the discussions point point cit- 
ing the multitude statements with which the writer full agreement. 

Rather, the questions raised will answered and clarifications offered. 

Specific Details.—The writer stands corrected Robertsonon the improper 
use the term “relaxation.” 

Constant mentions that the flows each pipe can checked during net- 
work analysis means Eq. For all the network runs used the pa- 
per, each junction and loop was checked insure that the algebraic sums 
all the corresponding flows and head losses were zero, with due allowance for 
attainable precision meter readings. For the initial run each network 
series, each and every pipe was checked via the equivalent Eq. There- 
after, every second third run all the pipes were completely checked 
again insure that all “fluistors” were operating properly. This stand- 
ing operating procedure the Philadelphia Water Department; leaves little 
chance for error. 

The error analysis Constant valuable supplement. Some indication 
relative error should have been included the paper. 

Constant contends that the procedures presented are principally applicable 
field operations. Radziul and Celenza have clearly demonstrated the value 
these procedures for operational purposes. However, equally not more 
important use the system analysis phase design, will explained 
subsequently. 

Bitoun has referred the number conditions required define Fig. 
and for determination and Eq. The “existing” district Fig. was 
studied secure more exact understanding its characteristics because 
extensive, synchronized field tests were not practicable due its large size. 
This preliminary study was great value subsequent studies modifica- 
tions required conditions. The field data were included Fig. for 
verification. More than one run was made for the Fig. presentation prove 
unequivocally that Eq. was valid under the prescribed conditions. was felt 
that single run would not adequate for that purpose with such complex 

network, even though only one run required define must em- 
phasized that major consideration the designof public utilities the pro- 
vision adequate facilities for future needs. Seldom will existing system 


January 1960, McPherson (Proc. Paper 2339). 
Prof, Engrg., Civ. Engrg. Dept., Univ. Urbana, 
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capable handling anticipated higher future demands. Also, modern com- 
putational devices make possible more exact and realistic studies. Indeed, 
the inadequacy many existing facilities can attributed the comparatively 
primitive computational procedures available and vogue during the era 
which they were designed. 

Radziul and Celenza have clearly demonstrated the efficacy using Eqs. 
and determine the characteristics existing prototype systems. The 
modifications planned for these two Figs. and 6(a)], either pro- 
gress scheduled, will vastly change their existing (Tables and charac- 
teristics. Because variations piping and source and storage facilities can- 
not simulated under prescribed loading conditions with prototype system, 
they must simulated means some type model. For this simulation 
one network analysis required define Eq. and two analyses are 
required define and Eq. Obviously, matters not what computing 
device used analyzing particular loading run, paraphrase one Lo- 
max’s observations. 

Lomax states that was “not stated, but implied that the source 
magnitudes must proportional well the load magnitudes.” This would 
true only the specific instance network with single source and 
equalizing storage: Eq. would apply for all system demands; each local de- 
mand load would fixed fraction here, because equals the 
several loads would also equal fixed fractionof Qp- network fed two 
sources hydraulically quite similar toa single-storage, single-source com- 
bination (see Figs. and when the equalizing storage contributing the 
demand less than Two more true sources, whether they pumps 
gravity storage, can generalized use Eq. The writer has made 
effective use this procedure inthe system analysis for asingle-district large 
city that has two gravity storage reservoirs dissimilar levels, pumping 
station remotecorner with second station ina different quadrant needed 
for the future. 

The system analysis the single-source case without equalizing storage 
obviously routine matter. With multiple sources, with single multiple 
sources combined with equalizing storage, the number hydraulic parameters 
that must satisfied increase. The capacity water levels 
storage, together with pump characteristic curves, must reconciled with 
the network head loss characteristics under variety known anticipated 
service conditions. The demonstration storage-network-pump balancing 
that will presented relatively simple example final equilibrium 
phase system analysis. The paper was addressed tothe issue gener- 
alization network characteristics. preparing the paper was feared that 
detailed treatment system analysis would detract from and overextend the 
subject, but evident from some the discussions that the paper should 
have been extended somewhat for more firm appreciation the potential 
uses the concepts 

Applicability Proportional Load Principle Non-Hydraulic 
Trunk and Maier and Miller have clearly indicated that the proportional load 
principle little value the analysis gas system networks. Trunk re- 
ports (correspondence dated March 23, 1960) that prior acquisition 
Analyzer 1953, his staff “resorted the trick thinking net- 
work asa single equivalent main,” when complete analysis was not warranted 
via manual computation. The relation between the equivalent pipe method and 
Eq. will discussed later section. 


ars 
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Mine ventilation networks have means special Mc- 
Network Analyzer the Bureau Mines, Health and Safety Re- 
search and Testing Center, Pittsburgh, The resistance 
mine airways proportional the square the volume rate flow. Al- 
though the basic procedure for analysis mine networks many respects 
similar that for other fluid media, loading conditions are rather unique. 
Harris, Acting Chief, Ventilation Group, Bureau Mines, Pittsburgh, 
describes the situation thusly (correspondence dated May 19, 
1960): 


“We use form Eq. however, the application ‘proportional 
loading’ mine ventilation would very limited. Where minor fan ad- 
justments (speed blade settings) are made with other changes 
the system form Eq. and proportional loading would and ap- 
plied. Where adjustments distributionare entailed, use Eq. would 
limited major airways which the resistance factor remains con- 
stant. Other factors would enter into the application.” 


all network analyses, accurate results can achieved only through use 
accurate pressure and quantity field survey data. 

Alan Byers, Head, Electronics Branch, The Franklin Institute, Philadel- 
phia, Pennsylvania (correspondence dated May 26, 1960) has carefully studied 
the paper: 


“My conclusion, after such study, was that the paper does not present 
analysis useful for load flow studies power system practice. 
order avoid having this negative opinion unchecked, engineers the 
Philadelphia Electric Company were asked study the paper also. Their 
electrical engineers system planning find that there nothing useful 
them the paper. They also had engineers fromtheir steam and gas 
distribution systems, which might present closer analogies, 
read the paper. Here the paper was found general interest, but not 
useful any specific way. seemsthat eachof these fields engineers 
have developed analytical methods welladapted totheir problems. Doubt- 
less, drawing the correct analogies, the methods developed the pa- 
per might also used these other fields. But there seems 
advantage gained.” 


Reversing the these considered comments, apparent that 
the engineer engaged water distribution network analyses should exercise 
caution intransferring adapting his use have 
been found satisfactory the electric power systems, gas, steam, 
mine ventilation networks. The word “water” deliberately was not incor- 
porated the title the paper because was mistakenly assumed that Eq. 
might prove useful the analysis non-hydraulic systems. 

has described linear resistance analogue, and the 
non-linear type, for gas network analysis. 

Fire Flow Analysis via Proportional Load Hydraulics.—The statement: 
“However, fire load cannot regarded representative proportional 


Network Analyzer for Solving Problems,” 
McElroy, Bur. Mines Information Circular 7704, November, 

“Calculators Solve Flow Problems,” Clennon, Gas Octo- 
ber, 1951, pp. 68-70. 

“Direct Reading Electric Analogue Computer Helps Solve Distribution Problems,” 
Gas Age, July 17, 1952, pp. 31-35, 
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load,” was apparently misconstrued the writers two the discussions. 
was not intended imply that fire flows cannot analyzed via the propor- 
tional loading principle. Linaweaver, Geyer, and Wolff have correctly explained 
that fire flow load per non-proportional, just the load imposed 
elevated storage while being filled. elementary network shown Fig. 
illustrate this point. The fire flow termed and the head loss from 
the fire flow point termed The data for the four runs plotted 
Fig. are given Table apparent that Eq. appropriate for the 
analysis, long each non-fire load fixed proportion the total non- 
fire load demand. 


TABLE SIMULATED NETWORK FIRE FLOW RUNS (SEE FIG, 10) 


Analyzer Runs 


Using the standard fire flow extrapolation formulato anticipate run F-4 from 
runs F-1 and F-2: 


0.54 


and run F-4 from runs F-1 and F-3: 


versus the true F-4 value 8.0 mgd. 

The elementary network given will usedas second example. 
Table are given the appropriate details for solution via Eq. Note that 
the true higher fire flow 1,000 whereas for the standard formula 
740 gpm. 

These examples tend exaggerate, some extent, discrepancies between 
the standard formula and true values. Only the most remote point the net- 
work was considered. may noted also that the larger the drop the ref- 
erence flow (39.7 ft, 71.9 and 2.3 ft) the better the formula approximation. 
understood that the National Board Fire Underwriters (NBFU) nor- 
mally requires psi drop before sample considered adequate 


Discussion “Friction Losses Lines with Service Connections,” Tong, 
Proceedings, ASCE, Vol. 86, No. November, 1960, 160. 


med 7.96 12.07 15.90 
mgd 8.07 8.09 7.90 
mgd 4.00 6.00 8.00 
46.4 47.6 47.8 45.8 
feet 17.9 57.6 89.8 129.1 
0.386 1.21 1.88 
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for rating purposes. Anderson reminds that for adequate design, fire flows 
must available under maximum day demand conditions. Whereas the stand- 
ard formula yields extrapolated flows that are conservative (that is, low), the 
NBFU seldom extends field test data over large increment beyond the test 
flows estimate capacities psi residual. 

should noted that the standard formula paradoxical, inasmuch 
strictly applicable only when the fire flowapproximates proportional sys- 


PIPE COEFFICIENTS 


NETWORK FIRE FLOW RUNS (SEE TABLE 


tem load. This approximationis fairly exact the extreme: when the fire flow 
insignificant part the total system demandor when practically the 
entire system demand. The first instance local phenomenon because the 
effect the additional system head loss attributable transport the fire 
flow through the network relatively small. the second instance the fire 


+0.4 
+0.2 
-0.2 
127 
1.80 
1.75 
1.70 
3.60 
LOG 
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flow occasions the entire system loss and the standard formula then equiv- 

the preceding two examples there equalizing storage. With equal- 
izing storage there would flow out the reservoir(s) and the fire flow 
might regarded additional, filling reservoir without fixed con- 
trolled water level. The use Eq. for fire flows with without equalizing 
storage conjunction with machine network analyses would pointless be- 
cause the design fire flow alone required and normally can determined 
directly. However, evident thata field test agiven site would require 
extrapolation the required flow, standard formula extrapolated flow based 
intermediate fire demand would more realistic terms ultimate 


NBFU 
has presented clear description the procedure for group 


hydrant flow test. 
Unique Solution Feature Networks.—Bitoun has referred demonstra- 


tion uniqueness given The demonstration inferential 
TABLE 9,—FIRE FLOW DATA FOR SMALL NETWORK GIVEN TONG 


feet 


14,0 


proof, limited consideration boundary conditions too generally 
treated. summarized: “The unicity the solution depends uniquely 
the fact that conduit characteristics are increasing functions.” This say 
that addition being continuous function the head loss characteristic 
pipes bilateral, inasmuch the magnitude the loss independent the 
direction flow. (The continuous, bilateral nature resistive electrical im- 
makes possible the useof A.C. electrical analogues inthe analysis 
fluid flow networks.) rigorous proof for appears possible through 
use special formulations matrix Attainment rigorous proof 
uniqueness for network equations with non-integer unlikely; 
inductive proof uniqueness follows. 

Consider any balanced distribution network for which the flow each main 
branch (both ends terminating intersection junctions) may described by: 


Because the network questionis balanced, the boundary conditions have been 


“Basic Water Works Manual,” Sweitzer, Amer. Concrete Pipe Assoc., 


Chicago, 1958. 
“Circuits Electrical Engineering,” Vail, Prentice-Hall, Inc., New 


York, 


Amer., 1950. 
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zero, 


i=1 


which the number separate flows junction (necessarily three 
more number) and the total number junctions; the algebraic sum 
head losses around each loop zero, 


i=b 
i=1 ---L 


which the number main branches loop and the total number 
loops. 

total number main branches plus one. Immediately the loads are as- 
signed values the branch flows remain the only unknowns. Because the input 
equals the sum the loads the input junction equation superfluous. The 
number equations L+J=1 equals the number unknowns and, therefore, 
solution always possible. 

For the balanced network under consideration Eqs. 13, 14, and have been 
exactly satisfied. Rewriting Eq. 


reduces its original bilateral non-linear form bilateral linear equivalent. 
The sign for each Eq. must determined from the sign the corre- 
sponding balancing branch Eq. (note that when the sign 
intrinsically positive Eq. 13). Because the magnitude and 
fixed and constant, each has singular constant value without sign, 
equal the magnitude only Qm-1), follows that the transformation 
into linear equivalents can have bearing the integrity Eqs. and 15. 
the equivalent linear form Eq. together with Eqs. and con- 
stitute system with unique solution; the solution for the original system 
equations also must unique. 

Determinants will used investigate the simultaneous solution the 
junction and loop equations for the balanced network with equivalent linear 
branches. For the jth root the solution 


which flow, defines the determinant containing the equation 
constants and isthe determinant the equation coefficients. the con- 
stant for each loop equation would zero (except when booster pump 
interposed ina loop). However, loads and inputs are non-zero constants 
the junction equations for active network andconsequently always un- 
equal zero. Because unequal zero, the system equations non- 
homogeneous and required that unequal zero for unique solu- 
tion.33 will equal zero this application only when any two its rows 


“The Mathematical Solution Engineering Problems,” Salvadori and 
Miller, McGraw-Hill Book New York, 1948, 
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columns have proportional elements. The order 
requiring representation for each and every branch flow each its rows. 
For network with two more loops distinct and independent set zero 
values characterizes each equation, precluding the possibility proportional 
determinant elements, insuring non-zero value for (and and conse- 
quently unique solution. The unique solution will yield single, real value 
for the root characterising each branch flow. From preceding arguments the 


TABLE EQUIVALENT LINEAR NETWORK DETERMINANTS 
(Balanced network Case Fig. 


For the system seven equations seven unknowns, where 16), the 
equations terms each branch are: 


Con- 


For branch 


+0.0180 


From which Therefore only one value can exist for gpm). 


corresponding solution with Eqs. 13, 14, and also must unique. Table 
has been prepared the non-zero, nonhomogeneous character the 
determinants, using the Case balanced network data Fig. 

Only single source has been discussed, but the proof applies well 
multiple sources with without storage. For the special case the 
preceding proof tacitly rigorous. 
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appears that Camp and were the first adapt linear resistance 

electric computing boards the analysis water distribution networks, 
means expression with the same meaning Eq. 16. and error 
procedure was used find the set values that would satisfy the corre- 
sponding values, using appropriate conversion factors reduce electrical 
hydraulic units. Much later, convergence procedures were 
that vastly reduced the number trials. Strangely, solution unique- 
ness was Offered. Analysis electrical analogy satisfaction 
Kirchhoff’s laws (1845); with A.C. circuits, vector rather than algebraic sums 
must used. Fich and give valuable discussion the elements 
network topology and also cite special types electronic network compo- 
nents that are associated with homogeneous systems equations solution 
indeterminate form may exist). Van Valkenburg38 gives brief attention 
resistive network analysis, but his book and number others reviewed 
special effort has been made prove general solution uniqueness, perhaps 
because the comparatively small importance purely resistive cases and 
emphasis general circuit characteristics. 

may noted passing that electrical resistive analogues, using either 
linear non-linear resistance elements, can manipulated satisfy dif- 
ferent values for any all pipes the network. Reasons why the use 
mixed values normally unrealistic have been mentioned previously 
(3). Robertson has presented very convincing, comprehensive arguments for 
the use single value for analyses via digital computer. Lomax has 
called attention tothe fact that laminarflow may occur some pipes net- 
work, for which but this consideration should not critical other 
than very smallarterial networks for unusually low minimum demand rates. 

Bitoun’s demonstration the validity Eq. valuable contribution. 
Several decades ago one the standard design procedures consisted re- 
ducing network the equivalent pipe method single pipe. order 
this all loads had consolidated into single load. The method was 
particularly useful for small network fire flow analyses. While this equivalent 
pipe approach repeatedly mentioned the literature, appears that Do- 
land19 was the first suggest the simplified version expressing the flows 
each pipe percentage the input. This principle has been extended 
without proof multiple loads inelementary network examples several text 
books. While writing the paper, the writer considered developing Eq. via 
equivalent pipes but decided against it. order fact that Eq. 
restricted single value and proportional loads, the simpler case 
without equalizing storage was offered three levels complexity in- 
troduction. Because all discussions relating the equations indicate 
clear understanding these limitations, the approach used appears have 
been successful. The demonstration the validity Bitoun could 


Analysis Water Distribution Systems Means Electric Net- 
work Camp and Hazen, Journal, N.E.W.W.A., Vol. 48, Decem- 
ber, 1934, pp. 

“Network-Flow Analysis Speeded Modified Electrical Analogy,” Perry, 
Jr., Vierling, and Kohler, E.N.R., September 22, 1949, pp. 

“Use Alternating-Current Network Calculator Distribution System Design,” 
Suryaprakasam, Reid, and Geyer, Journal, A.W.W.A., Vol. 42, 
December, 1950, 1154-1164, 

A-C Circuits,” Fich and Potter, Prentice-Hall, Inc., New 
York, 

“Network VanValkenburg, Prentice-Hall, Inc., New York, 1955, 
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extended tothe storage case However, boundary value approach will 
not lead full the latter. the writer the more 
complex case because his repeated attempts arrive solution direct 
derivation proved unsuccessful, Eq. would not have been developed, albeit 
empirically via similarity hypothesis. 

Returning the preceding proof solution uniqueness, each and con- 
stant the system simultaneous linear equations for multiplied 
any particular non-zero constant, the solutionis unchanged except that the roots 
are equal the original roots times the constant. Bitoun’s demonstration for 
Eq. indispensable for non-integer values 

System Analysis via Proportional Load Hydraulics.—The information 
Table given response the request Linaweaver, Geyer, and Wolff 
for clarification Table flow ratios. The corresponding hydraulic grade 


9.4 

6.6 

6.6 

18.9 +9.7 
17.4 
26.9 


(For minus sign indicates flow out storage indicates flow into stor- 
age.) (Non-proportional load: Chestnut Hill; both series and 1.2 mgd; both series 
and 0.8 mgd.) 


data for this Roxborough H.S. test series are given Table 12. Forthis series, 
realistic water levels andstorage rates were assigned for the Rox. standpipes, 
taken thereafter reference control. The other water levels are from 
the network analyses, using this datum. The pump total dynamic head require- 
ments for uniform daily pumping listed Table would approximated 
quite well the provisional pumps described the bottom the table [see 
also Figs. The W.O.L. standpipe water levels are reasonable, 
realistic and satisfactory for first phase design except for runs and 6B. 
However, these two runs are for outage W.O.L.P.S. the maximum day, 
emergency rare combination circumstances. The sys- 
tem layout Fig. was the first studied, and not the final arrangement adop- 
ted. 

Linaweaver, Geyer, and Wolff have emphasized, the computations pre- 
sented for the Roxborough H.S. District are merely first trial feasibility 


‘ 
No. mgd mgd mgd =QpR mgd mgd mgd 
(1) (2) mgd (3) (10) stand, 


DISCUSSION 
study, and thorough analysis would still required”----“accounting for 
changes levels tanks and standpipes.” The procedure for balancing that 
they describe could not have been applied effectively for this particular study 
because the proposed W.O.L. pumping station and standpipe, new pumps for 
the Rox. pumping station, and new major arterial mains were all unresolved 
design factors the beginning the study. should noted that Table 
could extendedover considerable number other combinations means 
Eq. for balance study investigate other combinations inputs. 
Constant the opinion that the developments offered the paper are 
little value tothe person concerned with improvements for particular distri- 
bution system, apparently because Eq. restricted toa given loading distri- 


TABLE 12,—WATER LEVELS AND PUMP REQUIREMENTS FOR TABLE ROX, 


Rox, Stand- 
pipes (Sp), 
(Con- 
trol) Feet 


Feet 


Provisional Pumps for Table Study; Acceptable for Above, Special, and Future Demand 
Requirements: 


Rated Capacity, Mgd Rated T.D.H., Shutoff, 
No. No, No, Feet Feet 


110 160 
170 260 


bution, network configuration, and set flow coefficients. given system de- 
sign can reduced the hydraulic characteristics four basic components: 
network, pump performance, pumping station and suction source, and storage. 
These must integrated and reconciled with the design loadings. the net- 
work proper the only basic component subject modification the design, 
Constant’s position well This would rather special design case. 
Nevertheless, Eq. could still inestimable value balancing final and 
alternative network schemes. Entirely too many designs have been based 
fragmentary system studies. comprehensive system analysis required 
insure that the benefits claimed for the design recommendations offered are 
physically attainable. From this standpoint difficult see how “the oper- 
ations end distribution problems” can considered independently from 


Run Standpipe T.D.H. 
Feet Feet 
509.0 110 514.0 171 519.6 
512.2 509.2 109 514.0 170 
511.8 511.6 109 163 516.0 
511.8 513.0 190 534.6 
509.0 505.9 122 523.9 
511.8 509.9 115 519.6 
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design alleged operable. Even perfunctory system analysis least 
four normal consumption demand rates should studied. Only the network 
analyses for two demand rates are necessary obtain the constants Eq. 
The network hydraulics for other rates can then computed using the equa- 
tion, precluding the necessity for more than two runs for that particular net- 
work, Use Eq. particularly effective when two more basic compo- 


TABLE STORAGE-NETWORK-PUMP BALANCING USING 
PROPORTIONAL LOAD PRINCIPLE; BELMONT 


via 
Eq. (2), 
Level 


Storage Storage 


Suction 
Water 


Feet 


(3) 


Notes: For million gal storage 25-ft, over interval. 
Storage water level (6:30 a.m.) hours taken 355.00 (Full). 
See Table 
See Fig. 6(a) and Table 


nents the system are subject modification the design, permitting in- 
several alternatives inthe timethat would otherwise expended 
only one combination. 

Proof the operability given design best demonstrated means 
some sort system balance. noted Linaweaver, Geyer, and Wolff, 


Time Pumping (By Cal, Actual 
hours Inter- feet 
nal, 
(1) (2) (4) (5) (6) (7) (8) (9) 
8-9 10,2 9.52 -0.35 7.70 111.8 
10-11 10.5 242.8 9.56 -0.49 7.38 111.6 
13-14 242.5 9.80 -0.52 7.68 109.5 109.7 
16-17 244.5 9.57 +0.61 344,80 11,21 111.5 
17-18 245.1 9.32 345.60 12.79 113.4 
21-22 247.8 8.30 +1.67 351,83 16.65 120.7 120.7 
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detailed hour-by-hour balance not always necessary desired. However, 
the more detailed the balancing the more indispensable Eq. becomes, pro- 
vided proportional loading can assumed. 

The use Eq. the balancing system can best illustrated 
means detailed example. storage-network-pump balance the Belmont 
H.S. District (Fig. 6(a), Table has been presented (5), but without benefit 
Eq. The system balancing outlined Table uses the same demands, 
storage, network and pump previously noted example (5). 
the details Table 13, the pumping station suction loss was taken con- 
stant (for 9.3 mgd, that is, the ave.). this way the clear well water 
level and suction loss could combined terms “pumping station suc- 
tion water level.” Other than this abbreviation, details are the same pre- 
viously given the writer (5), but the balance more exact. 

Each hourly demand Table was balanced trial sequence. Each 
trial value leads uniquely different value approximately 
sixty trial combinations were needed for the preparation Table 13. Without 
Eq. considerable number would have been required for 
comparable exactitude. Recall that only two network analyses are needed 
determine 

realistic system balance possible only when the basic system compo- 
nents are compatible. The Philadelphia Water Department developing 
system balance digital computer program that will incorporate Eq. and sub- 
routines for pump characteristics loss. The pump- 
ing station hydraulics will because and sta- 
tion losses are non-linear independent design factors. The comparative in- 
fluence diverse andstorage capacities inthe balancing 
for given networks will investigated under variety demand schedules. 
The principal objective the study the determination criteria for 
optimum equalizing storage andpumping facilities. was 
study would not feasible the absence the generalized network param- 
eter, Eq. 

Corrections and Additions.—A recent article should added 
the series given the list references (6). 

The last equation the paper should corrected stated Constant. 
The word “now” should substituted for “not” the first line 59, 
1960. Inclusion the following phrase enclosed parentheses would clarify 
the sentence page 127, HY9, 1960: “For systems without storage with 
storage ator nearthe pumping station, proportional load hydraulics (via Eq. 
accurately------ Reference appeared Vol. 33, No. pages 234 236, 
February, 1941. the line preceding Eq. page HY7 1960, “In 
any loop---” should changed read, “At any junction---.” 


“Value and Use System Records Long-Range Planning,” Kincaid, 
Journal, A.W.W.A., Vol. 52, February, 1960, pp. 277-283, 
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DRAG AND LIFT SPHERES WITHIN CYLINDRICAL 


Closure Donald Young 


DONALD analysis the results this investi- 
gation terms the conventional critical indeed interesting. 
Although, noted, his analysis speculative, does demonstrate how lift 
and drag data the type obtained may utilized sediment transport prob- 
lems. The writer agrees that for the relatively large values d/D used 
this investigation the contraction the flow certainly significant, and that 
the study should extended values d/D that would more char- 
acteristic sediment transport problems. pointed out the paper, the 
same technique was used inthis investigation can used extend the range 
values the parameters. However, the particle size becomes smaller 
will become increasingly difficult determine the conditions impending 
motion. 

With regard Laursen’s comment concerning the location the lift should 
noted that far the equilibrium equations for the particle are concerned 
the resultant fluid force acting the particle can resolved into components 
any point along the line actionof the force. The writer has arbitrarily re- 
solved the resultant fluid force into lift that passes through the center the 
sphere and drag that does not. would, course, incorrect arbitrar- 
ily assume that both components pass through the center the sphere. Re- 
gardless where the line actionof the lift located, the condition still re- 
mains that impending motion, direction normal the pipe wall, the re- 
sultant fluid force must pass through the sphere order for ro- 
tation not occur. 

The angle referred Laursen can readily obtained from and 
and the values and Cy, given Table For impending motion 
and therefore 


The remarks Indri concernthe difference flow conditions that exist be- 
tween the solid bodies falling ina still fluid and the bodies held 
equilibrium moving fluid. This distinction that frequently over- 
looked. Certainly the velocity profile the stream precedingthe solid has 
effect the force exerted the solid. Also, any turbulence the approach 
stream possibly significant. This problem that has concerned investi- 
gators working with wind tunnel models for some time. 


June 1960, Donald Young Paper 2526), 
Assoc. Prof., Dept. Theoretical and Applied Mechanics, Iowa State 
Science and Tech., Ames, 
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For the investigation under discussion the experimental conditions utilized 
were similar those normally encountered studies related the incipient 
motion suspensions, that is, the particles were rest and the fluid moving. 
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MODELS PRIMARILY DEPENDENT THE REYNOLDS 


Closure Simmons, Jr. 


SIMMONS, JR. clarification concepts and nomenclature pre- 
sented Priest, Kolf, and Reitmeyer are important those working with 
closed conduit flow problems. Priest correctly points out that Reynolds num- 
ber applies viscous motion and not directly turbulent motion. Kolf and 
Reitmeyer likewise point out distinction andaffirm the useof Euler relation- 
ships for analyzing the model data. Chapple’s contribution related and 
interesting problem well illustrates facet Reynolds number significance. 

The paper emphasized that through understanding the Reynolds num- 
ber, models large closed-conduit systems, which viscosity might thought 
predominant factor, may made relatively insensitive viscosity. 
This predicated the model being operated Reynolds numbers greater 
than some critical value below which viscous effects are significant. Such 
models are indeed analyzed with parameters the Euler variety that are di- 
mensionally similar but conceptually different from the Froude ones. Thus, 
the statement that data may analyzed using the same relations for 
models based the Froude law” misleading. should amended read 
data may analyzed using Euler relations.” pointed out, this 
conceptual difference involves change model techniques data analysis. 
point deserving further attention that determining Reynolds numbers 
high enough viscous forces can disregarded. definition, involves 
characteristic length and flow velocity the fluid system. These will have 
different significance different systems. For example, for orifice flow ina 
full pipeline based orifice diameter and flow velocity about 105 
appears sufficient.3 Similarly draft tube study based inlet con- 
ditions, about 105 was doubt many other instances 
satisfactory test results R’s less that can cited. However, 
minimum near minimum values are used, these critical values must 
determined through experiment. Unless worthwhileto make such tests, 
seems convenient and work inthe range the 106 value suggested 
the paper. 


June 1960, Simmons, Jr. Paper 2531). 

Engr., Div. Engrg. Lab., Reclamation, Denver, Colo. 

Model Studies Martin Dam Draft Tubes,” Kindsvater and 
Randolph, Jr., Transactions, ASCE, Vol, 120, 1955. 
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WATER EDDY FORCES OSCILLATING CYLINDERS* 


Discussion Donald VanSickle 


DONALD A.M. ASCE.—This paper most welcome addi- 
tion the authors’ series papers concerning wave forces The 
accumulation such data they have presented essential before design 
criteria can established that correctly account for the many factors involved. 

The writer particularly interested inthe relationship forces 
piles with tandem spacings because, another investigation, studied 
this effect for H-section piles uniform flow. interest tocompare the 
authors’ test results oscillating cylinders, their previous results cylin- 
ders with uniform the writer’s results H-piles uniform flow, 
and the results NACA investigation cylinders uniform 

facilitate this comparison, the various test results are summarized 
Fig. which the data are presented ratios the drag onthe downstream 
pile tandem pair tothe drag for single pile for various center-to-center 
pile spacings. was necessary ratios forthe NACA tests using 
the given drag coefficients and drag coefficient 1.2 for cylinder. 
For the authors’ tests oncylinders withuniform motion, the ratios were deter- 
mined from the mean values the drag the various velocities. 
The results for the H-piles were already the formof drag coefficients. The 
authors’ results oscillating cylinders are presented the ratios drag 
forces; however, the depths and velocities for the tandem cylinders corre- 
spond those for the single cylinder, the ratios drag forces should corre- 
spond the ratios drag coefficients. should mentioned that the center- 
to-center spacings for the square H-piles are pile widths, whereas those for 
the cylinders are cylinder diameters. 

examination Fig. reveals the following points: 


All the essentially the same shape. The ratioof drag co- 
efficients increases with increasedspacing about four diameters spacing. 
For spacings excess four diameters, the ratio drag coefficients re- 
mains practically constant. would appear, then, that center-to-center spa- 
cing approximately four diameters the critical spacing far drag 
forces are concerned, 


November 1960, Alan Laird, Charles Johnson, and Robert Walker 
Paper 2652), 

Hydr, Engr., Turner Collie Engrs., Inc., Houston, Tex. 

Walker, Proceedings, ASCE, Vol. 85, No. March, 1959. 

Experimental Investigation the Drag Characteristics Model H-Section 
Bridge Pier,” Donald VanSickle, thesis presented the Univ. Texas Austin, 
1958 partial fulfillment the requirements for the degree Master Science. 

“The Interference Between Struts Various Combinations,” David Biermann 
and Herrnstein, Jr., NACA Report No, 468, 1933. 
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The results the authors’ previous tests tandemcylinders with uni- 
form motion are substantial agreement with the results the NACA tests. 
Biermann and suggest that the negative drag forces 
close spacings are caused “the vortices produced the front cylinder 
(which) probably encircle the rear cylinder, impinging the back surface 
with sufficient force produce forward reaction.” The directions the 
eddies shown indicate that this certainly reasonable explanation. 
Both sets test results further indicate, for spacings greater than four di- 
ameters, leveling off the drag coefficient for cylinder with neighbor 
from 0.25 0.40 the drag coefficient for single cylinder. Biermann and 
advance the hypothesis that “the boundary layer the second 
cylinder madeturbulent the rapid changes velocity and direction the 
flow behind the first cylinder.” The well-known fact that the drag coefficient 


Oscillating cylinders, 
by authors 


by Van Sickle ® 


Cylinders with uniform motion, 


0.4 authors® 
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Ratio of 


Center-to-center spacing, pile diameters 


FIG. DRAG COEFFICIENT RATIO WITH PILE SPACING 


cylinder decreases approximately one-fourth its normal value when 
the boundary layer becomes turbulent!1 lends much support tothis hypothesis. 

The results the study tandem H-piles are shown shaded area 
Fig. The lower line the area represents Froude number 0.2, and 
the upper line represents Froude number 0.6. avoid complicating Fig. 
the curves for intermediate values Froude number are not shown; how- 
ever, the Froude number increases from 0.2 0.6, the curves shift upward 
fairly uniform steps. The shape these curves very similar the shape 
the curves uniform flow. The negative drag coefficients for 
close spacings can accounted for suggested point Although, for spa- 
cings greater thanfour diameters, similar leveling off the drag coefficient 


Hunter Rouse, ed, Engineering Hydraulics, John Wiley Sons, New York, 1950, 
124, 


Cylinders uniform flow, 
| V=64 | 
: | | 
| 


DISCUSSION 203 


for H-pile with neighbor occurs values 0.60 0.80 the coefficient for 
single H-pile, convenient explanation available. The point separation 
for the H-pile essentially fixed the sharp flange edge, and shift down- 
stream, with reduction drag coefficient, occurs when the flow becomes 
turbulent. 

contrast these observations for tandem piles uniform flow, the 
authors’ data for oscillating cylinders indicate that oc- 
Furthermore, while the drag coefficient levels off for spacings greater 
than four diameters, does the same value for single cylinder. 
other words, the authors point out, for oscillating cylinders the effect the 
upstream cylinder negligible for spacings greater than four diameters. 


occurred the writer that because the authors’ results were presented 
ratio the force the cylinder with neighbor the force ona single 
cylinder, the leveling off the ratio value 1.0 could explained the 
oscillatory motion, for some unknown reason, caused the boundary layer for 
the single cylinder turbulent. The occurrence turbulent flow the 
downstream cylinder the tandem pair would then produce the turbulent drag 
force, and drag ratio 1.0. The authors point out, however, that “the drag 
coefficients corrected infinite cylinder values were the usual 1.3 1.4 for 
the vertical cylinders extending through the water surface.” Because this 
the drag coefficient for laminar boundary layer conditions, evident that 
the previously described explanation cannot apply. 

Another alternative that the downstream cylinder oscillating tan- 
dem pair does not encounter turbulent whereas the downstream 
cylinder uniform flow does. There would seem noobvious reason why 
such condition should exist. 

would not seem factor because the authors’ pre- 
vious work, using essentially the same equipment and procedures, gave results 
substantial agreement with the NACA tests. 

There remains the further possibility that the hypothesis advanced Bier- 
mann and Herrnstein explain the reduced drag the downstream cylinder 
not valid. The fact that their hypothesis cannot explain the reduced drag 
the downstream H-pile adds weight this argument. Such conclusion, how- 
ever, does not explain why piles uniform flow have reduced drag the down- 
stream pile spacings excess four diameters, whereas cylinders os- 
cillating motion not have reduced drag for these spacings. 

the authors know any reason, could they present reasonable hy- 
pothesis, why this situation should exist? 

Could they state whether not they actually encountered negative drag 
forces any their tests oscillating cylinders? 
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FISH PASSAGE THROUGH HYDRAULIC TURBINES* 


ASCE.—The research the passage young fish through 
hydraulic turbines, Von Gunten and his associates, one the most im- 
portant projects arising from the 3.5 million dollar fisheries-engineering re- 
search program initiated the Corps United States Army Engineers. 

These recent tests furnish further evidence that the Corps Engineers’ 
program, participated all the major fisheries organizations the Pacific 
Northwestern states, has paid handsome dividends. wealth data has been 
published giving criteria for the safe passage fish past dams.3,4,5 The 
writer believes that this program and further research will very likely show 
that rivers can developed for multi-purposes, including the preservation 
anadromous fish. 

The preliminary tests the passage fingerlings through small Francis 
and Kaplan runners diameter inthe laboratory showed that cavitating 
conditions turbine will kill fish. The rate mortality related the 
magnitude partial vacuum present the water passages. 

The full scaletests Cushman dam confirmed laboratory findings showing 
that significant mortality does not occur among young fish subjected ab- 
rupt reductionin water pressure from psi atmospheric. The tests 
showed significant improvement survival rates resulting from the reduc- 
tion cavitation the turbine whenthe tailwater level rose ft. The maxi- 
mum mortality occurred when the turbine operated under partial loading, 
condition known result worsening cavitation phenomena. 

The writer had the opportunity witnessing the Cushman tests and was par- 
ticularly interested observing the fish for evidence the bends. Experi- 
ments the University British Columbia (Vancouver, C.) have shown 
that the bends are likely occur among fish subjected partial vacuum 
high intensity for periods time considerably greater than the time pas- 
sage through the Cushman turbines. The writer saw the bends 
among the Cushman fish after the tests and believes that few fish were killed 
seriously injured the effects vacuum unaccompanied cavi- 


May 1960, Glenn Von Gunten Paper 2812), 
Head, Dept. Engrg., Univ. British Columbia, Vancouver, Canada. 
“Progress Report Fisheries-Engineering Research Program,” Corps Engrs., 
Army, Pacific Div., November, 
Report Fisheries-Engineering Research Program,” Corps Engrs., 
Army, Pacific Div., July, 1960, 

Annotated Bibliography Fish Passing Facilities, Committee Power Project 
Planning and Design, Power Div., ASCE, September, 1960, 
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One interesting feature the Cushman tests was the efficiency the large 
net used totrap fingerlings. The mouth the net completely covered the draft 
tube exit and trapped approximately 90% the fish passing through the turbine 
during the tests. Ivan Donaldson and his colleagues are commended 
for their pioneering work the design large nets for this purpose. 

The research Von Gunten and his associates points further program 
designed solve the problem successful passage fish through turbines. 
One important aspect this problem the relationships be- 
tween turbine specific speed, the cavitation parameter sigma, and fish mor- 
tality. The Cushman turbines are placed elevation about above tail- 
water level high tide. future plant similar Cushman, sigma could 
increased more than 50% lowering the turbines about ft. One wonders 
about the extent improvement fish survival rates with such substantial 
increase the value sigma. 

the light the author’s findings, would interest would dis- 
cuss the changes operating conditions called for existing plants the 
Columbia River system. During the migration period, nowstandard prac- 
tice avoid turbine operation under low inplants which young 
fish pass through turbines their way the sea? 

wouldalso interest tohave Von Gunten’s comments the possibility 
safe passage migrants through the Grand Coulee and Chief Joseph tur- 
bines. would now appear that the rehabilitation the salmonruns the up- 
per Columbia River within the realm future possibilities. 
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“Uniform Water Conveyance Channels Alluvial Material” Daryl 
Simons and Maurice Albertson, May 1960, Proc. Paper 2484. 


64. Under the heading “Design Procedures: Estimating Knowing 
Vice Versa” lines and change the word wetter wetted. 
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